Search for: All records

Abstract Simulating present‐day solid Earth deformation and volatile cycling requires integrating diverse geophysical data sets and advanced numerical techniques to model complex multiphysics processes at high resolutions. Subduction zone modeling is particularly challenging due to the large geographic extent, localized deformation zones, and the strong feedbacks between reactive fluid transport and solid deformation. Here, we develop new workflows for simulating 3‐dimensional thermal‐mechanical subduction and patterns of volatile dehydration at convergent margins, adaptable to include reactive fluid transport. We apply these workflows to the Hikurangi margin, where recent geophysical investigations have offered unprecedented insight into the structure and processes coupling fluid transport and solid deformation across broad spatiotemporal scales. Geophysical data sets constraining the downgoing and overriding plate structure are collated with the Geodynamic World Builder, which provides the initial conditions for forward simulations using the open‐source geodynamic modeling software code ASPECT. We systematically examine how plate interface rheology and hydration of the downgoing plate and upper mantle influence Pacific–Australian convergence and seismic anisotropy. Models prescribing a plate boundary viscosity of Pa s best reproduce observed plate velocities, and changing the configuration of the Pacific–Australia plate boundary directly influences the modeled plate motions. Models considering hydrated olivine fabrics best reproduce observations of seismic anisotropy. Predicted patterns of slab dehydration and mantle melting correlate well with observations of seismic attenuation and arc volcanism. These results suggest that hydration‐related rheological heterogeneity and related fluid weakening may strongly influence slab dynamics. Future investigations integrating coupled fluid transport and global mantle flow will provide insight into the feedbacks between subduction dynamics, fluid pathways, and arc volcanism. 

Abstract Subduction zones are home to multiple geohazards driven by the evolution of the regional tectonics, including earthquakes, volcanic eruptions and landslides. Past evolution builds the present‐day structure of the margin, while the present‐day configuration of the system determines the state‐of‐stress in which individual hazardous events manifest. Regional simulations of subduction zones provide a tool to synthesize the tectonic history of a region and investigate how geologic features lead to variations in the state of stress across the subduction system. However, it is challenging to design regional models that provide a force‐balance that is consistent with the large‐scale motion of surrounding tectonic plates while also not over‐constraining the solution. Here, we present new models for the Cascadia subduction zone that meet these criteria and demonstrate how the motion of the subducting Juan de Fuca plate can be used to determine the along‐strike variations in the viscous (long‐term) coupling across the plate boundary. All successful models require lower viscous coupling in the northern section of the trench compared to the central and southern sections. However, due to uncertainties in the geometry of the Cascadia slab, we find that there is a trade‐off between along‐strike variation in viscous coupling and slab shape. Better constraints on the slab shape, and/or use of other observations are needed to resolve this trade‐off. The approach presented here provides a framework for further exploring how geologic features in the overriding plate and the properties of the plate boundary region affect the state‐of‐stress across this and other subduction zones. 

Abstract Mineral phase transitions can either hinder or accelerate mantle flow. In the present day, the formation of the bridgmanite + ferropericlase assemblage from ringwoodite at 660 km depth has been found to cause weak and intermittent layering of mantle convection. However, for the higher temperatures in Earth's past, different phase transitions could have controlled mantle dynamics. We investigate the potential changes in convection style during Earth's secular cooling using a new numerical technique that reformulates the energy conservation equation in terms of specific entropy instead of temperature. This approach enables us to accurately include the latent heat effect of phase transitions for mantle temperatures different from the average geotherm, and therefore fully incorporate the thermodynamic effects of realistic phase transitions in global‐scale mantle convection modeling. We set up 2‐D models with the geodynamics softwareAspect, using thermodynamic properties computed by HeFESTo, while applying a viscosity profile constrained by the geoid and mineral physics data and a visco‐plastic rheology to reproduce plate‐like behavior and Earth‐like subduction morphologies. Our model results reveal the layering of plumes induced by the wadsleyite to garnet (majorite) + ferropericlase endothermic transition (between 450 and 590 km depth and over the 2000–2500 K temperature range). They show that this phase transition causes a large‐scale and long‐lasting temperature elevation in a depth range of 500–650 km depth if the potential temperature of the mantle is higher than 1800 K, indicating that mantle convection may have been partially layered in Earth's early history. 

Abstract Determining the fate of subducted oceanic crust is critical for understanding material cycling through Earth’s deep interior and sources of mantle heterogeneity. A key control on the distribution of subducted slabs over long timescales is the bridgmanite to post-perovskite phase transition in the lowermost mantle, thought to cause rheological weakening. Using high-resolution computational models, we show that the ubiquitous presence of weak post-perovskite at the core-mantle boundary can facilitate or prevent the accumulation of basaltic oceanic crust, depending on the amount of weakening and the crustal thickness. Moderately weak post-perovskite ( ~ 10–100× weaker) facilitates segregation of crust from subducted slabs, increasing basalt accumulation in dense piles. Conversely, very weak post-perovskite (more than 100× weaker) promotes vigorous plumes that entrain more crustal material, decreasing basalt accumulation. Our results reconcile the contradicting conclusions of previous studies and provide insights into the accumulation of subducted crust in the lowermost mantle throughout Earth’s history. 

