Breaking Earth's shell into a global plate network.

The initiation mechanism of Earth's plate tectonic cooling system remains uncertain. A growing consensus suggests that multi-plate tectonics was preceded by cooling through a single-plate lithosphere, but models for how this lithosphere was first broken into plates have not converged on a mechanism or a typical early plate scale. A commonality among prior efforts is the use of continuum mechanics approximations to evaluate this solid mechanics problem. Here we use 3D spherical shell models to demonstrate a self-organized fracture mechanism analogous to thermal expansion-driven lithospheric uplift, in which globe-spanning rifting occurs as a consequence of horizontal extension. Resultant fracture spacing is a function of lithospheric thickness and rheology, wherein geometrically-regular, polygonal-shaped tessellation is an energetically favored solution because it minimizes total crack length. Therefore, warming of the early lithosphere itself-as anticipated by previous studies-should lead to failure, propagating fractures, and the conditions necessary for the onset of multi-plate tectonics.

Fracture spacing in layered materials and pattern transition from parallel to polygonal fractures.

We perform three-dimensional simulations of fracture growth in a three-layered plate model with an embedded heterogeneous layer under horizontal biaxial stretch (representing stretch from directional to isotropic) by the finite element approach. The fractures develop under a quasistatical, slowly increasing biaxial strain. The material inhomogeneities are accounted for by assigning each element a failure threshold that is defined by a given statistical distribution. A universal scale law of fracture spacing to biaxial strain in terms of principal stress ratio is well demonstrated in a three-dimensional fashion. The numerically obtained fracture patterns show a continuous pattern transition from parallel fractures, laddering fracture to polygonal fractures, which depends strongly on the far-field loading conditions in terms of principal stress ratio lambda = sigma(2)/sigma(1), from uniaxial (lambda = 0), anisotropic (0 < lambda < 1) to isotropic stretch (lambda = 1). We find that, except for further opening of existing fractures after they are well-developed (saturation), new fractures may also initiate and propagate along the interface between layers, which may serve as another mechanism to accommodate additional strain for fracture saturated layers.