Subduction zones are pivotal for the recycling of Earth’s outer layer into its interior. However, the conditions under which new subduction zones initiate are enigmatic. Here, we constructed a transdisciplinary database featuring detailed analysis of more than a dozen documented subduction zone initiation events from the last hundred million years. Our initial findings reveal that horizontally forced subduction zone initiation is dominant over the last 100 Ma, and that most initiation events are proximal to pre-existing subduction zones. The SZI Database is expandable to facilitate access to the most current understanding of subduction zone initiation as research progresses, providing a community platform that establishes a common language to sharpen discussion across the Earth Science community.

No abstract

No abstract

We propose a model for increasing liquid saturation in a granular packing, which can account for liquid redistribution at saturation levels beyond the well-studied capillary bridge regime. The model is capable of resolving and combining capillary bridges, menisci, and fully saturated pores to form local liquid clusters of any shape. They can exchange volume due to the local Laplace pressure gradient via a liquid film on the surfaces of grains. Local instabilities such as Haines jumps trigger the discontinuous evolution of the liquid front. The applicability of the model is demonstrated and compared to benchmark experiments on the level of individual liquid structures as well as on larger systems.

In computational geodynamics, the Finite Element (FE) method is frequently used. The method is attractive as it easily allows employment of body-fitted deformable meshes and a true free surface boundary condition. However, when a Lagrangian mesh is used, remeshing becomes necessary at large strains to avoid numerical inaccuracies (or even wrong results) due to severely distorted elements. For this reason, the FE method is oftentimes combined with the particle-in-cell (PIC) method, where particles are introduced which track history variables and store constitutive information. This implies that the respective material properties have to be interpolated from the particles to the integration points of the finite elements. In numerical geodynamics, material parameters (in particular the viscosity) usually vary over a large range. This may be due to strongly temperature-dependent rheologies (which result in large but smooth viscosity variations) or material interfaces (which result in viscosity jumps). Here, we analyze the accuracy and convergence properties of velocity and pressure of the hybrid FE-PIC method in the presence of large viscosity variations. Standard interpolation schemes (arithmetic and harmonic) are compared to a more sophisticated interpolation scheme which is based on linear least squares interpolation for two types of elements (Q 1 P 0 and Q 2 P À1 ). In the case of a smooth viscosity field, the accuracy and convergence is significantly improved by the new interpolation scheme. In the presence of viscosity jumps, the order of accuracy is strongly decreased.

The presence of water in minerals generally alters their physical properties. Ringwoodite is the most abundant phase in the lowermost mantle transition zone and can host up to 1.5–2 wt% water. We studied high‐pressure lattice thermal conductivity of dry and hydrous ringwoodite by combining diamond‐anvil cell experiments with ultrafast optics. The incorporation of 1.73 wt% water substantially reduces the ringwoodite thermal conductivity by more than 40% at mantle transition zone pressures. We further parameterized the ringwoodite thermal conductivity as a function of pressure and water content to explore the large‐scale consequences of a reduced thermal conductivity on a slab's thermal evolution. Using a simple 1‐D heat diffusion model, we showed that the presence of hydrous ringwoodite in the slab significantly delays decomposition of dense hydrous magnesium silicates, enabling them to reach the lower mantle. Our results impact the potential route and balance of water cycle in the lower mantle.

No abstract

No abstract