Tectonic Motions Research Articles

Robust Imaging of Fault Slip Rates in the Walker Lane and Western Great Basin From GPS Data Using a Multi‐Block Model Approach

AbstractThe Walker Lane (WL) in the western Great Basin (GB) is an active plate boundary system accommodating 10%–20% of the relative tectonic motion between the Pacific and North American plates. Its neotectonic framework is structurally complex, having hundreds of faults with various strikes, rakes, and crustal blocks with vertical axis rotation. Faults slip rates are key parameters needed to quantify seismic hazard in such tectonically active plate boundaries but modeling them in complex areas like the WL and GB is challenging. We present a new modeling strategy for estimating fault slip rates in complex zones of active crustal deformation using data from GPS networks. The technique does not rely on prior estimates of slip rates from geologic studies, and only uses data on the surface trace location, dip, and rake. The iterative framework generates large numbers of block models algorithmically from the fault database to obtain many estimates of slip rates for each fault. This reduces bias from subjective choices about how discontinuous faults connect and interact to accommodate strain. Each model iteration differs slightly in block boundary configuration, but all models honor geodetic and fault data, regularization, and are kinematically self‐consistent. The approach provides several advantages over bespoke models, including insensitivity to outlier data, realistic uncertainties, explicit mapping of off‐fault deformation, and slip rates that are more objective and independent of geologic slip rates. Comparisons to the U.S. National Seismic Hazard Model indicate that ∼80% of our geodetic slip rates agree with their geologic slip rates to within uncertainties.

Python Software Tool for Diagnostics of the Global Navigation Satellite System Station (PS-NETM)–Reviewing the New Global Navigation Satellite System Time Series Analysis Tool

The time series of GNSS coordinates contain signals caused by the age-related movement of tectonic plates, the deformation of the Earth’s surface, as well as errors at different time scales from sub-daily tidal deformation to the long-term deformation of the surface load. Depending on the nature of the signal, specific approaches are used for both the visual interpretation and pre-processing of time series and their statistical analysis. However, none of the present software analyzes the nature of the residual errors but assumes their random nature and obedience to the classical normal distribution. One of the methods for analyzing the time series of coordinates with residual, unaccounted-for systematic errors is the non-classical error theory of measurements. The result of this work is a developed software solution for analyzing the time series of GNSS coordinates to test their normality, or in other words, to test whether a particular GNSS station is subject to the influence of small, unaccounted-for errors. Conclusions: After testing our software on four reference stations in Europe, we concluded that none of the chosen stations followed the normal law of distribution; thus, it is vital to perform such tests before conducting any experiments on the time series from reference stations.

A Review of the Contribution of Satellite Altimetry and Tide Gauge Data to Evaluate Sea Level Trends in the Adriatic Sea within a Mediterranean and Global Context

The relatively new sea level satellite altimetry and secular coastal tide gauge data made the reconstruction of sea levels on regional and global scales possible about one century back. Due to better estimations of the Earth’s crustal, glacial, tectonic, and other possible motion biases in tide gauge data, some additional improvements can be expected in sea level reconstructions, analysis, and predictions. A more detailed review of published sea level-related results was conducted for the Eastern Adriatic coast, including the operation of the tide gauge network and data processing, crustal movement estimations, and the establishment of a new reference height system in Croatia, based on five tide gauge sea level data. It was shown that sea level variation and trend-related indicators are spatially homogeneous, especially on a sub-Adriatic scale. The regional Adriatic Sea mean sea level rise rate of +2.6 mm/year for the satellite altimetry era (1993–2019) is less than the global mean sea level (GMSL) rise rate of +3.3 mm/year for the period of 1993–2022. Several empirical methods for GMSL projections and expected IPCC (Intergovernmental Panel on Climate Change) assessments until the end of the 21st century are considered.

Primena vremenski promenljivih graničnih uslova za brzinu konvergencije ploča pri numeričkom modelovanju obdukcije ofiolita

2D numerical geodynamic modelling applied in this paper gives a better insight into the process of obduction of two distinct ophiolite belts during the closure of a single ocean. Closure of the ocean is achieved by imposed convergence of the tectonic plates. In the created numerical themomechanical models, the influence of the variable velocity of the tectonic plate motions is analyzed. Convergence rate variation is modelled with the goal of achieving obduction and the permanent emplacement of ophiolites via complex evolution of the plate velocities. The basis of the model is the closure of an ocean through a flat intra-oceanic subduction, during which the plates change intensity and direction of their motion. Divergent phase of the evolution leads to the short-lived spreading of the new oceanic domain in from which, after the final closure, one of the two ophiolite belts originates.

