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.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


This content will become publicly available on November 1, 2026

Title: The Effect of Composition on Shear Localization in Planetary Lithospheres
Earth's particular style of plate‐tectonics—characterized by localized deformation along dynamic plate boundaries and long‐lived stable plate interiors—appears to be unique among rocky objects in the solar system. However, it is entirely unknown how common plate tectonics and related lithospheric phenomena are among the vast population of exoplanets discovered astronomically or assumed to exist throughout the Universe. In this study, we explore the effect of planetary composition on mylonitization—a set of microphysical processes that is commonly associated with shear localization and plate boundary deformation on Earth. A model for planet compositions, based on stellar spectroscopy, is used to define a plausible range of theoretical mineral abundances in the mantles of rocky Earth‐sized exoplanets. These mineral abundances, along with experimental rock rheology, are used to model microphysical evolution with two‐phase mixing. The model is then used to determine the effect of composition on the time‐scales for shear zone formation. We demonstrate that lithospheres composed of sub‐equal proportions of two mineral phases will form shear zones over relatively short time‐scales, a more favorable condition for forming Earth‐like plate boundaries. In contrast, lithospheres that are nearly monomineralic may require unrealistically long time‐scales to form plate boundary shear zones. Using this approach, we identify specific nearby stars with the optimal range of compositions to be targeted by future astronomical missions, including the Habitable Worlds Observatory.  more » « less
Award ID(s):
2224725
PAR ID:
10647886
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
American Geophysical Union
Date Published:
Journal Name:
Journal of Geophysical Research: Planets
Volume:
130
Issue:
11
ISSN:
2169-9097
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. (, Philosophical transactions - Royal Society. Mathematical, physical and engineering sciences)
    For plate tectonics to operate on a planet, mantle convective forces must be capable of forming weak, localized shear zones in the lithosphere that act as plate boundaries. Otherwise, a planet's mantle will convect in a stagnant lid regime, where subduction and plate motions are absent. Thus, when and how plate tectonics initiated on Earth is intrinsically tied to the ability of mantle convection to form plate boundaries; however, the physics behind this process are still uncertain. Most mantle convection models have employed a simple pseudoplastic model of the lithosphere, where the lithosphere "fails" and develops a mobile lid when stresses in the lithosphere reach the prescribed yield stress. With pseudoplasticity high mantle temperatures and high rates of internal heating, conditions relevant for the early Earth, impede plate boundary formation by decreasing lithospheric stresses, and hence favor a stagnant lid for the early Earth. However, when a model for shear zone formation based on grain size reduction is used, early Earth thermal conditions do not favor a stagnant lid. While lithosphere stress drops with increasing mantle temperature or heat production rate, the deformational work, which drives grain size reduction, increases. Thus the ability of convection to form weak plate boundaries is not impeded by early Earth thermal conditions. However, mantle thermal state does change the style of subduction and lithosphere mobility; high mantle temperatures lead to a more sluggish, drip-like style of subduction. This "sluggish lid" convection may be able to explain many of the key observations of early Earth crust formation processes preserved in the geologic record. Moreover, this work highlights the importance of understanding the microphysics of plate boundary formation for assessing early Earth tectonics, as different plate boundary formation mechanisms are influenced by mantle thermal state in fundamentally different ways. 
    more » « less
  2. ; ; ; (, Geophysical Journal International)
    SUMMARY The occurrence of plate tectonics on Earth is rooted in the physics of lithospheric ductile weakening and shear-localization. The pervasiveness of mylonites at lithospheric shear zones is a key piece of evidence that localization correlates with reduction in mineral grain size. Most lithospheric mylonites are polymineralic and the interaction between mineral phases, such as olivine and pyroxene, especially through Zener pinning, impedes normal grain growth while possibly enhancing grain damage, both of which facilitate grain size reduction and weakening, as evident in lab experiments and field observations. The efficacy of pinning, however, relies on the mineral phases being mixed and dispersed at the grain scale, where well-mixed states lead to greater mylonitization. To model grain mixing between different phases at the continuum scale, we previously developed a theory treating grain-scale processes as diffusion between phases, but driven by imposed compressive stresses acting on the boundary between phases. Here we present a new model for shearing rock that combines our theory for diffusive grain mixing, 2-D non-Newtonian flow and two-phase grain damage. The model geometry is designed specifically for comparison to torsional shear-deformation experiments. Deformation is either forced by constant velocity or constant stress boundary conditions. As the layer is deformed, mixing zones between different mineralogical units undergo enhanced grain size reduction and weakening, especially at high strains. For constant velocity boundary experiments, stress drops towards an initial piezometric plateau by a strain of around 4; this is also typical of monophase experiments for which this initial plateau is the final steady state stress. However, polyphase experiments can undergo a second large stress drop at strains of 10–20, and which is associated with enhanced phase mixing and resultant grain size reduction and weakening. Model calculations for polyphase media with grain mixing and damage capture the experimental behaviour when damage to the interface between phases is moderately slower or less efficient than damage to the grain boundaries. Other factors such as distribution and bulk fraction of the secondary phase, as well as grain-mixing diffusivity also influence the timing of the second stress drop. For constant stress boundary conditions, the strain rate increases during weakening and localization. For a monophase medium, there is theoretically one increase in strain rate to a piezometric steady state. But for the polyphase model, the strain rate undergoes a second abrupt increase, the timing for which is again controlled by interface damage and grain mixing. The evolution of heterogeneity through mixing and deformation, and that of grain size distributions also compare well to experimental observations. In total, the comparison of theory to deformation experiments provides a framework for guiding future experiments, scaling microstructural physics to geodynamic applications and demonstrates the importance of grain mixing and damage for the formation of plate tectonic boundaries. 
    more » « less
  3. ; ; ; ; ; ; (, Basin Research)
    Abstract Rift initiation within cold, thick, strong lithosphere and the evolving linkage to form a contiguous plate boundary remains debated in part owing to the lack of time–space constraints on kinematics of basement‐involved faults. Different rift sectors initiate diachronously and may eventually link to produce a jigsaw spatial pattern, as in the East African rift, and along the Atlantic Ocean margins. The space–time distribution of earthquakes illuminates the geometry and kinematics of fault zones within the crystalline crust, as well as areas with pressurized magma bodies. We use seismicity and Global Navigation System Satellites (GNSS) data from the Turkana Rift Array Investigating Lithospheric Structure (TRAILS) project in East Africa and a new digital compilation of faults and eruptive centres to evaluate models for the kinematic linkage of two initially separate rift sectors: the Main Ethiopian Rift (MER) and the Eastern rift (ER). The ca. 300 km wide zone of linkage includes failed basins and linkage zones; seismicity outlines active structures. Models of GNSS data indicate that the ca. 250 km‐wide zone of seismically active en echelon basins north of the Turkana Depression is a zone, or block, of distributed strain with small counterclockwise rotation that serves to connect the Main Ethiopian and Eastern rifts. Its western boundary is poorly defined owing to data gaps in South Sudan. Strain across the northern and southern boundaries of this block, and an ca. 50 km‐wide kink in the southern Turkana rift is accommodated by en echelon normal faults linked by short strike‐slip faults in crystalline basement, and relay ramps at the surface. Short segments of obliquely oriented basement structures facilitate across‐rift linkage of faults, but basement shear zones and Mesozoic rift faults are not actively straining. This configuration has existed for at least 2–5 My without the development of localized shear zones or transform faults, documenting the importance of distributed deformation in continental rift tectonics. 
    more » « less
  4. ; (, Geochemistry, Geophysics, Geosystems)
    Earth's style of planetary heat transport is characterized by plate tectonics which requires rock strength to be reduced plastically in order to break an otherwise stagnant lithospheric lid, and for rocks to have a memory of past deformation to account for strain localization and the hysteresis implied by geological sutures. Here, we explore ∼107Rayleigh number, visco‐plastic, 3‐D global mantle convection with damage. We show that oceanic lithosphere‐only models generate strong toroidal‐poloidal power ratios and features such as a mix of long‐wavelength tectonic motions and smaller‐scale, back‐arc tectonics driven by downwellings. Undulating divergent plate boundaries can evolve to form overlapping spreading centers and microplates, promoted and perhaps stabilized by the effects of damage with long memory. The inclusion of continental rafts enhances heat flux variability and toroidal flow, including net rotation of the lithosphere, to a level seen in plate reconstructions for the Cenozoic. Both the super‐continental cycle and local rheological descriptions affect heat transport and tectonic deformation across a range of scales, and we showcase both general tectonic dynamics and regionally applied continental breakup scenarios. Our work points toward avenues for renewed analysis of the typical, mean behavior as well as the evolution of fluctuations in geological and model plate boundary evolution scenarios. 
    more » « less
  5. ; ; ; ; (, Journal of Geophysical Research: Solid Earth)
    Shear wave splitting of teleseismic core phases such as SKS is commonly used to constrain mantle seismic anisotropy, a proxy for convective deformation. In plate boundaries, sharp lateral variations of splitting measurements near transform faults are often linked to deformation within a lithospheric shear zone below, but potential seismic waveform effects from heterogeneous structure on small scales may influence the interpretation. Here, we explore possible finite frequency effects on shear wave splitting near fault zones in a fully three‐dimensional anisotropic setting. We find that shear zones wider than 80 km, a scale set by the Fresnel zone, can be clearly detected, but narrower zones are less distinguishable. Near the edge of the shear zone, the combined effect of anisotropy and scattering generates false splitting measurements with large delay times and fast axis orientation approaching the back‐azimuth, a bias which can only be identified when records from different back‐azimuths are analyzed together. This substantiates that back‐azimuthal variations of splitting can arise not just from vertical layering but also lateral changes of anisotropic media. We also test the effects of shear zone edge geometry, epicentral distance, filtering frequency, crustal thickness, and sediment cover. Our study delineates the ability of shear wave splitting to resolve and investigate fault zones, and emphasizes the importance of good azimuthal coverage to correctly interpret observed anisotropy. Based on revisiting previous shear wave splitting and lithospheric deformation studies, we infer that many crustal fault zones are underlain by lithospheric shear zones at least 20 km wide. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email