Search

It has recently been demonstrated that the interaction of a mantle plume with sufficiently old oceanic lithosphere can initiate subduction. However, the existence of large lithospheric heterogeneities, such as a buoyant plateau, in proximity to a rising plume head may potentially hinder the formation of a new subduction zone. Here, we investigate this scenario by means of 3-D numerical thermomechanical modeling. We explore how plume-lithosphere interaction is affected by lithospheric age, relative location of plume head and plateau border, and the strength of the oceanic crust. Our numerical experiments suggest four different geodynamic regimes: (a) oceanic trench formation, (b) circular oceanic-plateau trench formation, (c) plateau trench formation, and (d) no trench formation. We show that regardless of the age and crustal strength of the oceanic lithosphere, subduction can initiate when the plume head is either below the plateau border or at a distance less than the plume radius from the plateau edge. Crustal heterogeneity facilitates subduction initiation of old oceanic lithosphere. High crustal strength hampers the formation of a new subduction zone when the plume head is located below a young lithosphere containing a thick and strong plateau. We suggest that plume-plateau interaction in the western margin of the Caribbean could have resulted in subduction initiation when the plume head impinged onto the oceanic lithosphere close to the border between plateau and oceanic crust.

Carbon is a key control on the surface chemistry and climate of Earth. Significant volumes of carbon are input to the oceans and atmosphere from deep Earth in the form of degassed CO2 and are returned to large carbon reservoirs in the mantle via subduction or burial. Different tectonic settings (e.g. volcanic arcs, mid-ocean ridges, and continental rifts) emit fluxes of CO2 that are temporally and spatially variable, and together they represent a first-order control on carbon outgassing from the deep Earth. Changes in the relative importance of different tectonic settings throughout Earth's history has therefore played a key role in balancing the deep carbon cycle on geological timescales. Over the past ten years the Deep Carbon Observatory has made enormous progress in constraining estimates of carbon outgassing flux at different tectonic settings. Using plate boundary evolution modelling and our understanding of present-day carbon fluxes, we develop time series of carbon fluxes into and out of the Earth’s interior through the past 200 million years. We highlight the increasing importance of carbonate-intersecting subduction zones over time to carbon outgassing, and the possible dominance of carbon outgassing at continental rift zones, which leads to maxima in outgassing at 130 and 15 Ma. To a first-order, carbon outgassing since 200 Ma may be net positive, averaging ~50 Mt C yr-1 more than the ingassing flux at subduction zones. Our net outgassing curve is poorly correlated with atmospheric CO2, implying that surface carbon cycling processes play a significant role in modulating carbon concentrations and/or there is a long-term crustal or lithospheric storage of carbon which modulates the outgassing flux. Our results highlight the large uncertainties that exist in reconstructing the corresponding in- and outgassing fluxes of carbon. Our synthesis summarizes our current understanding of fluxes at tectonic settings and their influence on atmospheric CO2, and provides a framework for future research into Earth’s deep carbon cycling, both today and in the past.

We evaluate the spatial and temporal evolution of Earth's long-wavelength surface dynamic topography since the Jurassic using a series of high-resolution global mantle convection models. These models are Earth-like in terms of convective vigour, thermal structure, surface heat-flux and the geographic distribution of heterogeneity. The models generate a degree-2-dominated spectrum of dynamic topography with negative amplitudes above subducted slabs (i.e. circum-Pacific regions and southern Eurasia) and positive amplitudes elsewhere (i.e. Africa, north-western Eurasia and the central Pacific). Model predictions are compared with published observations and subsidence patterns from well data, both globally and for the Australian and southern African regions. We find that our models reproduce the long-wavelength component of these observations, although observed smaller-scale variations are not reproduced. We subsequently define "geodynamic rules" for how different surface tectonic settings are affected by mantle processes: (i) locations in the vicinity of a subduction zone show large negative dynamic topography amplitudes; (ii) regions far away from convergent margins feature long-term positive dynamic topography; and (iii) rapid variations in dynamic support occur along the margins of overriding plates (e.g. the western US) and at points located on a plate that rapidly approaches a subduction zone (e.g. India and the Arabia Peninsula). Our models provide a predictive quantitative framework linking mantle convection with plate tectonics and sedimentary basin evolution, thus improving our understanding of how subduction and mantle convection affect the spatio-temporal evolution of basin architecture.

Localization processes in the viscous lower crust generate ductile shear zones over a broad range of scales affecting long‐term lithosphere deformation and the mechanical response of faults during the seismic cycle. Here we use centimeter‐scale numerical models in order to gain detailed insight into the processes involved in strain localization and rheological weakening in viscously deforming rocks. Our 2‐D Cartesian models are benchmarked to high‐temperature and high‐pressure torsion experiments on Carrara marble samples containing a single weak Solnhofen limestone inclusion. The models successfully reproduce bulk stress‐strain transients and final strain distributions observed in the experiments by applying a simple, first‐order softening law that mimics rheological weakening. We find that local stress concentrations forming at the inclusion tips initiate strain localization inside the host matrix. At the tip of the propagating shear zone, weakening occurs within a process zone, which expands with time from the inclusion tips toward the matrix. Rheological weakening is a precondition for shear zone localization, and the width of this shear zone is found to be controlled by the degree of softening. Introducing a second softening step at elevated strain, a high strain layer develops inside the localized shear zone, analogous to the formation of ultramylonite bands in mylonites. These results elucidate the transient evolution of stress and strain rate during inception and maturation of ductile shear zones.

Rock deformation at depths in the Earth’s crust is often localized in high temperature shear zones occurring at different scales in a variety of lithologies. The presence of material heterogeneities is known to trigger shear zone development, but the mechanisms controlling initiation and evolution of localization are not fully understood. To investigate the effect of loading conditions on shear zone nucleation along heterogeneities, we performed torsion experiments under constant twist rate (CTR) and constant torque (CT) conditions in a Paterson-type deformation apparatus. The sample assemblage consisted of cylindrical Carrara marble specimens containing a thin plate of Solnhofen limestone perpendicular to the cylinder’s longitudinal axis. Under experimental conditions (900 °C, 400 MPa confining pressure), samples were plastically deformed and limestone is about 9 times weaker than marble, acting as a weak inclusion in a strong matrix. CTR experiments were performed at maximum bulk shear strain rates of ~ 2*10-4s-1, yielding peak shear stresses of ~ 20 MPa. CT tests were conducted at shear stresses of ~ 20 MPa resulting in bulk shear strain rates of 1-4*10-4s-1. Experiments were terminated at maximum bulk shear strains of ~ 0.3 and 1.0.Strain was localized within the Carrara marble in front of the inclusion in an area of strongly deformed grains and intense grain size reduction. Locally, evidences for coexisting brittle deformation are also observed regardless of the imposed loading conditions. The local shear strain at the inclusion tipis up to 30 times higher than the strain in the adjacent host rock, rapidly dropping to 5times higher at larger distance from the inclusion. At both investigated bulk strains, the evolution of microstructural and textural parameters is independent of loading conditions. Ourresults suggest that loading conditions do not significantly affect material heterogeneity-induced strain localization during its nucleation and transient stages.