Initiation of subduction is rarely encountered in modern tectonic environments due to its ephemeral and destructive nature. We report the geological and geophysical evidence indicating a transitional phase from buckling to embryonic subduction along the eastern Korean margin. The transition appears to be caused by compressional reactivation of the strike-slip boundary between the continental (Korean Peninsula) and oceanic (Ulleung Basin) crusts since the Early Pliocene. Evidence for compressional reactivation includes (1) a west-dipping major thrust and coincident crustal buckling of the Ulleung Basin; (2) an east-west structural asymmetry inferred from the gravity anomaly and P-wave tomography; and (3) ongoing crustal uplift and high-angle faults along the eastern Korean margin. The juxtaposition of underthrusting and buckling of the crust in the Ulleung Basin, and its associated ubiquitous reverse faulting on the eastern Korean margin, imply the potential development of a new subduction system along the western margin of the East Sea (Japan Sea). We propose that the East Sea comprises two incipient subduction margins (i.e., the Korean and Japanese sides), which are now competing to reach a self-sustaining subduction stage.

SUMMARY Viscoplasticity has been considered to be a dominant element in causing the nucleation of shear instability leading to lithospheric weakening. Here, we propose that a simple contrast in shear moduli may be sufficient for explaining the fast timescale asymmetric shear instability in a bimaterial setting. Not much attention has been paid to heterogeneous elasticity in geodynamical modelling because it is dominant only for short timescales. Up to now, no studies have been made on asymmetric shear instability induced by elastic modulus contrast. Thermal–mechanical numerical simulations based on high‐resolution (from 0.4 km × 0.4 km to 0.2 km × 0.4 km meshes) finite‐element methods were performed to understand the effects of shear modulus contrast on inducing asymmetric instabilities. Strain‐rate and stress‐dependent rheology are used with a wide range of activation energy 0–850 kJ mol–1 for all models. Numerical results with enough shear modulus contrast show asymmetric shear instability, which is generated around the interface and then propagates across the interface. Two parts of the lithosphere with different shear moduli (stiff for higher and soft for lower shear modulus lithospheres), which are simply connected to each other without a pre‐defined weak zone, were compressed at a constant rate of 2 cm yr–1. Having different shear modulus is justified by chemical heterogeneity of geological minerals and their pressure–temperature dependence. To explore the dynamical effects generated by the contrast in the elastic modulus, the shear modulus of the soft lithosphere is fixed at 32 GPa, whereas that of stiff lithosphere is increased systematically from 32 up to 640 GPa. We also examined the role of activation energy (0–850 kJ mol–1) on the geometrical pattern and the initiation time of asymmetric shear localization. The shear modulus contrast has to be close to two for triggering asymmetric shear instability and is found to be by far a more important controlling factor in causing shear instability than activation energy of the creep law. The instability develops rapidly between 250 000 and 500 000 yr after deformation begins, and thermal weakening in the shear zone is greater, when a stronger shear modulus contrast is prescribed. Our work suggests that initiation of lithosphere‐scale asymmetric instability would be faster than previous considerations. Our finding stresses that naturally occurring shear modulus contrast has also important impact on many geological problems related to bimaterial instability.

Carbon dioxide emissions from dolomite decarbonation play an essential role in the weakening of carbonate faults by lowering the effective normal stress, which is thermally activated at temperatures above 600–700 °C. However, the mechanochemical effect of low-crystalline ultrafine fault gouge on the decarbonation and slip behavior of dolomite-bearing faults remains unclear. In this study, we obtained a series of artificial dolomite fault gouges with systematically varying particle sizes and dolomite crystallinities using a high-energy ball mill. The laboratory-scale pulverization of dolomite yielded MgO at temperatures below 50 °C, indicating that mechanical decarbonation without significant heating occurred due to the collapse of the crystalline structure, as revealed by X-ray diffraction and solid-state nuclear magnetic resonance results. Furthermore, the onset temperature of thermal decarbonation decreased to ~400 °C. Numerical modeling reproduced this two-stage decarbonation, where the pore pressure increased due to low-temperature thermal decarbonation, leading to slip weakening on the fault plane even at 400–500 °C; i.e., 200–300 °C lower than previously reported temperatures. Thus, the presence of small amounts of low-crystalline dolomite in a fault plane may lead to a severely reduced shear strength due to thermal decomposition at ~400 °C with a small slip weakening distance.

No abstract

Buoyant trench sediment generates an isostatic restoring force that opposes the slab pull of subducting plates; however, the influence of this force on subduction dynamics is poorly understood. Here, we performed three‐dimensional free subduction simulations adopting a variety of sediment distributions along a trench to investigate the correlation between trench motions and heterogeneous buoyant forces. Two endmembers, sediment‐rich and sediment‐starved centers, induced convex and concave trenches, respectively. The trench curvatures obtained from the natural subduction zones are well constrained by our models over wide magnitudes (1–9 km) and wavelengths (100–400 km) of deficient and excess sediments. Conversely, a uniform sediment distribution leads to an extremely narrow range of trench curvatures. These results imply that trench sediment contributes significantly to deviated curvatures along trench strikes. Finally, we observed stress localization at shallow slab hinges that lie beneath abundant sediments, clarifying the relationship between sediment thickness and subduction earthquakes.