Progressive weakening within the overriding plate during dual inward dipping subduction

Dual inward dipping subduction often produces complex deformation patterns in the overriding plate. However, the geodynamic process of how dual inward dipping subduction relates to this deformation is still poorly understood. Here we apply a composite viscosity, dependent on multiple parameters, e.g., temperature, pressure, strain rate etc., in 2-D thermo-mechanical numerical modelling to investigate how dual inward dipping subduction modifies the rheological structure of the overriding plate. Three variables are investigated to understand what controls the maximum degree of weakening. We find that the initial length and thickness of the overriding plate are negatively correlated with the magnitude of viscosity reduction. While the initial thickness of the subducting plate positively relates to the magnitude of viscosity reduction. The progressive weakening can result in a variety of stretching states ranging from 1) little or no lithosphere thinning and extension, to 2) limited thermal lithosphere thinning, and 3) localised rifting followed by spreading extension. Compared with single sided subduction, dual inward dipping subduction further reduces the magnitude of viscosity of the overriding plate. It does this by creating a dynamic fixed boundary condition for the overriding plate and forming a stronger upwelling mantle flow underlying the overriding plate. Three types of feedback weakening cycles are recognised, among which the strain rate weakening mechanism plays the dominant role in lowering the viscosity of the overriding plate throughout the simulation. Strain rate weakening is also a precondition for initiating thermal weakening, strain localisation and lithosphere thinning.

Cooling of the continental plate during flat-slab subduction

Most flat-slab subduction regions are marked by an absence of arc volcanism, which is consistent with closure of the hot mantle wedge as the subducting plate flattens below the continent. Farther inland, low surface heat flow is observed, which is generally attributed to cooling of the continent by the underlying flat slab. However, modern flat slabs have only been in place for <20 Ma, and it is unclear whether there has been sufficient time for cooling to occur. We use numerical models to assess temporal variations in continental thermal structure during flat-slab subduction. Our models show that the flat slab leads to continental cooling on timescales of tens of millions of years. Cool slab temperatures must diffuse through the continental lithosphere, resulting in a delay between slab emplacement and surface cooling. Therefore, the timescales primarily depend on the flat-slab depth with shallower slabs resulting in shorter timescales. The magnitude of cooling increases for a shallow or long-lived flat slab, old subducting plate, and fast convergence rates. For regions with flat slabs at 45–70 km depth (e.g., Mexico and Peru), shallow continental cooling initiates 5–10 Ma after slab emplacement, and low surface heat flow in these regions is largely explained by the presence of the flat slab. However, for the Pampean region in Chile, with an ~100-km-deep slab, our models predict that conductive cooling has not yet affected the surface heat flow. The low heat flow observed requires additional processes such as advective cooling from the infiltration of fluids released through dehydration of the flat slab.

Temperature distribution in a subduction zone

A model of the thermal and hydrodynamic structure of the subduction zone is proposed. This model includes free convection flows in the asthenospheric layer and layer C (mantle transition zone). Temperature profiles in the subducting lithospheric plate, as well as in the continental limb of the subduction zone, are presented. The heat flux due to friction at the contact between the subducting plate and the continental limb significantly affects the heat transfer and, consequently, the temperature field formation in the subduction zone. The temperature level in the crustal layer of the submerging plate implies that there is no melting in the crustal layer.

Lithosphere–Atmosphere–Ionosphere Coupling Effects Based on Multiparameter Precursor Observations for February–March 2021 Earthquakes (M~7) in the Offshore of Tohoku Area of Japan

The purpose of this paper is to discuss the lithosphere–atmosphere–ionosphere coupling (LAIC) effects with the use of multiparameter precursor observations for two successive Japanese earthquakes (EQs) (with a magnitude of around 7) in February and March 2021, respectively, considering a seemingly significant difference in seismological and geological hypocenter conditions for those EQs. The second March EQ is very similar to the famous 2011 Tohoku EQ in the sense that those EQs took place at the seabed of the subducting plate, while the first February EQ happened within the subducting plate, not at the seabed. Multiparameter observation is a powerful tool for the study of the LAIC process, and we studied the following observables over a 3-month period (January to March): (i) ULF data (lithospheric radiation and ULF depression phenomenon); (ii) ULF/ELF atmospheric electromagnetic radiation; (iii) atmospheric gravity wave (AGW) activity in the stratosphere, extracted from satellite temperature data; (iv) subionospheric VLF/LF propagation data; and (v) GPS TECs (total electron contents). In contrast to our initial expectation of different responses of anomalies to the two EQs, we found no such conspicuous differences of electromagnetic anomalies between the two EQs, but showed quite similar anomaly responses for the two EQs. It is definite that atmospheric ULF/ELF radiation and ULF depression as lower ionospheric perturbation are most likely signatures of precursors to both EQs, and most importantly, all electromagnetic anomalies are concentrated in the period of about 1 week–9 days before the EQ to the EQ day. There seems to exist a chain of LAIC process (cause-and-effect relationship) for the first EQ, while all of the observed anomalies seem to occur nearly synchronously in time for the send EQ. Even though we tried to discuss possible LAIC channels, we cannot come to any definite conclusion about which coupling channel is plausible for each EQ.

