Dynamics of the mid-ocean ridge plate boundary: Recent observations and theory

The global mid-ocean ridge system is one of the most active plate boundaries on the earth and understanding the dynamic processes at this plate boundary is one of the most important problems in geodynamics. In this paper I present recent results of several aspects of mid-ocean ridge studies concerning the dynamics of oceanic lithosphere at these diverging plate boundaries. I show that the observed rift valley to no-rift valley transition (globally due to the increase of spreading rate or locally due to the crustal thickness variations and/or thermal anomalies) can be explained by the strong temperature dependence of the power law rheology of the oceanic lithosphere, and most importantly, by the difference in the rheological behavior of the oceanic crust from the underlying mantle. The effect of this weaker lower crust on ridge dynamics is mainly influenced by spreading rate and crustal thickness variations. The accumulated strain pattern from a recently developed lens model, based on recent seismic observations, was proposed as an appealing mechanism for the observed gabbro layering sequence in the Oman Ophiolite. It is now known that the mid-ocean ridges at all spreading rates are offset into individual spreading segments by both transform and nontransform discontinuities. The tectonics of ridge segmentation are also spreading-rate dependent: the slow-spreading Mid-Atlantic Ridge is characterized by distinct bulls-eye shaped gravity lows, suggesting large along-axis variations in melt production and crustal thickness, whereas the fast-spreading East-Pacific Rise is associated with much smaller along-axis variations. These spreading-rate dependent changes have been attributed to a fundamental differences in ridge segmentation mechanisms and mantle upwelling at mid-ocean ridges: the mantle upwelling may be intrinsically plume-like (3-D) beneath a slow-spreading ridge but more sheet-like (2-D) beneath a fast-spreading ridge.     

The Mid-Atlantic Ridge between 29°N and 31°30′N in the last 10 Ma

The segmentation of the Mid-Atlantic Ridge between 29°N and 31°30′ N during the last 10 Ma was studied. Within our survey area the spreading center is segmented at a scale of 25–100 km by non-transform discontinuities and by the 70 km offset Atlantis Transform. The morphology of the spreading center differs north and south of the Atlantis Transform. The spreading axis between 30°30′N and 31°30′N consists of enéchelon volcanic ridges, located within a rift valley with a regional trend of 040°. South of the transform, the spreading center is associated with a well-defined rift valley trending 015°. Magnetic anomalies and the bathymetric traces left by non-transform discontinuities on the flanks of the Mid-Atlantic Ridge provide a record of the evolution of this slow-spreading center over the last 10 Ma. Migration of non-transform offsets was predominantly to the south, except perhaps in the last 2 Ma. The discontinuity traces and the pattern of crustal thickness variations calculated from gravity data suggest that focused mantle upwelling has been maintained for at least 10 Ma south of 30°30′ N. In contrast, north of 30°30′N, the present segmentation configuration and the mantle upwelling centers inferred from gravity data appear to have been established more recently. The orientation of the bathymetric traces suggests that the migration of non-transform offsets is not controlled by the motion of the ridge axis with respect to the mantle. The evolution of the spreading center and the pattern of segmentation is influenced by relative plate motion changes, and by local processes, perhaps related to the amount of melt delivered to spreading segments. Relative plate motion changes over the last 10 Ma in our survey area have included a decrease in spreading rate from 32 mm a⁻¹ to 24 mm a⁻¹, as well as a clockwise change in spreading direction of 13° between anomalies 5 and 4, followed by a counterclockwise change of 4° between anomaly 4 and the present. Interpretation of magnetic anomalies indicates that there are significant variations in spreading asymmetry and rate within and between segments for a given anomaly time. These differences, as well as variations in crustal thickness inferred from gravity data on the flanks of spreading segments, indicate that magmatic and tectonic activity are, in general, not coordinated between adjacent spreading segments.     

Statistical analysis of isotopic ratios in MORB: the mantle blob cluster model and the convective regime of the mantle

