skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Continuum approximation of dyking with a theory for poro-viscoelastic–viscoplastic deformation
SUMMARY 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
Award ID(s):
2121666 2121568
PAR ID:
10414412
Author(s) / Creator(s):
; ; ; ;
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
  1. The strength of lithospheric plates is a central component of plate tectonics, governed by brittle processes in the shallow portion of the plate and ductile behavior in the deeper portion. We review experimental constraints on ductile deformation of olivine, the main mineral in the upper mantle and thus the lithosphere. Olivine deforms by four major mechanisms: low-temperature plasticity, dislocation creep, dislocation-accommodated grain-boundary sliding (GBS), and diffusion-accommodated grain-boundary sliding (diffusion creep). Deformation in most of the lithosphere is dominated by GBS, except in shear zones—in which diffusion creep dominates—and in the brittle-ductile transition—in which low-temperature plasticity may dominate. We find that observations from naturally deformed rocks are consistent with extrapolation of the experimentally constrained olivine flow laws to geological conditions but that geophysical observations predict a weaker lithosphere. The causes of this discrepancy are unresolved but likely reside in the uncertainty surrounding processes in the brittle-ductile transition, at which the lithosphere is strongest. ▪ Ductile deformation of the lithospheric mantle is constrained by experimental data for olivine. ▪ Olivine deforms by four major mechanisms: low-temperature plasticity, dislocation creep, dislocation-accommodated grain-boundary sliding, and diffusion creep. ▪ Observations of naturally deformed rocks are consistent with extrapolation of olivine flow laws from experimental conditions. ▪ Experiments predict stronger lithosphere than geophysical observations, likely due to gaps in constraints on deformation in the brittle-ductile transition. 
    more » « less
  2. 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. 
    more » « less
  3. Fatigue analysis in metallic frame structures can be challenging due to associated computational costs; if localized plasticity is involved, then the approach of three-dimensional (3D) continuum plasticity models for direct computation of stresses will be infeasible for the analysis of cyclic loading that would need to be modeled in medium- to high-cycle fatigue and vibratory fatigue applications. This difficulty is particularly accentuated in architected structures, for which high-resolution 3D finite element analysis (FEA) would be prohibitively expensive. In this work, we propose an alternative approach based on the use of novel elasto-plastic frame model with continuous flow (i.e. no sharp yield function) for modeling 3D frame and lattice structures. Rather than splitting the strains (as is done in classical plasticity) we split the deformation measures, extension, curvature and twist, into elastic and plastic components and postulate a rate type evolution rule for the plastic variables in terms of the stress resultants (axial force, bending moment, and torque). The combination of structural models together with the use of elasto-plastic operator split to solve the resulting boundary value problem allows for much faster determination of localized plasticity than continuum models can provide. The use of a continuous transition from elastic to rate-independent plasticity (as opposed to an abrupt change with classical plasticity models) allows us to capture localized microplasticity and determine resulting fatigue progression using a cycle-count-free, plastic work-based approach, formulated in terms of the curvatures and resultants. We demonstrate that (a) the model is able able to reproduce the response of 3D FEA with very few elements and (b) the model has the ability to rapidly predict the fatigue life under variable amplitude combined loading with relatively few frame elements. 
    more » « less
  4. 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
  5. 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