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.
To reach Earth’s surface, magma must ascend from the hot, ductile asthenosphere through cold and brittle rock in the lithosphere. It does so via fluid-filled fractures called dykes. While the continuum mechanics of ductile asthenosphere is well established, there has been little theoretical work on the cold and brittle regime where dyking and faulting occurs. Geodynamic models use plasticity to model fault-like behaviour; plasticity also shows promise for modelling dykes. Here we build on an existing model to develop a poro-viscoelastic–viscoplastic theory for two-phase flow across the lithosphere. Our theory addresses the deficiencies of previous work by incorporating (i) a hyperbolic yield surface, (ii) a plastic potential with control of dilatancy and (iii) a viscous regularization of plastic failure. We use analytical and numerical solutions to investigate the behaviour of this theory. Through idealized models and a comparison to linear elastic fracture mechanics, we demonstrate that this behaviour includes a continuum representation of dyking. Finally, we consider a model scenario reminiscent of continental rifting and demonstrate the consequences of dyke injection into the cold, upper lithosphere: a sharp reduction in the force required to rift.
more » « less- NSF-PAR ID:
- 10414412
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Geophysical Journal International
- Volume:
- 234
- Issue:
- 3
- ISSN:
- 0956-540X
- Page Range / eLocation ID:
- p. 2007-2031
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
Abstract Micromechanics‐based continuum damage criteria have previously been developed to simulate the initiation of ductile fracture in structural steels under conditions with large‐scale plasticity where conventional fracture mechanics indices are invalid. Such models have been combined with methods to simulate ductile crack growth for monotonic loading. In this study, a micromechanics‐based adaptive cohesive zone model for simulating ductile crack propagation under monotonic loading is extended to handle cyclic loading. The proposed model adaptively modifies the cohesive traction–separation relationship for crack opening and closure, as the loading reverses between tension and compression. The approach is implemented into the finite element analysis platform WARP3D, and results of simulations that use the model are compared with data from coupon‐scale tests. The results demonstrate that the proposed model can accurately simulate the effect of crack propagation on specimen response, as well as other key aspects of observed behavior, including crack face closure and crack tunneling.
-
more » « less
Abstract Flexure occurs on intermediate geologic timescales (∼1 Myr) due to volcanic‐island building at the Island of Hawaii, and the deformational response of the lithosphere is simultaneously elastic, plastic, and ductile. At shallow depths and low temperatures, elastic deformation transitions to frictional failure on faults where stresses exceed a threshold value, and this complex rheology controls the rate of deformation manifested by earthquakes. In this study, we estimate the seismic strain rate based on earthquakes recorded between 1960 and 2019 at Hawaii, and the estimated strain rate with 10−18–10−15s−1in magnitude exhibits a local minimum or neutral bending plane at 15 km depth within the lithosphere. In comparison, flexure and internal deformation of the lithosphere are modeled in 3D viscoelastic loading models where deformation at shallow depths is accommodated by frictional sliding on faults and limited by the frictional coefficient (
μ f), and at larger depths by low‐temperature plasticity and high‐temperature creep. Observations of flexure and the seismic strain rate are best‐reproduced by models withμ f = 0.3 ± 0.1 and modified laboratory‐derived low‐temperature plasticity. Results also suggest strong lateral variations in the frictional strength of faults beneath Hawaii. Our models predict a radial pattern of compressive stress axes relative to central Hawaii consistent with observations of earthquake pressure (P) axes. We demonstrate that the dip angle of this radial axis is essential to discerning a change in the curvature of flexure, and therefore has implications for constraining lateral variations in lithospheric strength. -
more » « less
Abstract Although the surface deformation of tectonic plate boundaries is well determined by geological and geodetic measurements, the pattern of flow below the lithosphere remains poorly constrained. We use the crustal velocity field of the Plate Boundary Observatory to illuminate the distribution of horizontal flow beneath the California margin. At lower-crustal and upper-mantle depths, the boundary between the Pacific and North American plates is off-centered from the San Andreas fault, concentrated in a region that encompasses the trace of nearby active faults. A major step is associated with return flow below the Eastern California Shear Zone, leading to the extrusion of the Mojave block and a re-distribution of fault activity since the Pleistocene. Major earthquakes in California have occurred above the regions of current plastic strain accumulation. Deformation is mechanically coupled from the crust to the asthenosphere, with mantle flow overlaid by a kinematically consistent network of faults in the brittle crust.
-
more » « less
Abstract The bottom of the lithosphere is characterized by a thermally controlled transition from brittle to ductile deformation. While the mechanical behavior of rocks firmly within the brittle and ductile regimes is relatively well understood, how the transition operates remains elusive. Here, we study the mechanical properties of pure olivine gouge from 100 to 500°C under 100 MPa pore‐fluid pressure in a triaxial deformation apparatus as a proxy for the mechanical properties of the upper mantle across the brittle‐ductile transition. We describe the mechanical data with a rate‐, state‐, and temperature‐dependent constitutive law with multiple thermally activated deformation mechanisms. The stress power exponents decrease from 70 ± 10 in the brittle regime to 17 ± 3 and 4 ± 2 in the semi‐brittle and ductile regimes, respectively. The mechanical model consistently explains the mechanical behavior of olivine gouge across the brittle‐ductile transition, capturing the gradual evolution from cataclasis to crystal plasticity.
