Strain weakening enables continental plate tectonics

International audienceMuch debate exists concerning the strength distribution of the continental lithosphere, how it controls lithosphere-scale strain localization and hence enables plate tectonics. No rheological model proposed to date is comprehensive enough to describe both the weakness of plate boundary and rigid-like behaviour of plate interiors. Here we show that the duality of strength of the lithosphere corresponds to different stages of microstructural evolution. Geological constraints on lithospheric strength and large strain numerical experiments reveal that the development of layers containing weak minerals and the onset of grain boundary sliding upon grain size reduction in olivine cause strain localisation and reduce strength in the crust and subcontinental mantle, respectively. The positive feedback between weakening and strain localization leads to the progressive development of weak plate boundaries while plate interiors remai

MeltMigrator : a MATLAB-based software for modeling three-dimensional melt migration and crustal thickness variations at mid-ocean ridges following a rules-based approach

Author Posting. © American Geophysical Union, 2017. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Geochemistry, Geophysics, Geosystems 18 (2017): 445–456, doi:10.1002/2016GC006686.MeltMigrator is a MATLAB®-based melt migration software developed to process three-dimensional mantle temperature and velocity data from user-supplied numerical models of mid-ocean ridges, calculate melt production and melt migration trajectories in the mantle, estimate melt flux along plate boundaries, and predict crustal thickness distribution on the seafloor. MeltMigrator is also capable of calculating compositional evolution depending on the choice of petrologic melting model. Programmed in modules, MeltMigrator is highly customizable and can be expanded to a wide range of applications. We have applied it to complex mid-ocean ridge model settings, including transform faults, oblique segments, ridge migration, asymmetrical spreading, background mantle flow, and ridge-plume interaction. In this technical report, we include an example application to a segmented mid-ocean ridge. MeltMigrator is available as a supplement to this paper, and it is also available from GitHub and the University of Maryland Geodynamics Group website.National Science Foundation Grant Number: OCE-0937277 and OCE-14582012017-07-2

A probabilistic damage model of stress-induced permeability anisotropy during cataclastic flow

Author Posting. © American Geophysical Union, 2007. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research 112 (2007): B10207, doi:10.1029/2006JB004456.A fundamental understanding of the effect of stress on permeability evolution is important for many fault mechanics and reservoir engineering problems. Recent laboratory measurements demonstrate that in the cataclastic flow regime, the stress-induced anisotropic reduction of permeability in porous rocks can be separated into 3 different stages. In the elastic regime (stage I), permeability and porosity reduction are solely controlled by the effective mean stress, with negligible permeability anisotropy. Stage II starts at the onset of shear-enhanced compaction, when a critical yield stress is attained. In stage II, the deviatoric stress exerts primary control over permeability and porosity evolution. The increase in deviatoric stress results in drastic permeability and porosity reduction and considerable permeability anisotropy. The transition from stage II to stage III takes place progressively during the development of pervasive cataclastic flow. In stage III, permeability and porosity reduction becomes gradual again, and permeability anisotropy diminishes. Microstructural observations on deformed samples using laser confocal microscopy reveal that stress-induced microcracking and pore collapse are the primary forms of damage during cataclastic flow. A probabilistic damage model is formulated to characterize the effects of stress on permeability and its anisotropy. In our model, the effects of both effective mean stress and differential stress on permeability evolution are calculated. By introducing stress sensitivity coefficients, we propose a first-order description of the dependence of permeability evolution on different loading paths. Built upon the micromechanisms of deformation in porous rocks, this unified model provides new insight into the coupling of stress and permeability.W.Z. was partially supported by the National Science Foundation under grants NSF-OCE0221436 and NSF-EAR 0510459, and the Department of Energy under grant #DEFGO200ER15058 (WHOI). LM was supported by the National Science Foundation under grant NSF-EAR0337678

The recent history of the Galapagos Triple Junction preserved on the Pacific plate

Author Posting. © The Author(s), 2013. This is the author's version of the work. It is posted here by permission of Elsevier B.V. for personal use, not for redistribution. The definitive version was published in Earth and Planetary Science Letters 371-372 (2013): 6-15, doi:10.1016/j.epsl.2013.04.018.At the Galapagos triple junction, the Cocos and Nazca plates are broken by a succession of transient rifts north and south of the Cocos-Nazca (C-N) Rift. Modeling has suggested that each rift initiated at the East Pacific Rise (EPR), its location controlled by the distance of the C-N Rift tip from the EPR. Evidence on the Pacific plate confirms that each transient rift formed a true RRR triple junction with the EPR and clarifies the history of the region. At ~1.5 Ma the triple junctions began jumping rapidly toward C-N Rift suggesting that the C-N Rift tip moved closer to the EPR. Pacific abyssal hills became broad and shallow indicating enhanced magma supply to the region. At ~1.4 Ma, the Galapagos microplate developed when extension became fixed on the southern transient rift to form the South scarp of the future Dietz rift basin. Lavas flooded the area and a Galapagos-Nazca magmatic spreading center initiated at the EPR. We suggest that a hotspot was approaching the southern triple junction from the west. The hotspot crossed to the Nazca plate ~1.25 Ma. Dietz seamount formed within the young spreading center, dikes intruded Dietz rift basin, and eruptions built volcanic ridges. Since ~0.8 Ma magmatic spreading has jumped northward twice, most recently to Dietz volcanic ridge. Amagmatic extension to the east has formed the large North scarp of Dietz rift basin. Northward jumping of the southern triple junction has maintained the microplate boundary close to the proposed hotspot.DKS was partially supported by NSF grant OCE-1028537, WZ by NSF grant EAR-1056317, and LM by NSF grant OCE-1060878