A statistical examination of isotopic distributions for MORB from various ocean ridges leads to the “blob cluster model”, in which the oceanic crust accreting at ridges results from the mixing of two components within the ascending mantle. These are (1) upper mantle material and (2) discrete rising blobs of more radiogenic material. The blobs are fractionated to a variable degree and are distributed in the upper mantle circulation in a manner that is related to the spreading rate.(1) Themean values of the isotopic distributions allow us to calculate the probabilities of the two types of material within the mantle. The results show that theproportion of asthenospheric material in the mixtureincreases with the spreading rate, in agreement with the hypothesis of blob dilution within the upper mantle convection.Mass fluxes can be estimated for the rising blobs from these probabilities, which depend on the respective concentrations in the sources of the two types of material. If the blobs originate in the lower mantle, this flux estimation would suggest that a significant part of the lower mantle has been injected into the upper mantle during earth history.(2) Thestandard deviations of the distributions depend on the “efficiency” of the mixing process:the more imbricated are the asthenospheric and blob materials in the mixture,the smaller is theisotopic spread. This efficiency parameter is shown to increase with the spreading rate, as already suggested by previous comparisons between the East Pacific Rise and the Mid-Atlantic Ridge. Moreover, this feature may also be correlated with other data such as ridge bathymetric variations.     

The mid-ocean ridge axial valley as a steady-state neck

In this paper the mid-ocean ridge axial valley is modelled as a steady-state lithospheric neck in which lithospheric stretching balances lithospheric accretion. Conversely, the axial high is a steady-state lithospheric bulge. The lithosphere is modelled as a thin plate with a Newtonian rheology. It is shown that an axial valley will occur if the rate of viscosity increase away from the ridge axis is faster than the rate at which accretion decreases. An axial high will occur if the opposite condition holds. This is consistent with the observation that axial valleys occur at low spreading rates and axial highs at high spreading rates. By fitting our model to profiles across the Mid-Atlantic Ridge and the East Pacific Rise and assuming the lithospheric thickness at the ridge axis to be 5 km, we find accretion widths of 6–8 km. We find the width over which there is a significant increase in lithospheric viscosity to be also 6–8 km.     

Asymmetric mantle dynamics in the MELT region of the East Pacific Rise

The mantle electromagnetic and tomography (MELT) experiment found a surprising degree of asymmetry in the mantle beneath the fast-spreading, southern East Pacific Rise (MELT Seismic Team, Science 280 (1998) 1215–1218; Forsyth et al., Science 280 (1998) 1235–1238; Toomey et al., Science 280 (1998) 1224–1227; Wolfe and Solomon, Science 280 (1998) 1230–1232; Scheirer et al., Science 280 (1998) 1221–1224; Evans et al., Science 286 (1999) 752–756). Pressure-release melting of the upwelling mantle produces magma that migrates to the surface to form a layer of new crust at the spreading center about 6 km thick (Canales et al., Science 280 (1998) 1218–1221). Seismic and electromagnetic measurements demonstrated that the distribution of this melt in the mantle is asymmetric (Forsyth et al., Science 280 (1998) 1235–1238; Toomey et al., Science 280 (1998) 1224–1227; Evans et al., Science 286 (1999) 752–756) at depths of several tens of kilometers, melt is more abundant beneath the Pacific plate to the west of the axis than beneath the Nazca plate to the east. MELT investigators attributed the asymmetry in melt and geophysical properties to several possible factors: asymmetric flow passively driven by coupling to the faster moving Pacific plate; interactions between the spreading center and hotspots of the south Pacific; an off-axis center of dynamic upwelling; and/or anomalous melting of an embedded compositional heterogeneity (MELT Seismic Team, Science 280 (1998) 1215–1218; Forsyth et al., Science 280 (1998) 1235–1238; Toomey et al., Science 280 (1998) 1224–1227; Wolfe and Solomon, Science 280 (1998) 1230–1232; Evans et al., Science 286 (1999) 752–756). Here we demonstrate that passive flow driven by asymmetric plate motion alone is not a sufficient explanation of the anomalies. Asthenospheric flow from hotspots in the Pacific superswell region back to the migrating ridge axis in conjunction with the asymmetric plate motion can create many of the observed anomalies.     

Inversion of gravity anomalies over spreading oceanic ridges