Temporal and spatial distribution and genesis of primary indium deposits around the world

Because rare dispersed elements including indium play an increasingly extensive and indispensable role in the field of modern high technology, the exploration, development and research of its resources are in full swing. The abundance of indium in the earth's crust is extremely low, and its original independent deposits have not been found so far. Indium is mainly parasitic in bulk non-ferrous metal deposits in the form of associated by-products. Because of its origin, the genetic types of indium deposits are as diverse as those non-ferrous metal deposits it parasitizes. The distribution of indium in geological time and earth space is uneven. The formation of indium ore deposit is closely related to tectonic-magmatic activity. It is the product of plate tectonic movement, continental rifts and mantle plumes or mantle hotspots in the interior of ancient continents. The emergence of indium ore deposit is inseparable from magmatic activity, especially granite magmatic activity. Ore-forming materials mainly come from the upper mantle and/or deep crust. Indium deposits formed under high temperature conditions have enrichment mechanisms. In addition to the enrichment factors of traditional metal deposits, the enrichment of nano-scale indium minerals through the nano effect is also a very important enrichment method. The indium deposits formed at the end of the differentiation and evolution of granite magma may not be enriched in one go, but is the result of multiple enrichments.

ITRF2020 Plate Motion Model

AbstractA tectonic Plate Motion Model (PMM) is essential for geodetic applications, while contributing to the understanding of geodynamic processes affecting the Earth's surface. We introduce a PMM derived from the horizontal velocities of 518 sites extracted from the ITRF2020 solution. These sites were chosen away from plate boundaries, Glacial Isostatic Adjustment regions, and other deforming zones. Unlike the ITRF2014‐PMM, which showed no significant Origin Rate Bias (ORB), velocities used to determine the ITRF2020‐PMM exhibit a statistically significant ORB (0.74 ± 0.09 mm/yr along the Z‐component). Users are advised to add the estimated ORB to the horizontal velocities predicted by the ITRF2020‐PMM rotation poles for full consistency with the ITRF2020. However, the predicted vertical velocities resulting from the addition of the ORB should be discarded. The overall precision with which the ITRF2020 velocity field is represented by the rigid ITRF2020‐PMM is at the level of 0.25 mm/yr WRMS.

Mathematical fundamentals of spherical kinematics of plate tectonics in terms of quaternions

To be a quantitative and testable tectonic model, plate tectonics requires spherical geometry and spherical kinematics in terms of finite rotations conveniently parametrized by their angle and axis and described by unit quaternions. In treatises on “Plate Tectonics” infinitesimal, instantaneous, and finite rotations, absolute and relative rotations are said to be applied to model the motion of tectonic plates. Even though these terms are strictly defined in mathematics, they are often casually used in geosciences. Here, their definitions are recalled and clarified as well as the terms rotation, orientation, and location on the sphere. For instance, infinitesimal rotations refer to a mathematical limit, when the angle of rotation tends to zero. Their rules do not apply to finite rotations, no matter how small their finite angles of rotation are. Mathematical approaches applying appropriate and feasible assumptions to model spherical motion of tectonic plates over geological times of hundreds of millions of years are derived including (i) sequences of incremental finite rotations, (ii) sequences of accumulating successive concatenations of finite rotations, and (iii) continuous rotations in terms of fully transient quaternions. The incremental and the accumulating approaches provide complementary views. While the relative Euler pole appears to migrate in the latter, it appears fixed in the former. Path, mean, and instantaneous velocity of the migrating Euler pole are derived as well as the angular and trajectoral velocity of the rotational motion about it. The approaches are illustrated by a geological example with actual data and a numerical yet geologically inspired example with artificial data. The former revisits the three‐plate scenario with stationary axes of two “absolute” rotations implying transient “relative” rotations about a migrating Euler pole and employs a proper plate‐circuit argument to determine them numerically without resorting to approximations. The latter applies an involved interplay of incremental and accumulating modeling inducing split–join cycles to approximate sinusoidal trajectories as reported to record plates' motion during the Gondwana breakup.

Generation of Evolving Plate Boundaries and Toroidal Flow From Visco‐Plastic Damage‐Rheology Mantle Convection and Continents