Abstract Large igneous provinces (LIPs) have been linked to both surface and deep mantle processes. During the formation, tenure and break-up of the supercontinent Pangaea, there is an increase in emplacement events for both continental and oceanic LIPs. There is currently no clear consensus on the origin of LIPs, but a hypothesis relates their formation to crustal emplacement of hot plume material originating in the deep mantle. The interaction of subducted slabs with the lowermost mantle thermal boundary and subsequent return flow is a key control on such plume generation. This mechanism has been explored for LIPs below the interior of a supercontinent (i.e. continental LIPs). However, a number of LIPs formed exterior to Pangaea (e.g. Ontong Java Plateau), with no consensus on their formation mechanism. Here, we consider the dynamics of supercontinent processes as predicted by numerical models of mantle convection and analyse whether circum-supercontinent subduction could generate both interior (continental) and exterior (oceanic) deep mantle plumes. Our numerical models show that subduction related to the supercontinent cycle can reproduce the location and timing of the Ontong Java Plateau, Caribbean LIP and potentially the Shatsky Rise by linking the origin of these LIPs to the return flow that generated deep mantle exterior plumes. 

Abstract Mantle convection models based on geophysical constraints have provided us with a basic understanding of the forces driving and resisting plate motions on Earth. However, existing studies computing the balance of underlying forces are contradicting, and the impact of plate boundary geometry on surface deformation remains unknown. We address these issues by developing global instantaneous 3‐D mantle convection models with a heterogeneous density and viscosity distribution and weak plate boundaries prescribed using different geometries. We find that the plate boundary geometry of the Global Earthquake Model (GEM, Pagani et al., 2018,https://doi.org/10.1177/8755293020931866), featuring open plate boundaries with discrete lithospheric‐depth weak zones in the oceans and distributed crustal faults within continents, achieves the best fit to the observed GPS data with a directional correlation of 95.1% and a global point‐wise velocity residual of 1.87 cm/year. A good fit also requires plate boundaries being 3 to 4 orders of magnitude weaker than the surrounding lithosphere and low asthenospheric viscosities between 5 × 10¹⁷and 5 × 10¹⁸ Pa s. Models without asthenospheric and lower mantle heterogeneities retain on average 30% and 70% of the plate speeds, respectively. Our results show that Earth's plate boundaries are not uniform and better described by more discrete plate boundaries within the oceans and distributed faults within continents. Furthermore, they emphasize the impact of plate boundary geometry on the direction and speed of plate motions and reaffirm the importance of slab pull in the uppermost mantle as a major plate driving force. 

Abstract Mantle plumes are thought to recycle material from the Earth's deep interior. One constraint on the nature and quantity of this recycled material comes from the observation of seismic discontinuities. The detection of the X‐discontinuity beneath Hawaii, interpreted as the coesite‐stishovite transition, requires the presence of at least 40% basalt. However, previous geodynamic models have predicted that plumes cannot carry more than 15%–20% of high‐density basaltic material. We propose this contradiction can be resolved by taking into account the length scale of chemical heterogeneities. While previous modeling studies assumed mechanical mixing on length scales smaller than the model resolution, we here model basaltic heterogeneities with length scales of 30–40 km, allowing for their segregation relative to the pyrolitic background plume material. Our models show that larger basalt fractions than previously thought possible—exceeding 40%—can temporarily accumulate within plumes at the depth of the X‐discontinuity. Two key mechanisms facilitate this process: (a) The random distribution of basaltic heterogeneities induces large temporal variations in the basalt fraction with cyclical highs and lows. (b) The high density contrast between basalt and pyrolite below the coesite‐stishovite transition causes ponding and accumulation of basalt within the rising plume at that depth. Because the statistical effect dominates, large values of 35%–40% basalt are only sustained temporarily. These results further constrain the chemical composition of the Hawaiian plume. Beyond that, they provide a geodynamic mechanism that explains the seismologic detection of the X‐discontinuity and highlights how recycled material is carried toward the surface. 

SUMMARY Phase transitions play an important role for the style of mantle convection. While observations and theory agree that a substantial fraction of subducted slabs and rising plumes can move through the whole mantle at present day conditions, this behaviour may have been different throughout Earth’s history. Higher temperatures, such as in the early Earth, cause different phase transitions to be dominant, and also reduce mantle viscosity, favouring a more layered style of convection induced by phase transitions. A period of layered mantle convection in Earth’s past would have significant implications for the secular evolution of the mantle temperature and the mixing of mantle heterogeneities. The transition from layered to whole mantle convection could lead to a period of mantle avalanches associated with a dramatic increase in magmatic activity. Consequently, it is important to accurately model the influence of phase transitions on mantle convection. However, existing numerical methods generally preclude modelling phase transitions that are only present in a particular range of pressures, temperatures or compositions, and they impose an artificial lower limit on the thickness of phase transitions. To overcome these limitations, we have developed a new numerical method that solves the energy equation for entropy instead of temperature. This technique allows for robust coupling between thermodynamic and geodynamic models and makes it possible to model realistically sharp phase transitions with a wide range of properties and dynamic effects on mantle processes. We demonstrate the utility of our method by applying it in regional and global convection models, investigating the effect of individual phase transitions in the Earth’s mantle with regard to their potential for layering flow. We find that the thickness of the phase transition has a bigger influence on the style of convection than previously thought: with all other parameters being the same, a thin phase transition can induce fully layered convection where a broad phase transition would lead to whole-mantle convection. Our application of the method to convection in the early Earth illustrates that endothermic phase transitions may have induced layering for higher mantle temperatures in the Earth’s past.