Models of spreading ocean ridges are derived by Bayesian gravity inversion with geophysical and geodynamic a priori information. The aim is to investigate the influence of spreading rate, plate dynamics and tectonic framework on crust and upper mantle structure by comparing the Mid Atlantic Ridge (MAR), the Indian Ocean Ridge (IND) and the East Pacific Rise (PAC). They differ in mean spreading rate, dynamic settings, as attached slabs, and plume interaction. Topography or bathymetry, gravity, isostasy, seismology and geology, etc. are averaged along the ridges and guide the construction of initial 2D models, including features as mean plumes, i.e. averaged along the ridge. This is a gross simplification, and the results are considered preliminary.Three model types are tested: (a) the temperature anomaly; (b) asthenospheric rise into thickening lithosphere; (c) a crustal root as had been anticipated before seafloor spreading was discovered. Additional model components are a mean plume, a non-compensated ridge uplift, an under-compensated asthenospheric rise, e.g. of partially molten material, and seismic velocity models for P and S waves. Model type (c), tends to permute to model type (b) from thick crust to thin axial lithosphere. Model type (a) renders ‘realistic’ values of the thermal expansivity, but is insufficient to fit the gravity data; partial melt may disturb the simple temperature effect. A combination of (a) and (b) is most adequate. Exclusive seismic velocity models of S or P waves do not lead to acceptable densities nor to adequate gravity fitting. The different ridges exhibit significant differences in the best models: ATL and IND show an axial mass excess fostering enhanced ridge push, and ATL, in addition, suggests a mean plume input, while PAC shows an axial mass deficit reducing ridge push, most probably due to dominance of slab pull in the force balance.Goodness of the gravity fit alone is no justifiable criterion for goodness of model, indeed minor modifications to each model within the uncertainties of the assumptions can make the fit arbitrarily good. Goodness of model is quantified exclusively by a priori information.     

Subduction influence on magma supply at the East Scotia Ridge

Despite a spreading rate of 65–70 km Ma⁻¹, the East Scotia Ridge has, along most of its length, a form typically associated with slower rates of sea floor spreading. This may be a consequence of cooler than normal mantle upwelling, which could be a feature of back-arc spreading. At the northern end of the ridge, recently acquired sonar data show a complex, rapidly evolving pattern of extension within 100 km of the South Sandwich Trench. New ridge segments appear to be nucleating at or near the boundary between the South American and Scotia Sea plates and propagating southwards, supplanting older segments. The most prominent of these, north of 56°30′S, has been propagating at a rate of approximately 60 km Ma⁻¹ for at least 1 Ma, and displays a morphology unique on this plate boundary. A 40 km long axial high exists at the centre of this segment, forming one of the shallowest sections of the East Scotia Ridge. Beneath it, seismic reflection profiles reveal an axial magma chamber, or AMC, reflector, similar to those observed beneath the East Pacific Rise and Valu Fa Ridge. Simple calculations indicate the existence here of a narrow (<1 km wide) body of melt at a depth of approximately 3 km beneath the sea floor. From the topographic and seismic data, we deduce that a localised mantle melting anomaly lies beneath this segment. Rates of spreading in the east Scotia Sea show little variation along axis. Hence, the changes in melt supply are related to the unique tectonic setting, in which the South American plate is tearing to the east, perhaps allowing mantle flow around the end of the subducting slab. Volatiles released from the torn plate edge and entrained in the flow are a potential cause of the anomalous melting observed. A southward mantle flow may have existed beneath the axis of the East Scotia Ridge throughout its history.     

Topography and evolution of the East Pacific Rise between 5°S and 20°S

In the absence of convincing magnetic anomaly information, topographic profiles have been used to infer the tectonic history of the East Pacific Rise between 5°S and 20°S. Profiles projected at right angles to the rise crest show a sharp drop in elevation at roughly the same distance on either side of the crest. Profiles to the east of the rise also show a second topographic high at 95°W. Comparison of these profiles with empirical depth versus age curves for the North Pacific suggests that this rise, the Galapagos Rise, is the fossil East Pacific Rise which terminated close to 6 mybp by the spreading center jumping 900 km to the west. The extreme youthfulness of the present East Pacific Rise, the step structure of its flanks, and the similarity in age of the top of this step and the crest of the Galapagos Rise substantiate this interpretation. This jump coincided with a similar readjustment involving the Mathematicians ridge at 5° to 20°N and the opening of the Gulf of California.     