AbstractEarth's style of planetary heat transport is characterized by plate tectonics which requires rock strength to be reduced plastically in order to break an otherwise stagnant lithospheric lid, and for rocks to have a memory of past deformation to account for strain localization and the hysteresis implied by geological sutures. Here, we explore ∼107 Rayleigh number, visco‐plastic, 3‐D global mantle convection with damage. We show that oceanic lithosphere‐only models generate strong toroidal‐poloidal power ratios and features such as a mix of long‐wavelength tectonic motions and smaller‐scale, back‐arc tectonics driven by downwellings. Undulating divergent plate boundaries can evolve to form overlapping spreading centers and microplates, promoted and perhaps stabilized by the effects of damage with long memory. The inclusion of continental rafts enhances heat flux variability and toroidal flow, including net rotation of the lithosphere, to a level seen in plate reconstructions for the Cenozoic. Both the super‐continental cycle and local rheological descriptions affect heat transport and tectonic deformation across a range of scales, and we showcase both general tectonic dynamics and regionally applied continental breakup scenarios. Our work points toward avenues for renewed analysis of the typical, mean behavior as well as the evolution of fluctuations in geological and model plate boundary evolution scenarios.

On the distribution of tangential mass forces in the Earth's lithosphere

The study aims to determine and interpret the distribution of the global tangential mass force (TMF) vector field by azimuthal orientation and intensity. Using cluster and correlation analysis, we compared the direction of the TMF vector field with the direction of movement of permanent GNSS stations and the direction of movement of the GSRM model continental velocities from the Global Strain Rate Map Project. Methodology. The author continues their study of additional planetary stresses in the lithosphere caused by distributed mass forces. The forces in question may be linked to the repositioning of the Earth's lithosphere, which can create stresses aimed at aligning the distribution of lithospheric masses with the geoid's figure. This repositioning happens through the mechanism of gravitational forces and the principle of minimum potential energy. The presence of a deviation of the plum line from the normal to the surface of the solid Earth determines the appearance of TMF acting in the upper shell of the Earth. It is proposed to calculate the amplitudes and directions of the vectors of such TMF based on data regarding the difference in the parameters of two global ellipsoids that approximate the physical surface of the lithosphere and the geoid. Originality. For the modern era, the value of the angle of rotation between the smallest axis of the ellipsoid approximating the surface of the lithosphere and the axis of rotation of the Earth is 2.6°. The distribution of the TMF vector field is consistent with the contours of the continents, i.e., the arrows of the vectors indicate the directions of lateral movement of tectonic plates and the movement of continents during the Earth's evolution. As a result of the change in the orientation of the ellipsoid describing the lithosphere, an updated field of potential horizontal forces is formed, which, by the conservation of the momentum of motion, move lithospheric masses and generate stresses and deformations in the lithospheric shell. Since the TMF has different directions and intensities, a cluster analysis of the TMF distribution was performed. It revealed certain regularities in the distribution of these parameters. We also compared the directions of the TMF vector field with the directions of movement of permanent GNSS stations and the directions of movement of model velocities of the continents of the GSRM (digital model of the tensor field of the global velocity gradient). Scientific novelty. The study detailed the peculiarities of the connection between the directions of the TMF vector field, the directions of movement of permanent GNSS stations, and the ones of the model velocities of the GSRM continents. Studies of the TMF, which arise as a result of the reorientation of the thin solid shell of our planet, have shown that a deformation field of shear is formed on its surface. In our opinion, this is one of the likely factors of the process that triggers global movements of lithospheric blocks. As a result, the shape of the lithosphere is transformed, which is characterized by a change in the size of the axes of the ellipsoids describing the surface of the lithosphere and their orientation. Practical significance. The research results make it possible to more reliably interpret the peculiarities of the TMF distribution. These forces can trigger mechanisms for discharging accumulated stresses, which is important for studying seismicity.

Stability analysis of reference station and compensation for monitoring stations in GNSS landslide monitoring

The Real-Time Kinematic (RTK) positioning method of the Global Navigation Satellite System (GNSS) has been widely used for landslide monitoring. The stability of its reference station is crucial to obtain accurate and reliable monitoring results. Unstable reference stations due to the geological environment and human activities are difficult to detect and in practical applications often ignored. As a result, it affects the positioning solutions and subsequently the interpretation and detection of landslide motions, which must be addressed in GNSS landslide monitoring. To solve this problem, we propose using the Precise Point Positioning (PPP) technique to analyze the stability of the reference station by verifying its position. The deformations of the monitoring stations are then compensated. First, the reference station coordinates are obtained by the PPP technique and tectonic motion is considered in data processing. The change or breakout of the reference station position is then determined using a cumulative sum control chart method. Finally, each monitoring station’s displacements are compensated according to the displacements of the reference station. According to the results of the Tengqing landslide experiment, the PPP technique can be used in GNSS landslide monitoring to analyze the stability of reference stations. With PPP, millimeter-level accuracy for the coordinates of reference stations is achieved. Compared to the traditional deformation series, the compensated displacement series more reliably reflects the landslide motions. This study will increase the reliability of monitoring results and contribute to implementing GNSS in monitoring landslides.

