An official website of the United States government
Abstract The rheology of oceanic lithosphere is important to our understanding of mantle dynamics and to the emergence and manifestations of plate tectonics. Data from experimental rock mechanics suggest rheology is dominated by three different deformation mechanisms including frictional sliding, low‐temperature plasticity, and high‐temperature creep, from shallow depths at relatively cold temperatures to large depths at relatively high temperatures. However, low‐temperature plasticity is poorly understood. This study further constrains low‐temperature plasticity by comparing observations of flexure at the Hawaiian Islands to predictions from 3‐D viscoelastic loading models with a realistic lithospheric rheology of frictional sliding, low‐temperature plasticity, and high‐temperature creep. We find that previously untested flow laws significantly underpredict the amplitude and overpredict the wavelength of flexure at Hawaii. These flow laws can, however, reproduce observations if they are weakened by a modest reduction (25–40%) in the plastic activation energy. Lithospheric rheology is strongly temperature dependent, and so we explore uncertainties in the thermal structure with different conductive cooling models and convection simulations of plume‐lithosphere interactions. Convection simulations show that thermal erosion from a plume only perturbs the lithospheric temperature significantly at large depths so that when it is added to the thermal structure, it produces a small increase in deflection. In addition, defining the temperature profile by the cooling plate model produces only modest weakening relative to the cooling half‐space model. Therefore, variation of the thermal structure does not appear to be a viable means of bringing laboratory‐derived flow laws for low‐temperature plasticity into agreement with geophysical field observations and modeling.
more »
« less
Reconciling lithospheric rheology between laboratory experiments, field observations and different tectonic settings
SUMMARY Recent modelling studies have shown that laboratory-derived rheology is too strong to reproduce observations of flexure at the Hawaiian Islands, while the same rheology appears consistent with outer rise—trench flexure at circum-Pacific subduction zones. Collectively, these results indicate that the rheology of an oceanic plate boundary is stronger than that of its interior, which, if correct, presents a challenge to understanding the formation of trenches and subduction initiation. To understand this dilemma, we first investigate laboratory-derived rheology using fully dynamic viscoelastic loading models and find that it is too strong to reproduce the observationally inferred elastic thickness, Te, at most plate interior settings. The Te can, however, be explained if the yield stress of low-temperature plasticity is significantly reduced, for example, by reducing the activation energy from 320 kJ mol−1, as in Mei et al., to 190 kJ mol−1 as was required by previous studies of the Hawaiian Islands, implying that the lithosphere beneath Hawaii is not anomalous. Second, we test the accuracy of the modelling methods used to constrain the rheology of subducting lithosphere, including the yield stress envelope (YSE) method, and the broken elastic plate model (BEPM). We show the YSE method accurately reproduces the model Te to within ∼10 per cent error with only modest sensitivity to the assumed strain rate and curvature. Finally, we show that the response of a continuous plate is significantly enhanced when a free edge is introduced at or near an edge load, as in the BEPM, and is sensitive to the degree of viscous coupling at the free edge. Since subducting lithosphere is continuous and generally mechanically coupled to a sinking slab, the BEPM may falsely introduce a weakness and hence overestimate Te at a trench because of trade-off. This could explain the results of recent modelling studies that suggest the rheology of subducting oceanic plate is stronger than that of its interior. However, further studies using more advanced thermal and mechanical models will be required in the future in order to quantify this.
more »
« less
- Award ID(s):
- 1940026
- PAR ID:
- 10391238
- Date Published:
- Journal Name:
- Geophysical Journal International
- Volume:
- 228
- Issue:
- 2
- ISSN:
- 0956-540X
- Page Range / eLocation ID:
- 857 to 875
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract The Hawaiian Ridge has long been a focus site for studying lithospheric flexure due to intraplate volcano loading, but crucial load and flexure details remain unclear. We address this problem using wide‐angle seismic refraction and reflection data acquired along a ∼535‐km‐long profile that intersects the ridge between the islands of Maui and Hawai'i and crosses 80–95 Myr‐old lithosphere. A tomographic image constructed using travel time data of several seismic phases reveals broad flexure of Pacific oceanic crust extending up to ∼200–250 km either side of the Hawaiian Ridge, and vertically up to ∼6–7 km. TheP‐wave velocity structure, verified by gravity modeling, reveals that the west flank of Hawaii is comprised of extrusive lavas overlain by volcanoclastic sediments and a carbonate platform. In contrast, the Hāna Ridge, southeast of Maui, contains a high‐velocity core consistent with mafic or ultramafic intrusive rocks. Magmatic underplating along the seismic line is not evident. Reflectors at the top and bottom of the pre‐existing oceanic crust suggest a ∼4.5–6 km crustal thickness. Simple three‐dimensional flexure modeling with an elastic plate thickness,Te, of 26.7 km shows that the depths to the reflectors beneath the western flank of Hawai'i can be explained by volcano loading in which Maui and the older islands in the ridge contribute ∼43% to the flexure and the island of Hawai'i ∼51%. Previous studies, however, revealed a higherTebeneath the eastern flank of Hawai'i suggesting that isostatic compensation may not yet be complete at the youngest end of the ridge.more » « less
-
Subduction initiation remains one of the least understood global processes of plate tectonics. Prominent models have been cast in terms of two broad classes: spontaneous cases due to lithospheric gravitational instabilities and induced cases due to forced plate convergence. Yet gravitationally unstable lithosphere is old, strong, and difficult to begin to bend into a subduction zone and convergent forces necessary to begin subduction are often too large given the plates involved. These models also consider the asthenospheric mantle as passive, even though relative motion between slabs and the asthenosphere has long been regarded as a strong control on subduction dynamics. Here I propose that subduction-transform edge propagator (STEP) faults can initiate subduction depending on the absolute motion of the STEP fault with respect to the asthenosphere. STEP faults form where subduction zones end and the subducting plate tears forming a down flexed transcurrent plate boundary at the surface shearing against the adjacent rear arc lithospheric plate. However, STEP faults are not simple transcurrent boundaries. Absolute motion of the down flexed STEP fault edge with respect to the surrounding asthenosphere can produce a strong sea anchor force that either continues to bend the edge downward, initiating subduction, or opposes slab bending, inhibiting subduction. In the south Pacific, the southern end of the New Hebrides Trench and the northern end of the Tonga Trench are type-example STEP faults with opposite senses of dip but both moving northward with respect to the asthenosphere. The northward dipping New Hebrides STEP fault moves northward in a mantle reference frame creating a strong asthenospheric flow against the STEP fault edge, inducing active subduction at the Matthew-Hunter trench. In contrast, the Tonga STEP fault dips southward but also has a northward component of motion with respect to the mantle. Asthenosphere thus flows southward beneath the down flexed Tonga STEP fault edge opposing further bending. Subduction does not initiate at the Tonga STEP fault despite a ~100 Myr age contrast between the Pacific and north Fiji and Lau basin lithospheres. Since absolute plate motions reflect the sum of all forces acting on the entire lithospheric plate, a strong sea anchor mantle force may be generated at a STEP fault edge, initiating subduction (or inhibiting it), even where lithosphere is old, strong, and resists bending and without requiring large convergent forces between plates, overcoming these objections to previous models.more » « less
-
Abstract During 1975–1988, an academic research ship, R/VRobert D.Conrad, acquired more than 150,000‐line‐km of multichannel seismic reflection profile data from each of the world's main ocean basins and their margins. This extensive legacy seismic data set, which involved both single ship and two‐ship data acquisition, has been widely used by the marine geoscience community. We report on our experience in reprocessing seismic reflection profile data acquired duringConradcruise RC2308 to the Hawaiian Islands region in August/September 1982. We show that the application of modern, industry standard processing techniques, including filtering, de‐bubble, deconvolution, and migration, can significantly enhance 40+ year old legacy seismic reflection profile data. The reprocessed data reveals more precisely, and with much less scatter, the flexure of Cretaceous Pacific oceanic crust caused by the Pliocene‐Recent volcanic loads that comprise the Hawaiian Islands. A comparison of observed picks of top oceanic crust which has been corrected for the Hawaiian swell and the Molokai Fracture Zone with the calculations of a simple 3‐dimensional elastic plate (flexure) model reveals a best fit elastic plate thickness of the lithosphere,Te, of 26.7 km, an average infill density of 2,701 kg m−3, and a Root Mean Square difference between observations and calculations of 305 m. Tests show these results depend weakly on the load density assumed and that the average infill density is close to what would be predicted from an arithmetic average of the flanking moat infill density and the infill density that immediately underlies the volcanic edifice.more » « less
-
SUMMARY Tectonic plate motions predominantly result from a balance between the potential energy change of the subducting slab and viscous dissipation in the mantle, bending lithosphere and slab–upper plate interface. A wide range of observations from active subduction zones and exhumed rocks suggest that subduction interface shear zone rheology is sensitive to the composition of subducting crustal material—for example, sediments versus mafic igneous oceanic crust. Here we use 2-D numerical models of dynamically consistent subduction to systematically investigate how subduction interface viscosity influences large-scale subduction kinematics and dynamics. Our model consists of an oceanic slab subducting beneath an overriding continental plate. The slab includes an oceanic crustal/weak layer that controls the rheology of the interface. We implement a range of slab and interface strengths and explore how the kinematics respond for an initial upper mantle slab stage, and subsequent quasi-steady-state ponding near a viscosity jump at the 660-km-discontinuity. If material properties are suitably averaged, our results confirm the effect of interface strength on plate motions as based on simplified viscous dissipation analysis: a ∼2 order of magnitude increase in interface viscosity can decrease convergence speeds by ∼1 order of magnitude. However, the full dynamic solutions show a range of interesting behaviour including an interplay between interface strength and overriding plate topography and an end-member weak interface-weak slab case that results in slab break-off/tearing. Additionally, for models with a spatially limited, weak sediment strip embedded in regular interface material, as might be expected for the subduction of different types of oceanic materials through Earth’s history, the transient response of enhanced rollback and subduction velocity is different for strong and weak slabs. Our work substantiates earlier suggestions as to the importance of the plate interface, and expands the range of quantifiable links between plate reorganizations, the nature of the incoming and overriding plate and the potential geological record.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1093/gji/ggab382