A nonstationary model of the thermal regime of axial zones of mid-ocean ridges: Formation of crustal and mantle magma chambers

A nonstationary model of spreading with periodic intrusions of a molten material into an axial zone of a mid-ocean ridge (MOR) is applied to numerical analysis of the thermal state in MOR axial zones and the formation of crustal and mantle magma chambers in them. The model satisfactorily explains the positions, dimensions, and shapes of magma chambers, as well as variations in these parameters depending on the spreading rate, temperature, and composition of crustal and mantle rocks. The release and absorption of the latent heat of rock melting, hydrothermal heating of the crust, and variations in the solidus and liquidus temperatures of crustal and mantle rocks as a function of their composition are factors controlling the shape and position of crustal magma chambers.     

The Juan Fernandez microplate north of the Pacific-Nazca-Antarctic plate junction at 35°S

Results of the R/V “Thomas Washington” Pascua 3 expedition provide evidence for the existence of the Juan Fernandez microplate just north of the junction between the East Pacific Rise (EPR) and the Chile fracture zone. Prior to Pascua 3, the microplate in the region had been hypothesized from the pattern of seismicity. The eastern and western boundaries of the Juan Fernandez microplate are well defined and represent north-south trending spreading centers characterized by very slow and very fast rates of accretion respectively. In agreement with the rates, the eastern boundary is represented by a rift valley and the western boundary by an EPR-type axial ridge. The northern boundary of the Juan Fernandez microplate is a 100°-trending wide fracture zone complex which may have resulted from northward transform fault migration. The fracture zone fails to meet the zone of accretion at the Pacific-Nazca-Juan Fernandez triple junction. In this area the zone of accretion displays a double ridge with a large overlap. The southern boundary of the Juan Fernandez microplate is still poorly constrained. The plate geometry derived from SEABEAM differs from that derived by Anderson-Fontana et al. (1986) [14] from a plate motion inversion scheme using primarily earthquake first-motion solutions together with limited bathymetric and magnetic data.     

Back-arc basin basalt systematics

The Mariana, east Scotia, Lau, and Manus back-arc basins (BABs) have spreading rates that vary from slow (<50 mm/yr) to fast (>100 mm/yr) and extension axes located from 10 to 400 km behind their island arcs. Axial lava compositions from these BABs indicate melting of mid-ocean ridge basalt (MORB)-like sources in proportion to the amount added of previously depleted, water-rich, arc-like components. The arc-like end-members are characterized by low Na, Ti and Fe, and by high H₂O and Ba/La; the MORB-like end-members have the opposite traits. Comparisons between basins show that the least hydrous compositions follow global MORB systematics and an inverse correlation between Na8 and Fe8. This is interpreted as a positive correlation between the average degree and pressure of mantle melting that reflects regional variations in mantle potential temperatures (Lau/Manus hotter than Mariana/Scotia). This interpretation accords with numerical model predictions that faster subduction-induced advection will maintain a hotter mantle wedge. The primary compositional trends within each BAB (a positive correlation between Fe8, Na8 and Ti8, and their inverse correlation with H₂O(8) and Ba/La) are controlled by variations in water content, melt extraction, and enrichments imposed by slab and mantle wedge processes. Systematic axial depth (as a proxy for crustal production) variations with distance from the island arc indicate that compositional controls on melting dominate over spreading rate. Hydrous fluxing enhances decompression melting, allowing depleted mantle sources just behind the island arc to melt extensively, producing shallow spreading axes. Flow of enriched mantle components around the ends of slabs may augment this process in transform-bounded back-arcs such as the east Scotia Basin. The re-circulation (by mantle wedge corner flow) to the spreading axes of mantle previously depleted by both arc and spreading melt extraction can explain the greater depths and thinner crust of the East Lau Spreading Center, Manus Southern Rifts, and Mariana Trough and the very depleted lavas of east Scotia segments E8/E9. The crust becomes mid-ocean ridge (MOR)-like where the spreading axes, further away from the island arc and subducted slab, entrain dominantly fertile mantle.     

Sr and Nd isotopes in basalts from the East Pacific Rise: significance for mantle heterogeneity