Main Ordovician–Silurian events in the South Ural and event stratigraphy regional tasks

The main regional abiotic events of the Ordovician and Silurian of the Southern Urals are manifested in the restructuring of sedimentation in the paleobasin. Eustatic fluctuations and tectonic motions primarily cause changes in the lithological composition of strata. Geochemical events influence the isotopic composition of elements.
 Biotic events affect biodiversity in micro- and macrofauna complexes. The largest of the five Phanerozoic extinctions was the global extinction event at the Ordovician–Silurian boundary (LOME). The Ordovician biotic diversification event (GOBE) is divided into three phases. The Darrivillian level of manifestation of this event is noted in the Southern Urals.
 The location of the global biotic events GOBE, LOME, and the geochemical event HICE in the sections of the Southern Urals is substantiated. Based on the analysis of the connection between biotic and abiotic events in the Ordovician-Silurian deposits of the Southern Urals, the significant role of conodonts, graptolites, chitinosoans and their correlation potential were revealed. Insufficient knowledge of the sections makes it possible to correlate regional events with certain assumptions.
 In order to carry out global event correlations, in addition to the comprehensive biostratigraphic characteristics of the Ordovician and Silurian deposits of the Southern Urals, geochemical studies of a number of sections are necessary. The most promising sections for identifying event levels are the sections of the western slope of the Zilair synclinorium: Nabiullino, Maksyutovo, Gremyachiy, Kaga-V.Avzyan, Sermenevo.

Nano‐bioconvective anisotropic slip flow in anisotropic porous medium with Coriolis force effects

AbstractThe present paper examines the impact of the Coriolis force and anisotropic Darcian porous medium on the nonsteady convective flow of bio‐nanofluid along a three‐dimensional stretchable surface. The problem has applications in extrusion processes and the geographical movement of tectonic plates under a moveable ocean. Emphasis is placed on the combined effects of anisotropic porous medium, anisotropic momentum slip and temperature, and microorganism slips at the wall. The pertaining equations are condensed to a set of similarity equations before being simulated using the Keller‐box method, an efficient numerical technique. The study reveals that the friction factor along the y‐axis decreases as the rotation parameter increases. Moreover, the heat and nanoparticle volume fraction movement gradient decrease along the y‐axis when the Darcy number and stretching rate increase. However, these rates increase with a higher velocity slip. In addition, the density of microorganisms in the bio‐nanofluid decreases when there are higher levels of unsteadiness, rotation, and Darcy numbers, as well as with increased bioconvection parameters.

Modes of Mantle Convection, Their Stability, and What Controls Their Existence