Controls on melt migration and extraction at the ultraslow Southwest Indian Ridge 10°–16°E

Author Posting. © American Geophysical Union, 2011. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research 116 (2011): B10102, doi:10.1029/2011JB008259.Crustal thickness variations at the ultraslow spreading 10–16°E region of the Southwest Indian Ridge are used to constrain melt migration processes. In the study area, ridge morphology correlates with the obliquity of the ridge axis with respect to the spreading direction. A long oblique “supersegment”, nearly devoid of magmatism, is flanked at either end by robust magmatic centers (Joseph Mayes Seamount and Narrowgate segment) of much lesser obliquity. Plate-driven mantle flow and temperature structure are calculated in 3-D based on the observed ridge segmentation. Melt extraction is assumed to occur in three steps: (1) vertical migration out of the melting region, (2) focusing along an inclined permeability barrier, and (3) extraction when the melt enters a region shallower than ∼35 km within 5 km of the ridge axis. No crust is predicted in our model along the oblique supersegment. The formation of Joseph Mayes Seamount is consistent with an on-axis melt anomaly induced by the local orthogonal spreading. The crustal thickness anomaly at Narrowgate results from melt extracted at a tectonic damage zone as it travels along the axis toward regions of lesser obliquity. Orthogonal spreading enhances the Narrowgate crustal thickness anomaly but is not necessary for it. The lack of a residual mantle Bouguer gravity high along the oblique supersegment can be explained by deep serpentization of the upper mantle permissible by the thermal structure of this ridge segment. Buoyancy-driven upwelling and/or mantle heterogeneities are not required to explain the extreme focusing of melt in the study area.This work was supported by grants OCE‐ 0623188 and OCE‐0937277 from the National Science Foundation

Experimental evidence for melt partitioning between olivine and orthopyroxene in partially molten harzburgite

Author Posting. © American Geophysical Union, 2016. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Solid Earth 121 (2016): 5776–5793, doi:10.1002/2016JB013122.Observations of dunite channels in ophiolites and uranium series disequilibria in mid-ocean ridge basalt suggest that melt transport in the upper mantle beneath mid-ocean ridges is strongly channelized. We present experimental evidence that spatial variations in mineralogy can also focus melt on the grain scale. This lithologic melt partitioning, which results from differences in the interfacial energies associated with olivine-melt and orthopyroxene-melt boundaries, may complement other melt focusing mechanisms in the upper mantle such as mechanical shear and pyroxene dissolution. We document here lithologic melt partitioning in olivine-/orthopyroxene-basaltic melt samples containing nominal olivine to orthopyroxene ratio of 3 to 2 and melt fractions of 0.02 to 0.20. Experimental samples were imaged using synchrotron-based X-ray microcomputed tomography at a resolution of 700 nm per voxel. By analyzing the local melt fraction distributions associated with olivine and orthopyroxene grains in each sample, we found that the melt partitioning coefficient, i.e., the ratio of melt fraction around olivine to that around orthopyroxene grains, varies between 1.1 and 1.6. The permeability and electrical conductivity of our digital samples were estimated using numerical models and compared to those of samples containing only olivine and basaltic melt. Our results suggest that lithologic melt partitioning and preferential localization of melt around olivine grains might play a role in melt focusing, potentially enhancing average melt ascent velocities.National Science Foundation Grant Numbers: 1250338, 1551300; Basic Energy Sciences Grant Number: DEFG0207ER15916; Advanced Photon Source Grant Number: DE-AC02-06CH113572017-02-2

Geodynamic Models of Melt Generation and Extraction at Mid-Ocean Ridges

It is widely accepted that plate divergence at mid-ocean ridges drives mantle flow, mantle melting, and the formation of new oceanic crust. However, many of the details of this process remain obscure because of the inaccessibility of the mantle to direct observation. Thus, geodynamic models are needed to provide insight into the processes that control the formation of new crust and hydrothermal circulation at mid-ocean ridges. These models allow us to test governing parameters and investigate physical hypotheses and conceptual models derived from geological, geophysical, and geochemical observations. During the span of the Ridge 2000 Program, a new generation of models was developed to calculate the width of the melt region and the extent of melting beneath mid-ocean ridges, track the pathways along which melts may migrate, and predict melt and residual mantle compositions as the system evolves. Findings from these studies illustrate the importance of melt focusing for the efficient delivery of melt to the ridge axis, the complexities of migrating melt in the vicinity of ridge offsets, and the effect of mantle rheology in model calculations.Keywords: Spreading rate dependence, Atlantic Ridge, East Pacific rise, Primary magmas, Asthenospheric flow, Mid-Ocean Ridges, Crustal thickness, Gravity anomalies, Deformable porous media, Off axis volcanis

montesi/RHEOL: Rheol development snapshot

Use at your own risk: not every feature is implemented. user interface by question and answer only. MATLAB code for calculating strength profiles Rheologies dataset build by build_rock.m, store in MATLAB matrix Driver code is StrengthProfile.m; calls other codes and datasets Some features may be broken right now. Especially, multiphase mixing is not currently operating. Keep one rock per layer