Isotopic data for Sr and Nd from fresh glassy East Pacific Rise basalts suggest that this part of the suboceanic mantle is characterized by subtle but distinct large-scale regional isotopic variability which may reflect differences between cells of the convecting mantle. In spite of a systematic N—S change in spreading rate of a factor of three along the sampled portion of the EPR, no correlation is observed between spreading rate and range of isotopic composition, indicating that the regional variations override homogenization effects which may be correlated with rate of magma generation and hence spreading rate. There is no clear signature in our data of effects from the postulated global “Dupal Anomaly” [30,31]. However, for a restricted ridge segment at the latitude of Easter Island, anomalously high⁸⁷Sr/⁸⁶Sr and low¹⁴³Nd/¹⁴⁴Nd occur, coupled with high incompatible element concentrations. These features are most easily understood as being the result of inclusion of a “plume” component in these ridge basalts.     

Isotope and trace element geochemistry of young Pacific seamounts: implications for the scale of upper mantle heterogeneity

Basalts from young seamounts situated within 6.8 m.y. of the East Pacific Rise, between 9° and 14°N latitude, display significant variations in ¹⁴³Nd/¹⁴⁴Nd (0.51295–0.51321), ⁸⁷Sr/⁸⁶Sr (0.7025–0.7031), and(La/Sm)_N (0.415–3.270). Nd and Sr isotope ratios are anti-correlated and form a trend roughly parallel to the “mantle array” on a¹⁴³Nd/¹⁴⁴Nd vs.⁸⁷Sr/⁸⁶Sr variation diagram. Nd and Sr isotope ratios display negative and positive correlations, respectively, with(La/Sm)_N. The geochemical variations observed at the seamounts are nearly as great or greater than those observed over several hundred kilometers of the Reykjanes Ridge, or at the islands of Iceland or Hawaii.

Samples from one particular seamount, Seamount 6, display nearly the entire observed range of chemical variations, offering an ideal opportunity to constrain the nature of heterogeneities in the source mantle. Systematics indicative of magma mixing are recognized when major elements, trace elements, trace element ratios, and isotope ratios are compared with each other in all possible permutations. The source materials required to produce the end-member magmas are: (1) a typical MORB-source-depleted peridotite; and (2) a relatively enriched material which may represent ancient mantle segregations of basaltic melt, incompletely mixed remnants of subducted ocean crust, or metasomatized peridotite such as that found at St. Paul's Rocks or Zabargad Island. Due to the proximity of the seamounts to the East Pacific Rise (EPR), the source materials are thought to comprise an intimate mixture in the mantle immediately underlying the seamounts and the adjacent EPR. Lavas erupted at the ridge axis display a small range of isotopic and incompatible trace element compositions because the large degrees of melting and presence of magma chambers tend to average the chemical characteristics of large volumes of mantle.

If the postulated mantle materials, with large magnitude, small-scale heterogeneities, are ubiquitous in the upper mantle, chemical variations in basalts ranging from MOR tholeiites to island alkali basalts may reflect sampling differences rather than changes in bulk mantle chemistry.     

Linear volcanism in French Polynesia

The island chains of French Polynesia form subparallel line segments whose southeasterly extensions are perpendicular to the East Pacific Rise, the site of present sea-floor spreading in the eastern Pacific Ocean. Samples collected from island members of the Society and Austral Islands chains are used, together with previously reported age determinations for the Marquesas and Pitcairn-Gambier Islands, in a geochronological study of the southeastward migration of volcanism in each of those four lineaments. The suggestion from geomorphologic evidence that island ages increase to the northwest within each island chain, is confirmed by K---Ar whole-rock ages. The linear volcanism which built the islands of French Polynesia began in the Miocene and continues today.Rates of migration of volcanism are calculated from the nearly linear relationship between average island ages and distance from the southeast ends of the four island lineaments. The four rates are indistinguishable, within limits of detection, at 11 ± 1 cm/year. These rates are consistent with the model of rigid Pacific plate movement over four fixed sources of volcanism, be they dynamic as in the hot spot/plume models or passive as in models of propagating lithospheric fractures. If it is accepted that these volcanic sources trace the motion of the lithosphere over the mantle and thus define the “absolute” frame of reference for plate movement, Pacific plate motion may be fixed to the geometry and volcanic migration rates of French Polynesia. This allows calculation of the absolute motion of all other plates, providing an accurate relative motion model is known (Minster et al., 1974). Such a calculation predicts that Africa is virtually stationary and that the Mid-Atlantic Ridge and East Pacific Rise are moving slowly to the west.     