AbstractThe motion of Earth's tectonic plates is the surface expression of mantle convection beneath. Analytical convection models have attempted to relate observables, such as plate velocities and surface heat flow, with the thermo‐mechanical state of the mantle, and remain deeply influential in global geophysics. While such models tend to focus on describing the mantle's behavior today, there is evidence that suggests they may not be the best description for an earlier, hotter mantle. Early in Earth's history, higher temperatures may have led to different convective regimes that include active‐lid (today's form of convection), sluggish‐lid, and stagnant‐lid convection. In this study, we adopt and extend the analytical theory laid out by Crowley & O'Connell (2012, https://doi.org/10.1111/j.1365-246X.2011.05254.x) that self‐consistently characterizes the first two of these convection regimes. For a given thermo‐mechanical state, the theory predicts one to multiple solutions that each represent distinct modes of mantle convection. Here, we derive new scaling laws that connect these modes to the mantle's Rayleigh number and identify a new fundamental, dimensionless number, which we term the O’Connell number (named after Richard “Rick” J. O'Connell (1941–2015)), that describes a plate's resistance in context of the overall convection of the system. In addition, we identify two key bifurcations that bracket the Rayleigh‐O'Connell number space within which multiple solutions may exist. Finally, we remove the assumption of steady‐state in their original framework in order to perform a linear stability analysis to characterize the relative stability of each convection mode.

Biogeographic influences on the evolution and historical dispersal of the Australo‐Pacific Dacini fruit flies (Tephritidae: Dacinae)

AbstractFruit flies (Tephritidae: Dacini) are a frugivorous insect group that exhibit high endemic diversity in the rainforests of Australia and the western Pacific. In this region, biogeography has been influenced by tectonic plate movements and cycles of isolation and re‐connection of landmasses and rainforest habitats during glacial periods. However, how such factors have influenced the speciation and historical dispersal of the regional Dacini is largely unknown. To address this, we use a dated phylogeny to reconstruct the biogeographical history of the tribe. We found the Dacini radiated eastward into the Pacific islands largely from sources in New Guinea. We also found evidence for historical dispersal from both Australia and New Guinea into New Caledonia, a pathway unique to this island compared with neighbouring islands. There was also evidence for multiple, bidirectional dispersal events between Papua New Guinea and Australia, likely facilitated by the cyclically exposed Torres Strait land bridge. Cape York in far northern Australia was likely the only entry point for species dispersing into Australia; there was no evidence for entry of flies into Australia directly from West Papua or Wallacea. Several lineages radiated after entering Australia, such as members of the Bactrocera dorsalis species group. Within Australia, speciation was not associated with the biogeographic barriers known to have impacted other rainforest fauna in eastern Australia. Overall, we demonstrate that isolation between islands and large landmasses is important in the evolution of the Australo‐Pacific Dacini, but the reason for their extensive radiation within Australia and Papua New Guinea remains unclear.

Evidence for enhanced aridification since 13 Ma in the Qom back-arc basin, Central Iran

The Iranian Plateau was once occupied by the Neotethys Sea, which connected with the Indian Ocean in the southeast and the proto-Mediterranean in the west during most of the Cenozoic era. However, the Neotethys Sea has undergone a stepwise retreat, and the Tethyan Seaway also closed in the Miocene driven by both Arabia-Eurasia collision and global eustatic sea-level drop. Changes in the land-sea distribution and uplift of the Zagros Mountain belt must have affected paleoclimatic changes in Iran. To date, a high-resolution paleoclimatic record is lacking for Iran. In this paper, we present new environmental magnetic and color parameter data for terrestrial sediments in the Qom back-arc basin, Central Iran from the middle to late Miocene (17–11.5 Ma). All of our parameters indicate an abrupt climatic change at ∼13 Ma. We suggest that the less reddish color and a shift in environmental magnetic parameters after ∼13 Ma implies weakened chemical weathering under a colder/drier climate. Enhanced aridification was mainly controlled by closure of the Tethyan Seaway and shrinkage of the Neotethys Sea after ∼13 Ma, which greatly reduced moisture transport to the downwind Central Iran by westerlies, whereas the rain shallow effect of the initially uplifted Zagros Mountain belt and reduced seawater evaporation related to lower sea-surface temperatures after the Mid-Miocene Climatic Transition only played minor roles. Therefore, our results provide a window to discuss relationships between plate tectonic movements, land/sea distribution, and climatic change in a tectonically active region during the late Cenozoic.

High‐Resolution Mantle Flow Models Reveal Importance of Plate Boundary Geometry and Slab Pull Forces on Generating Tectonic Plate Motions

AbstractMantle 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 × 1017 and 5 × 1018 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.

Online data service for geologic formations (Lexicons) of China, India, Vietnam and Thailand with one‐click visualizations onto East Asia plate reconstructions

AbstractPaleogeography is the merger of sediment and volcanic facies, depositional settings, tectonic plate movements, topography, climate patterns and ecosystems through time. The construction of paleogeographic maps on tectonic plate reconstruction models requires a team effort to compile databases, data sharing standards and map projection methods. Two goals of the Paleogeography Working Group of the International Union of Geological Sciences (IUGS) program for Deep‐Time Digital Earth (DDE) are: (1) to interlink online national lexicons for all sedimentary and volcanic formations, and develop online ones for nations that currently lack these; (2) to target specific regions and intervals for testing/showcasing paleogeography output from the merger of these databases. Following those goals, we developed and applied new cloud‐based lexicon data services and interactive visualization techniques to regions in East Asia. This has been a successful collaboration among computer engineers and plate modellers and has involved stratigraphers in India (ONGC), China (Chengdu Univ. Tech., and Chinese Acad. Geol. Sci.), United States (Paleomap Project, and Purdue Univ.), Australia (GPlates visualization team) and Vietnam (Vietnam Nat. Univ.). Independent online lexicons with map‐based and stratigraphy‐based user interfaces have been developed (as of the date of this submission in March 2022) for all Proterozoic to Quaternary formations on the Indian Plate (over 800) and Vietnam (over 200), the majority of the Devonian through Neogene of China (ca. 2000) and partially for Thailand. A multi‐database search system returns all geologic formations of a desired geologic time from these four independent databases. With one click, users can plot the regional extent of one or of all of those regional formations on different plate reconstruction models of that desired age, and these polygons are filled with the lithologic facies pattern.

Probabilistic Evaluation of Seismic Performance of Steel Buildings with Torsional Irregularities in Plan and Soft Story under Mainshock-Aftershock Sequence

Due to the tectonic plate movement, numerous aftershocks may occur following an earthquake. There are several real mainshock-aftershock ground motion records in which the peak ground acceleration of aftershock is greater than that of the mainshock. Furthermore, irregularities such as soft story and torsional irregularity may change the structural behavior under the influence of earthquakes. This study investigates the fragility curves associated with the three-, five-, and eight-story models with steel moment-resisting frames under the main earthquake and mainshock-aftershock sequence to probabilistically evaluate the aftershock impacts on the steel structures with irregularities in the plan and height. The seismic fragility curves were calculated for four damage levels by selecting the relative displacement capacity at seismic performance levels of slight, moderate, extensive, and complete damage from the US HAZUS code. The analysis was done using structural reliability relationships and incremental dynamic analysis with the OpenSees software platform. The analytical results showed that the structures with soft-story and torsional irregularity were more vulnerable with an increasing number of stories. Also, aftershocks were found to have a more destructive effect on low-rise models. The probability of structural collapse in a given peak ground acceleration for a damaged state due to the mainshock-aftershock sequence is higher than that of the mainshock sequence.

Activated volcanism of Mount Fuji by the 2011 Japanese large earthquakes

The relation between earthquakes and volcanic eruptions, each of which is manifested by large-scale tectonic plate and mantle motions, has been widely discussed. Mount Fuji, in Japan, last erupted in 1707, paired with a magnitude (M)-9-class earthquake 49 days prior. Motivated by this pairing, previous studies investigated its effect on Mount Fuji after both the 2011 M9 Tohoku megaquake and a triggered M5.9 Shizuoka earthquake 4 days later at the foot of the volcano, but reported no potential to erupt. More than 300 years have already passed since the 1707 eruption, and even though consequences to society caused by the next eruption are already being considered, the implications for future volcanism remain uncertain. This study shows how volcanic low-frequency earthquakes (LFEs) in the deep part of the volcano revealed unrecognized activation after the Shizuoka earthquake. Our analyses also show that despite an increase in the rate of occurrence of LFEs, these did not return to pre-earthquake levels, indicating a change in the magma system. Our results demonstrate that the volcanism of Mount Fuji was reactivated by the Shizuoka earthquake, implying that this volcano is sufficiently sensitive to external events that are considered to be enough to trigger eruptions.

The complex Rodrigues triple junction migration since ca. 8 Ma: A response to episodic Amsterdam–St. Paul hotspot tail capture by the Southeast Indian Ridge?

AbstractThe mutual intersection of three plate boundaries, so-called triple junctions, has long been recognized as crucial boundaries for unraveling the spatiotemporal motion of tectonic plates. Yet, the dynamic and tectonic processes ruling their migration remain enigmatic. At the Rodrigues triple junction, the Southwest Indian Ridge lengthens northeastward in response to the unsteady linearity of the Southeast Indian Ridge–Central Indian Ridge dueling ridge system. This mechanism is episodically compensated by the southeastward propagation of the Central Indian Ridge and alternating stages of recession and progression of the Southeast Indian Ridge. This has led to an apparent length constancy of the first Southeast Indian Ridge segment over the past ca. 8 Ma. These multiple ridge propagation episodes, from northwest, southwest, and southeast, result in a net northeastward migration of the Rodrigues triple junction. Here, we suggest that its migration since ca. 8 Ma is a consequence of short-term readjustments of its plate boundaries induced by transitory motion changes of the Capricorn plate, driven by episodic push forces exerted from ephemeral captures of the Amsterdam–St. Paul plume tail by the Southeast Indian Ridge at the southeastern part of the Capricorn plate (77.3°E–78.6°E).