The Salton Trough is one of the few regions on Earth where rifting is subaerial instead of submarine. We use the relative attenuation of teleseismic P phases recorded by the Salton Trough Seismic Imaging Project to investigate lithospheric and asthenospheric structures that form during extension. Map‐view analysis reveals stronger attenuation within the Salton Trough than in the adjacent provinces. We then construct tomographic models for variations in seismic attenuation with depth to discriminate between crustal and mantle signals with a damped least squares approach and a Bayesian approach. Synthetic tests show that models from damped least squares significantly underestimate the strength of attenuation and cannot separate crustal and mantle signals even if the tomographic models are allowed to be discontinuous at the lithosphere‐asthenosphere boundary (LAB). We show that a Bayesian approach overcomes these problems when inverting the same synthetic data sets and that shallow and deep signals are more clearly separated when imposing a discontinuity. With greater than 95% confidence, the results reveal first, that attenuation occurs primarily beneath the LAB; second, that the width of the attenuative region is narrower than the rift at 120 km depth; and third, that the strength of attenuation requires that the attenuative feature represents a melting‐column similar to those beneath mid‐ocean ridges. The narrow width of the melting column below the volatile‐free solidus is inconsistent with models for passive upwelling, where flow is driven only by rifting. Instead, we attribute the generation of incipient oceanic crust to mantle upwelling focused by buoyancy into a narrow diapir.
Geophysical observations at some mid-ocean ridges document an across-axis asymmetry in indicators of magma production. Other observations are interpreted as showing non-monotonic variations in the depth of the lithosphere–asthenosphere boundary. These patterns are inconsistent with the classical models of mantle corner flow and half-space cooling. To investigate this discrepancy, we use models of coupled magma/mantle dynamics beneath mid-ocean ridges in which phase densities are determined by melt–residue partitioning of iron and magnesium, and bulk density is affected by residual porosity. Our models predict that emergent gradients in density drive ridge-local convection. In particular, we show that convective upwelling is enhanced by porous buoyancy and suppressed by compositional buoyancy. Despite this suppression, models that include both compositional and porous buoyancy are more sensitive to long-wavelength mantle heterogeneity than models with porous buoyancy alone. This sensitivity enables models to readily form across-axis asymmetry of upwelling. In some cases, it leads to lithospheric delamination and time-dependent, small-scale convection. We conclude that melting-induced buoyancy effects may explain the magmatic asymmetry and variations in lithospheric thickness that are inferred from observations.
more » « less- NSF-PAR ID:
- 10371579
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Geophysical Journal International
- Volume:
- 231
- Issue:
- 3
- ISSN:
- 0956-540X
- Page Range / eLocation ID:
- p. 2055-2078
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
Abstract Geodynamic simulations underpin our understanding of upper‐mantle processes, but their predictions require validation against observational data. Widely used geophysical datasets provide limited constraints on dynamic processes into the geological past, whereas under‐exploited geochemical observations from volcanic lavas at Earth's surface constitute a valuable record of mantle processes back in time. Here, we describe a new peridotite‐melting parameterization, BDD21, that can predict the incompatible‐element concentrations of melts within geodynamic simulations, thereby providing a means to validate these simulations against geochemical datasets. Here, BDD21's functionality is illustrated using the Fluidity computational modeling framework, although it is designed so that it can be integrated with other geodynamic software. To validate our melting parameterization and coupled geochemical‐geodynamic approach, we develop 2‐D single‐phase flow simulations of melting associated with passive upwelling beneath mid‐oceanic ridges and edge‐driven convection adjacent to lithospheric steps. We find that melt volumes and compositions calculated for mid‐oceanic ridges at a range of mantle temperatures and plate spreading rates closely match those observed at present‐day ridges with the same conditions. Our lithospheric step simulations predict spatial and temporal melting trends that are consistent with those recorded at intraplate volcanic provinces in similar geologic settings. Taken together, these results suggest that our coupled geochemical‐geodynamic approach can accurately predict a suite of present‐day geochemical observations. Since our results are sensitive to small changes in upper‐mantle thermal and compositional structure, this novel approach provides a means to improve our understanding of the mantle's thermo‐chemical structure and flow regime into the geological past.
-
more » « less
Abstract Gravity, magnetic, and bathymetry data collected along a continuous 1,400‐km‐long spreading‐parallel flow line across the Mid‐Atlantic Ridge indicate significant tectonic and magmatic fluctuations in the formation of oceanic crust over a range of time scales. The transect spans from 28 Ma on the African Plate to 74 Ma on the North American plate, crossing the Mid‐Atlantic Ridge at 35.8°N. Gravity‐derived crustal thicknesses vary from 3–9 km with a standard deviation of 1.0 km. Spectral analysis of bathymetry and residual mantle Bouguer anomaly show a diffuse power at >1 Myr and concurrent peaks at 390, 550, and 950 kyr. Large‐scale (>10 km) mantle thermal and compositional heterogeneities, variations in upper mantle flow, and detachment faulting likely generate the >1 Myr diffuse power. The 550‐ and 950‐kyr peaks may reflect the presence of magma solitons and/or regularly spaced ~7.7 and 13.3 km short‐wavelength mantle compositional heterogeneities. The 390‐kyr spectral peak corresponds to the characteristic spacing of faults along the flow line. Fault spacing also varies over longer periods (>10 Myr), which we interpret as reflecting long‐lived changes in the fraction of tectonically versus magmatically accommodated extensional strain. A newly discovered off‐axis oceanic core complex (Kafka Dome) found at 8 Ma on the African plate further suggests extended time periods of tectonically‐dominated plate separation. Fault spacing negatively correlates with gravity‐derived crustal thickness, supporting a strong link between magma input and fault style at mid‐ocean ridges.
-
Global Trends of Axial Relief and Faulting at Plate Spreading Centers Imply Discrete Magmatic Eventsmore » « less
Abstract Observed variations in across‐axis topographic relief and faulting style at spreading centers have been challenging to explain. Axial highs are seen at fast‐spreading centers, while valleys occur for slow‐spreading centers. Fault offsets range from tens of meters at fast‐spreading ridges to tens of kilometers at some slow‐spreading ridges. Models that fit the axial relief fail to produce observed fault patterns, while models that fit the fault patterns fail to produce observed variations in axial relief. A recent mechanical analysis (Liu & Buck, 2018,
https://doi.org/10.1016/j.epsl.2018.03.045 ) suggests that including the effect of many discrete diking events can result in a gradual change in axial relief with crustal thicknesses. To compare this mechanical model directly with observations requires us to couple it with a two‐dimensional thermal model. This allows us to estimate the axial lithospheric thickness consistently as a function of the spreading rate and crustal thickness. For thinner axial lithosphere the model predicts an axial high with relief supported by low‐density material beneath the axial lithosphere. For axial lithospheric thickness between approximately one half and approximately three fourths of the crustal thickness, the axial depth decreases with magma supply increase. For thicker axial lithosphere the axial valley relief is controlled by axial brittle lithospheric thickness and near‐axis lithospheric geometry. We compared model predictions to data by compiling observations on axial relief and faulting mode for all spreading centers where seismic crustal thickness has been measured. Good fit to the data is obtained for model parameters giving dike widths in the axial lithosphere close to a meter. -
more » « less
Abstract We studied long‐term evolution of nontransform discontinuities (NTDs) on the Mid‐Atlantic Ridge from 0‐ to ~20‐ to 25‐Ma crust using plate reconstructions of multibeam bathymetry, long‐range HMR1 sidescan sonar, residual mantle Bouguer gravity anomaly (RMBA), and gravity‐derived crustal thickness. NTDs have propagated north and south with respect to flowlines of relative plate motion and both rapidly and slowly compared to the half spreading rate; at times they have been quasi‐stable. Fast, short‐term (<2 Myr) propagation is driven by reduced magma supply (increased tectonic extension) in the propagating ridge tip when NTD ridge‐axis offsets are small (≲5 km). Propagation at larger offsets generally is slower and longer term. These NTDs can show classic structures of rift propagation including inner and outer pseudofaults and crustal blocks transferred between ridge flanks by discontinuous jumps of the propagating ridge tip. In all cases crustal transfer occurs within the NTD valley. Aside from ridge‐axis offset, the evolution of NTDs appears to be controlled by three factors: (1) gross volume and distribution of magma supplied to ridge segments as controlled by 3‐D heterogeneities in mantle fertility and/or dynamic upwelling; this controls fundamental ridge segmentation. (2) The lithospheric plumbing system through which magma is delivered to the crust. (3) The consequent focusing of tectonic extension in magma‐poor parts of spreading segments, typically at segment ends, which can drive propagation. We also observe long‐wavelength (5‐10 Myr) RMBA asymmetry between the conjugate ridge flanks, and we attribute this to asymmetric distribution of density anomalies in the upper mantle.