Fundamentals of ridge crest topography

A linear relationship between the sea floor depth and the square root of age has been found for ocean lithosphere spreading from mid-ocean ridges. The asymptotic solution of depth as a function of age for the thermally contracting lithosphere predicts a linear dependence of depth ontwith a proportionality involving the initial lithosphere temperature, the thermal diffusivity, and the isostatic expansion coefficient averaged to include any temperature dependent phase changes. Empirical depth observations, when plotted as a function of the square root of age, bear out this prediction well, but there is a variation in the gradient,ht, along the ridge on a fine scale (up to 20% over 200 km). This implies a fundamental variation of the contraction parameter over the same scale, most probably of compositional origin. Details of a more complete cooling model near the ridge crest, including a crust of different thermal parameters than those of the mantle, predict a crestal height about 0.2 km below that of the simplified model. Individual profiles from the southeast Pacific show no such crestal deviation, and it is concluded that by quickly cooling the new crust, hydrothermal circulation may remove any effects of the crust which would be seen in the topography of a lithosphere cooled totally by conduction. The straightness of depth versust for older ocean data (to 80 m.y.) precludes any basal isothermal boundary shallower than 100 km.     

Paleobathymetry of the crest of spreading ridges related to the age of ocean basins

Measurements of the crest of the spreading ridge in the young ocean basins of the Afar region and Gulf of Aden and in the mature Indian, Atlantic, and Pacific Oceans show that the depth of the ridge crest is correlated (r = 0.99) with the logarithm of the age of the ocean basin. Ridge crests in a very young basin (Afar) are at sea level, at about 1.5 km in young basins (Gulf of Aden), and at about 2.6 km in mature basins (Indian, Atlantic, Pacific). A new curve that relates crestal depth and age of the ocean basin is coupled with the existing depth/age curve for oceanic crust in a comprehensive scheme which can be used for relating depth and age of oceanic crust.     

Mantle flow and melting beneath oceanic ridge–ridge–ridge triple junctions

Plate boundary geometry likely has an important influence on crustal production at mid-ocean ridges. Many studies have explored the effects of geometrical features such as transform offsets and oblique ridge segments on mantle flow and melting. This study investigates how triple junction (TJ) geometry may influence mantle dynamics. An earlier study [Georgen, J.E., Lin, J., 2002. Three-dimensional passive flow and temperature structure beneath oceanic ridge-ridge-ridge triple junctions. Earth Planet. Sci. Lett. 204, 115–132.] suggested that the effects of a ridge–ridge–ridge configuration are most pronounced under the branch with the slowest spreading rate. Thus, we create a three-dimensional, finite element, variable viscosity model that focuses on the slowest-diverging ridge of a triple junction with geometry similar to the Rodrigues TJ. This spreading axis may be considered to be analogous to the Southwest Indian Ridge. Within 100 km of the TJ, temperatures at depths within the partial melting zone and crustal thickness are predicted to increase by ~ 40 °C and 1 km, respectively. We also investigate the effects of differential motion of the TJ with respect to the underlying mantle, by imposing bottom model boundary conditions replicating (a) absolute plate motion and (b) a three-dimensional solution for plate-driven and density-driven asthenospheric flow in the African region. Neither of these basal boundary conditions significantly affects the model solutions, suggesting that the system is dominated by the divergence of the surface places. Finally, we explore how varying spreading rate magnitudes affects TJ geodynamics. When ridge divergence rates are all relatively slow (i.e., with plate kinematics similar to the Azores TJ), significant along-axis increases in mantle temperature and crustal thickness are calculated. At depths within the partial melting zone, temperatures are predicted to increase by ~ 150 °C, similar to the excess temperatures associated with mantle plumes. Likewise, crustal thickness is calculated to increase by approximately 6 km over the 200 km of ridge closest to the TJ. These results could imply that some component of the excess volcanism observed in geologic settings such as the Terceira Rift may be attributed to the effects of TJ geometry, although the important influence of features like nearby hotspots (e.g., the Azores hotspot) cannot be evaluated without additional numerical modeling.     