The effect of sediments on the dynamics and accretionary style of subduction margins

Subduction zones represent the only major pathway by which continental material can be returned to the Earth's mantle. Constraining the sediments mass flux through subduction zones is important to the understanding of both petrogenesis of continental crust, and the recycling of volatiles and continental material back into the mantle over long periods of geologic time. When sediments are considered, convergent margins appear to fall into one of two classes: accretionary and erosive. Accretionary margins are dominated by accretion of thick piles of sediments (> 1 km) from the subducting plate, while tectonic erosion is favored in regions where the sedimentary cover is < 1 km. However, as data help define geometry of the global subduction system, the consequences of the two styles of margins on subduction dynamics remain poorly resolved. In this study, we run systematic 2-D numerical simulations of subduction to investigate how sediment fluxes influence subduction dynamics and plate coupling. We vary the thickness and viscosity of the sediment layer entering subduction, the thickness of the upper plate, and the driving velocity of the subducting plate (i.e., kinematic boundary conditions). Our results show three modes of subduction interface: a) Tectonic erosion margin (high viscosity sediment layer), b) Low angle accretionary wedge margin (low viscosity, thin sediment layer), and c) High angle accretionary wedge margin (low viscosity, thick sediment layer). We find that the properties of the sediment layer modulate the extent of viscous coupling at the interface between the subducting and overriding plates. When the viscous coupling is increased, an erosive style margin will be favored over an accretionary style. On the other hand, when the viscous coupling is reduced, sediments are scrapped-off the subducting slab to form an accretionary wedge. Diagnostic parameters are extracted automatically from numerical simulations to analyze the dynamics and differentiate between these modes of subduction margin. Models of tectonic erosion margins show small radii of curvature, slow convergence rates and thin subduction interfaces, while results of accretionary margins show large radii of curvature, faster convergence rates and dynamic accretionary wedges. These diagnostics parameters are then linked with observations of present-day subduction zones.

ON HEAT SOURCE IN SUBDUCTION ZONE

The subduction of an oceanic plate is studied as the motion of a high-viscosity Newtonian fluid. The subducting plate spreads along the 670-km depth boundary under the influence of oppositely directed horizontal forces. These forces are due to oppositely directed horizontal temperature gradients. We consider the flow structure and heat transfer in the layer that includes both the oceanic lithosphere and the crust and moves underneath a continent. The heat flow is estimated at the contact between the subducting plate and the surrounding mantle in the continental limb of the subduction zone. Our study results show that the crustal layer of the subducting plate can melt and a thermochemical plume can form at the 670-km boundary. Our model of a thermochemical plume in the subduction zone shows the following: (1) formation of a plume conduit in the crustal layer of the subducting plate; (2) formation of a primary magmatic chamber in the area wherein the melting rate equals the rate of subduction; (3) origination of a vertical plume conduit from the primary chamber melting through the continent; (4) plume eruption through the crustal layer to the surface, i.e. formation of a volcano. Our experiments are aimed to model the plume conduit melting in an inclined flat layer above a local heat source. The melt flow structure in the plume conduit is described. Laboratory modeling have revealed that the mechanisms of melt eruption from the plume conduit differ depending on whether a gas cushion is present or absent at the plume roof.

Subduction history of the Caribbean from upper-mantle seismic imaging and plate reconstruction

The margins of the Caribbean and associated hazards and resources have been shaped by a poorly understood history of subduction. Using new data, we improve teleseismic P-wave imaging of the eastern Caribbean upper mantle and compare identified subducted-plate fragments with trench locations predicted from plate reconstruction. This shows that material at 700–1200 km depth below South America derives from 90–115 Myr old westward subduction, initiated prior to Caribbean Large-Igneous-Province volcanism. At shallower depths, an accumulation of subducted material is attributed to Great Arc of the Caribbean subduction as it evolved over the past 70 Ma. We interpret gaps in these subducted-plate anomalies as: a plate window and tear along the subducted Proto-Caribbean ridge; tearing along subducted fracture zones, and subduction of a volatile-rich boundary between Proto-Caribbean and Atlantic domains. Phases of back-arc spreading and arc jumps correlate with changes in age, and hence buoyancy, of the subducting plate.