Petrogenesis of diabase from accretionary prism in the southern Qiangtang terrane,central Tibet: Evidence from U–Pb geochronology,petrochemistry and Sr–Nd–Hf–O isotope characteristics

The Bangong–Nujiang suture (BNS) between the Lhasa and Qiangtang terranes is an important boundary and its petrogenesis is controversial. Diabase from the accretionary prism in the southern Qiangtang terrane yields a zircon U–Pb age of 181.3 ± 1.4 Ma. All the diabases show tholeiitic basalt compositions, gentle enrichment patters of light rare earth elements (REE), variable enrichment in incompatible element concentrations (e.g. Th and Rb), and no anomaly in high field strength elements (e.g. Nb and Ta), similar to that of enriched mid‐ocean ridge basalt (E‐MORB). They have relatively homogeneous whole rock Nd (εNd(t) = 7.3–9.1) and zircon Hf–O isotopic compositions (εHf(t) = 14.8–16.1, and δ¹⁸O = 4.57–6.12‰), possibly indicating melting of the depleted mantle and no significant crustal contamination during the petrogenesis. The element variations suggest that the diabases were formed by plume–ridge interaction at a mid‐ocean ridge within the Bangong–Nujiang ocean.     

Thermal and rheological state of the lithosphere and specific features of structuring in the rift zone of the Reykjanes Ridge (from the results of numerical and experimental modeling)

Specific features of the bottom topography structure and the character of morphostructural segmentation of the rift zone of the Reykjanes Ridge change substantially along the ridge strike with increasing distance from Iceland’s hotspot. A clearly pronounced regularity of changes is observed in the rift zone’s morphology from the axial uplift (in the northern part of the ridge) to the rift valleys (in the southern part of the ridge) through an intermediate or transitional type of morphology. The results of numerical modeling showed that changes in the rift zone’s morphology along the Reykjanes Ridge strike are largely caused by changes in the degree of mantle heating and depend on the intensity of magma supply. It is shown that under conditions of ultraslow spreading, it is these parameters that control the presence or absence of crustal magma chambers, as well as the thickness of the effectively-elastic layer of the axial lithosphere. The experimental modeling of topography-forming deformations and structuring on the Reykjanes Ridge showed that under oblique extension, specific features of the formation of axial fractures and the character of their segmentation mainly depend on the thickness of the axial lithosphere, its heating zone width, and the kinematics of spreading. The experiments also showed that the tendency of fractures to develop obliquely to the extension axis is caused by the action of the inclined zone of the location of the deformation, and shear deformations play a substantial role in the lithosphere’s destruction as the inclination angle increases.     

Tracing the history of submarine hydrothermal inputs and the significance of hydrothermal hafnium for the seawater budget—a combined Pb-Hf-Nd isotope approach

Secular variations in the Pb isotopic composition of a mixed hydrogenous-hydrothermal ferromanganese crust from the Bauer Basin in the eastern Equatorial Pacific provide clear evidence for changes in hydrothermal contributions during the past 7 Myr. The nearby Galapagos Rise spreading center provided a strong hydrothermal flux prior to 6.5 Ma. After 6.5 Ma, the Pb became stepwise more radiogenic and more similar to Equatorial Pacific seawater, reflecting the westward shift of spreading to the presently active East Pacific Rise (EPR). A second, previously unrecognized enhanced hydrothermal period occurred between 4.4 and 2.9 Ma, which reflects either off-axis hydrothermal activity in the Bauer Basin or a late-stage pulse of hydrothermal Pb from the then active, but waning Galapagos Rise spreading center.Hafnium isotope time-series of the same mixed hydrogenous-hydrothermal crust show invariant values over the past 7 Myr. Hafnium isotope ratios, as well as Nd isotope ratios obtained for this crust, are identical to that of hydrogenous Equatorial Pacific deep water crusts and clearly indicate that hydrothermal Hf, similar to Nd, does not travel far from submarine vents. Therefore, we suggest that hydrothermal Hf fluxes do not contribute significantly to the global marine Hf budget.