Large Slip Research Articles

Boundary Element Method for Viscous Flow through Out-phase Slip-patterned Microchannel under the Influence of Inclined Magnetic Field

Abstract In this paper, we investigate the impact of an inclined magnetic field of uniform intensity on viscous, incompressible pressure-driven Stokes flow through a slip-patterned, rectangular microchannel using the boundary element method based on the stream function-vorticity variables approach. The present investigation focuses only on the out-phase slip patterning of the microchannel walls. We address two scenarios of slip patterning, specifically large and fine slip patterning, which are determined by the periodicity of the patterning. We utilized the no-slip and Navier’s slip boundary conditions in an alternative manner on the walls. The Stokes equations govern the viscous flow through a microchannel. We assume a very small magnetic Reynold’s number to eliminate the equation of induced magnetic field in the present study. We analyzed the impact of considered dimensionless hydrodynamic parameters, including the Hartman number (Ha), inclination angle (θ) of the magnetic field, and the slip length (ls ) on fluid dynamics. In the case of fine slip, we observed significant variations in both velocity and pressure gradient, in contrast to large slip patterning. Fine slip patterning significantly increases the shear stress at slip regimes, while large slip periodicity significantly reduces it at no-slip regimes. The present investigation has several notable implications, such as regulation and advancement of mixing and heat transmission within microfluidic systems.

Evaluation of the interface shear strength of concrete containing treated and untreated reclaimed asphalt pavement aggregates

Concrete production significantly contributes to carbon emissions and depletes natural resources, promoting increased interest in substituting traditional materials with recycled alternatives. This shift necessitates a thorough examination of recycled materials’ performance in concrete applications. This study aims to provide a better understanding of the shear transfer behavior of concrete containing reclaimed asphalt pavement (RAP) aggregates. To mitigate the potential negative impact of RAP aggregates, mechanical treatment was applied. The effects of RAP inclusion and mechanical treatment on the shear transfer strength of reinforced concrete were investigated. A total of twenty-four push-off specimens were tested to investigate the effects of aggregate replacement and the clamping reinforcement ratio on specimen behavior. Digital image correlation (DIC) was employed to monitor the strains and displacements, allowing for detailed tracking of the interface slip, crack width development, and reinforcement strain. The findings revealed that replacing natural aggregates with RAP aggregates resulted in lower compressive strengths and lower ultimate push-off strengths of the resulting concrete at equivalent effective water-to-cement (w/c) ratios. Specimens with higher reinforcement ratios retained residual strengths up to 77 % of the ultimate strength, even at relatively large slip and crack widths. The clamping reinforcement ratio was identified as the key factor influencing the ultimate and residual strengths for both types of concrete. The experimental results were compared to the strengths calculated by the ACI, PCI, AASHTO, and CSA design equations to evaluate their applicability when different aggregate types are used.

The rupture process of the Hualien M7.3 sequence on April 3, 2024

The Hualien M7.3 earthquake on April 3, 2024, was a significant and strong earthquake in Taiwan, China in the past two decades. The rupture process of the main shock and strong aftershocks is of great significance to the subsequent seismic activity and seismogenic tectonic research. Based on local strong-motion data, we used the IDS (Iterative Deconvolution and Stacking) method to obtain the rupture process of the mainshock and two strong aftershocks on the 23rd. The rupture of the mainshock was mainly unilateral, lasting 31 s, with a maximum slip of 2 m, and the depth of the large slip zone is about 41-49 km. There is a clear difference between the rupture depth of the main shock and the two strong aftershocks. The depths of the large slip zones of the latter two are 3-9 km and 8-10 km, respectively. There is also a significant difference in the seismogenic fault between the mainshock and the aftershocks, and we believe that there are two seismogenic fault zones in the study area, the deep and the shallow fault zone. The slip of the deep faults activates the shallow faults.

Fault rock properties and conditions produce variance in slip during earthquake rupture propagation at the Nankai Trough

Although drilled samples of fault rocks have yielded information on frictional features of shallow subduction zones, the relationship of rupture propagation to the levels of friction and pore-fluid pressure remains uncertain. To investigate this topic, we performed dynamic rupture simulations along the megasplay fault that slipped during the 1944 Mw 8.0 Tonankai earthquake in the Nankai Trough. We used actual data from friction experiments on rocks from the fault segment and pre-existing pore pressures deduced from geophysical surveys for the shallow portion of 0–10 km depth along the fault. Simulations of low friction (friction coefficient ca. 0.04) produced large slip (about 30 m), whereas simulations using higher friction (friction coefficient ca. 0.2) suppressed the rupture. In simulations with low friction in which the pore-fluid pressure was nearly equal to the lithostatic stress, the slip decreased to about 25 m. However, when the simulations included slip-strengthening at shallow depth and higher friction, the slip still reached roughly 20 m. Such variability in slip during rupture propagation is caused by differences in the friction features and fluid pressure conditions of fault rocks, in which the friction features might be related to the mineral composition. Spatiotemporal heterogeneity in fault-rock type and their physical and hydraulic properties may fundamentally produce the complexity and variability of earthquake rupture propagation along the Nankai plate-subduction boundary.

Majority of Ruptures in Large Continental Strike-Slip Earthquakes Are Unilateral: Permissive Evidence for Hybrid Brittle-to-Dynamic Ruptures

Abstract Finite-element models of neotectonics require transform faults to rupture seismically even where preseismic shear stresses are low, presumably by dynamic-weakening mechanisms. A long-standing objection is that, if a rupture initiated at an asperity with high static friction stresses, which then transitioned to low dynamic-weakening stresses, local stress drop would be near total and on the order of 80 MPa, which is 4×–40× greater than observed. But the 5 Mw ≥ 7.8 transform earthquakes since 2000 initially ruptured on the branch faults of small net slip (Stein and Bird, 2024). If the slip initiates on a branch fault with different slip physics and no dynamic weakening, this solves the stress-drop problem. We propose that most large shallow earthquakes are hybrid ruptures, which begin on branch faults of small slip with high shear stresses, and then continue propagating on a connected dynamically weakened fault of large slip, even where shear stresses are low. One prediction of this model is that most large shallow ruptures should be unilateral. We test this prediction against the 100 largest (m ≥ 6.49) shallow continental strike-slip earthquakes 1977–2022, using information from the Global Centroid Moment Tensor and International Seismological Centre catalogs. The differences in time and location between the epicenter and the epicentroid define a horizontal “migration” velocity vector for the evolving centroid of each rupture. Early aftershock locations are summarized by a five-parameter elliptical model. Using the geometric relations between these (and mapped traces of active faults) and guided by a symmetrical decision table, we classified 55 ruptures as apparently unilateral, 30 as bilateral, and 15 as ambiguous. Our finding that a majority (55%–70%) of these ruptures are unilateral permits the interpretation that a majority of ruptures are hybrids, both in terms of geometry (branch fault to transform) and in terms of the physics of their fault slip.

Dynamic Rupture and Strong Ground-Motion Simulations of the 8 January 2022 Ms 6.9 Qinghai Menyuan Earthquake

Abstract The 2022 Ms 6.9 Qinghai Menyuan, China, earthquake is the most destructive earthquake to have occurred near the Lenglongling fault at the western segment of the Qilian–Haiyuan fault since 2016 Ms 6.4 Menyuan earthquake. The 2022 earthquake generated surface rupture measuring about 30 km with an unexpected maximum offset larger than 2.6 m in the epicentral area, and severely damaged the local infrastructure and transportation. To analyze the possible causes of the large surface slip and to reveal the rupture process, we modeled the dynamic rupture and strong ground motion of the 2022 Menyuan earthquake using the curved-grid finite-difference method. In the simulation, the geometry of the fault is constructed based on the observed trace of the surface ruptures. The background tectonic stress field is assumed to be uniform, and the slip-weakening law with the constant friction coefficients is adopted. Our modeling results showed that the rupture model with a focal depth of 6 km and a rupture width of 10 km provides a good fit to the observed surface slips and the field records. We also investigated the effects of the focal depth and the rupture size on the surface slips. It is found that under the same conditions, the dynamic rupture models with a larger rupture size generated greater coseismic slips at the surface. However, only the model with a relatively smaller rupture width produced an Mw∼6.7 event similar to the Menyuan earthquake. In contrast, the influence of the focal depth is less significant. The decrease of the focal depth only leads to a slight increase in surface slip. Our results illustrated that a surface-breaking rupture with a relatively narrow width may physically control the general characteristics of the earthquake. This study provides a new insight into the rupture dynamics of the 2022 Menyuan earthquake.

Localization of Deformation on Faults Driven by Fluids During the L’Aquila 2009 Earthquake

AbstractCoseismic rupture and aftershock development on a fault plane are complex and heterogeneous processes. The Mw 6.1 L’Aquila 2009 normal faulting earthquake is a perfect case to explore how fault geometry and rheology influence the rupture process and aftershocks distribution. In this study, we use for the first time a dense set of earthquake data to obtain enhanced images of the causative normal fault structure to the kilometer scale. The hypocenter of the emergent onset of the mainshock took place within a low Vp/Vs volume, while the large coseismic slip occurred a few kilometers above, as the rupture propagated through a high Vp and high Vp/Vs fluid‐filled rock volume. The increase of Vp/Vs in the fault hanging wall during the sequence suggests a strong dehydration in the earthquake asperity, with an upward fluid pressure migration along the fault toward the host rock volume. We propose that the localization of deformation on the fault plane is favored by high fluid pressure, while the spreading of aftershocks on a wide volume around the fault is driven by the depletion of fluids from the slipped portion of the fault plane and migration to small segments within the fault host rocks.

Shaking Table Test Study of Dynamic Response and Failure Mode of Liquefiable Layer Slope Under Earthquake

As a kind of special terrain, the landslide disaster generated by liquefiable layer slopes under earthquake has become a major engineering challenge due to its large scale and long slip distance. In order to study the seismic response and damage mode of liquefiable layer slope, this paper follows the research line of “geological generalization, physical modeling, and result analysis”, and takes the liquefiable layer slope in the upper reaches of the Yellow River Class secondary terrace as the object of study, generalizes the physical model of the slope, and carries out shaking table test. Based on the PGA amplification coefficient, Fourier analysis and HHT time–frequency characterization of the model slope, it was found that the PGA amplification coefficient increases gradually along the slope height, reaches the maximum value at the top of the liquefiable layer and then decreases gradually, which indicates that the liquefiable sand layer has an obvious energy dissipation effect, and on the horizontal direction of model slope is a tendency to the surface effect; the seismic waves at the discontinuous interface change drastically, and the Hilbert time–frequency spectrum transforms from multiple peaks to a single peak; with the increase of the intensity, the intrinsic frequency of the overall model decreases, and the high-frequency component within the liquefiable sand layer decreases from 5–15[Formula: see text]Hz to 0–5[Formula: see text]Hz, indicating that the liquefiable layer has a filtering effect; the damage process of the liquefiable layer slope is the tensile crack at the top of the slope — seismic subsidence at the top of the slope — the shear yielding at the angle of the slope — shear surface penetration at the slope face — overall slope instability and flow-slip damage. The research results will provide a reference for the study of the disaster mechanism of the liquefiable layer slope.

Roughness Evolution of Granite Flat Fracture Surfaces during Sliding Processes

Roughness is an essential factor affecting the shear process of discontinuous surfaces, and the evolution of roughness is closely related to the mechanical behavior of discontinuous surfaces. In this paper, with the help of granite specimens, a direct shear test was carried out on flat fracture surfaces obtained by sawing in order to study the evolution of roughness with shear slip. During the tests, the roughness evolution was evaluated using the arithmetic mean, root mean square and power spectral density of the roughness. The variation in these parameters all indicate that the friction surface with large slip tends to be rougher, at least under the loading conditions in this paper. And the increase in normal force will enhance this process, while the loading rate seems to have little effect on the roughness evolution. Finally, the analysis of the power spectral density shows that the roughness evolution in the spatial frequency of the profile line is mainly reflected in the middle– and low–frequency part, while the high–frequency part corresponding to the microscopic roughness body does not change much throughout the shear process.

Fracture mechanism of Ti60 alloy at high temperature in very high cycle fatigue regime and fatigue life prediction

AbstractThis study aimed to investigate the failure mechanism of Ti60 titanium alloy in the high cycle fatigue (HCF) and very high cycle fatigue (VHCF) at 500°C. Ti60 specimens were characterized before and after the 500°C VHCF test, and the fracture surfaces were observed. The results show that there are three different fatigue failure mechanisms at 500°C: (i) equiaxed primary α phase (αp) cleavage, forming small unsmooth facets and rough fracture area (RA), and fusion of non‐adjacent facets to grow into the main crack. (ii) αp grains gather into large grain with triple junction, and large grain slip to form large fusion facet. (iii) Oxide shedding and then cracks form. A dynamic recurrent neural network model is used to predict the fatigue life of Ti60 alloy, 91% of the overall predictions were within the scatter band of 3.0, and 100% were within the scatter band of 5.0.

Background and clustering characteristics of recent seismicity in Southwestern China

SUMMARY This study analysed seismicity in southwestern China (1 January 2008 to 30 June 2021) using the earthquake catalogue compiled by the China Earthquake Network Center and four different space–time Epidemic-Type Aftershock Sequence models: the 2-D point-source (PS) model, the 2-D finite-source (FS) model, the 3-D PS model and the 3-D FS model. Our objective was to understand the features of the background seismicity and the patterns of earthquake clusters to better evaluate the regional seismic hazard. We carefully investigated the aftershock sequences that followed 7 of the 10 MS ≥ 6.0 earthquakes that have struck this region since the occurrence of the 2008 Wenchuan MS 8.0 earthquake [i.e. the Panzhihua (31 August 2008; MS 6.0), Yaoan (9 July 2009; MS 6.0), Lushan (20 April 2013; MS 7.0), Ludian (3 August 2014; MS 6.5), Jinggu (7 October 2014; MS 6.6), Kangding (11 November 2014; MS 6.3) and Yangbi (21 May 2021; MS 6.4) earthquakes]. Our results revealed the following. (1) The background seismicity level for natural earthquakes is usually stable but can experience sudden change due to major events, such as the 2014 Ludian MS 6.5, and the 2014 Jinggu MS 6.6 events. Such changes in the background rate can reach 50 per cent. (2) Reservoir-induced earthquakes substantially increase the level of regional seismicity, indicating that they cannot be ignored when analysing natural seismicity and evaluating regional earthquake hazards. (3) Events triggered directly by the main shock occur mostly in regions adjacent to areas with large coseismic slip, showing a pattern complementary to the main shock ruptures.

Path Tracking Control Based on T-S Fuzzy Model for Autonomous Vehicles with Yaw Angle and Heading Angle

Existing vehicle-road models used for road tracking do not take into account the side slip angle, which leads to a reduction in road tracking accuracy in scenarios where the vehicle is at a large side slip angle, such as an emergency lane change. Consequently, this study presents a path-tracking control technique based on the T-S fuzzy model of heading angle vehicle autonomy. In this paper, based on the yaw angle-based vehicle tracking model, a heading angle-based tracking model considering the side slip angle is constructed. Second, since the vehicle speed varies with time, this paper selects the membership function of the vehicle speed to establish the T-S fuzzy model of autonomous vehicle based on the yaw angle and heading angle, respectively, and ensures the robustness and stability over the whole parameter space by the linear parameter variation robust H∞ controller. Then, cost functions based on the yaw angle and heading angle augmented error systems are created separately to optimize the system’s overall performance. Ultimately, simulation and experimentation confirm that the algorithm for control, which is based on the fuzzy model of the heading angle vehicle, has superior autonomous trajectory performance.

A finite geometry, inertia assisted coarsening-to-complexity transition in homogeneous frictional systems

The emergence of statistical complexity in frictional systems (where nonlinearity and dissipation are confined to an interface), manifested in broad distributions of various observables, is not yet understood. We study this problem in velocity-driven, homogeneous (no quenched disorder) unstable frictional systems of height H. The latter are described at the continuum scale within a realistic rate-and-state friction interfacial constitutive framework, where elasto-frictional instabilities emerge from rate-weakening friction. For large H, such frictional systems were recently shown to undergo continuous coarsening until settling into a spatially periodic traveling solution. We show that when the system’s height-to-length ratio becomes small — characteristic of various engineering and geophysical systems —, coarsening is less effective and the periodic solution is dynamically avoided. Instead, and consistently with previous reports, the system settles into a stochastic, statistically stationary state. The latter features slip bursts, whose slip rate is larger than the driving velocity, which are non-trivially distributed. The slip bursts are classified into two types: predominantly non-propagating, accompanied by small total slip and propagating, accompanied by large total slip. The statistical distributions emerge from dynamically self-generated heterogeneity, where both the non-equilibrium history of the interface and wave reflections from finite boundaries, mediated by material inertia, play central roles. Specifically, the dynamics and statistics of large bursts reveal a timescale ∼H/cs, where cs is the shear wave-speed. We discuss the robustness of our findings against variations of the frictional parameters, most notably affecting the magnitude of frictional rate-weakening, as well as against different interfacial state evolution laws. Finally, we demonstrate a reverse transition in which statistical complexity disappears in favor of the spatially periodic traveling solution. Overall, our results elucidate how relatively simple physical ingredients can give rise to the emergence of slip complexity.

New source model for the 1771 Meiwa tsunami along the southern Ryukyu Trench inferred from high-resolution tsunami calculation

The 1771 Meiwa tsunami which struck the southern Ryukyu Islands (Sakishima Islands) had greater than 22 m run-up height, leaving about 12,000 casualties in its wake. At many places, the tsunami inundation or lack of inundation is well recorded in historical documents. Several tsunami source models have been proposed for this event using historical records as constraints of tsunami calculations. Nevertheless, the source model remains under discussion. This study re-evaluated the tsunami wave source model of the 1771 Meiwa tsunami using high-resolution (10 m mesh) bathymetric and topographical data for tsunami calculation, the latest historical record dataset, and seismological knowledge. Results demonstrated that a tsunami earthquake along the southern Ryukyu Trench was the likely cause of the 1771 event. However, it is noteworthy that assumption of a large slip with 30 m is necessary for a shallow and narrow region (fault depth = 5 km, fault width = 30 km, Mw = 8.49) of the plate boundary in the Ryukyu Trench, which is far larger than previously thought. This requirement of very large initial water level change at the source might involve not only the fault rupture along the plate boundary but also deformation by splay faults, inelastic deformation of unconsolidated sediments near the trench axis, and/or giant submarine landslides. Results also show that the effects of fault parameters on the run-up were quite different depending on the offshore coral reef width. This phenomenon strongly constrained the fault width to 30 km. Our tsunami ray tracing analysis further revealed the effects of bathymetry on tsunami propagation. It is noteworthy that meter-long huge tsunami boulders tend to be distributed along the specific coasts at which the tsunami was concentrated by bathymetric effects. This finding suggests that past tsunamis, including the 1771 event, might have affected the specific coral reefs on Sakishima Islands repeatedly, which is crucially important for understanding the heterogeneous distribution of tsunami boulders. This feature might also be useful to elucidate the effects of large tsunamis on the corals and reefs because a direct comparison of coral reefs that are damaged and not damaged by tsunami waves is testable in narrow areas in the case of the Sakishima Islands.

Numerical modelling of earthquake sequences involving valving and pumping of fluids

SUMMARY Faults in the upper crust are sometimes thought to act as self-sealing valves, episodically releasing highly overpressured fluids trapped at greater depth during earthquakes. They are also often considered to be capable of actively pumping fluids into or out of faults in response to coseismic volumetric strain, thermal pressurization or other mechanisms. In this study, we investigate how these different types of behaviour (i.e. valve versus pump) are manifested in earthquakes. We do this using a 2-D plane strain model where frictional sliding on a thrust fault that is fed by a fluid source at its base is coupled to porous flow, thermal pressurization and strong variations in permeability. Our results show that thermal pressurization leads to dramatic dynamic weakening that produces earthquakes that propagate as slip pulses producing large stress drop, large slip and high slip velocities. On the other hand, valve-type behaviour typically produces smaller, less energetic earthquakes that commonly arrest before rupturing the entire fault. In some valve models, we observe complex compound ruptures and swarm activity, which is linked to the ascent of a propagating fluid pressure pulse driven by a large increase in permeability during sliding. Both pump and valve mechanisms can produce anomalously weak faults, though they are each associated with distinctly different fluid pressure and strength evolution over the seismic cycle and during rupture. Our models highlight the complex way in which fluids may interact with earthquakes, especially if valve and pump models coexist.

The Stability of Compressed Air Storage Underground Gas Storage Chamber and Steel Lining Structure Analysis

According to the address characteristics and structural characteristics of an underground artificial chamber gas storage, a structural model of an underground chamber including steel lining, flexible concrete, concrete lining, and surrounding rock is established, and two limit states of large slip and no contact between steel lining and flexible concrete are simulated respectively. The stability of the surrounding rock and lining structure are calculated and analyzed by using large finite element software FLAC3D and ANSYS respectively. The results show that the maximum deformation of the surrounding rock is 1.57cm, the maximum stress of steel lining is 262.08MPa, and the changes of plastic zone, displacement, and steel lining stress of surrounding rock can meet the design requirements of this project.

The evolution of stresses and shear resistance on rough faults at large slip

The deviation of natural faults from planarity results in geometric asperities and a heterogeneous stress field, affecting the shear resistance and slip behavior. This study numerically examines the effects of macro-scale roughness, local friction, and wear on the evolution of the fault stresses and shear resistance with fault slip. In contrast to previous analytical and numerical studies, our quasistatic simulations allow for slip at the scale of roughness wavelengths, large roughness levels, and fault opening. The simulations indicate that analytical solutions, which predict a linear increase in shear resistance with slip, significantly overestimate the additional shear resistance from roughness, i.e., the roughness drag, at slips larger than several percent of the minimum roughness wavelength. Starting with mated fault surfaces, the average shear resistance on the fault initially rises rapidly, but the rate of increase diminishes significantly as slip increases and the faults open. While initially showing a larger increase in shear resistance than self-similar faults, at large slip, self-affine faults exhibit significantly smaller shear resistance, and the background stresses are bounded without the need for off-fault inelastic deformation. Although smaller than predicted by analytical solutions, the simulation results show that the additional resistance from macro-scale geometrical irregularities is important, enabling rough faults to operate under shear resistance significantly larger than that from local friction only. The impact of macro-scale roughness increases with the roughness amplitude and the friction coefficient, while wear reduces fault stresses and average shear resistance. The combined effect of local friction and roughness is nonlinear, involving a direct contribution to shear resistance from roughness, as observed in prior studies, and an additional contribution from the increasing average normal stress with slip, which amplifies the resistance from friction and becomes more pronounced at large slip.

Liquefaction-induced large deformation method with automatic time-step mapping and interfacial interpolation improvement: Case study of the San Fernando dam

This study constructed the u-p equation for saturated soil in the meshfree global weak form within the arbitrary Lagrangian-Eulerian (ALE) framework to solve the problem of liquefaction-induced large deformation induced by strong seismic loading. The interpolation improvement in the near-interface zone was achieved using domain truncation optimization and the local nodal refinement algorithm, and the relative error of the shape function was proposed as the judgment criterion for radial basis function (RBF) field variable mapping. Finally, a high precision meshfree liquefaction-induced large deformation method (MFLLDM) with automatic time step mapping and interfacial zone interpolation improvement was established. The MFLLDM achieved the efficient dynamic updating of soil stress state and pore pressure information in space during the deformation process and more accurately depicted the localized large shear deformation characteristics of the soil in the near-interface zone. The proposed method was integrated into a custom-built large-scale finite element method (FEM) computing system (GEODYNA), based on object-oriented and super element program design technology. This enabled the efficient coupled analysis of MFLLDM in the liquefied-induced large deformation zone and FEM in the small deformation zone. Numerical consolidation examples were used to validate the accuracy, convergence, and applicability of the soil elasto-plastic model for MFLLDM. The method was then used to simulate the liquefaction damage of the San Fernando dam, which reproduced the development process of large slip deformation and high-water pressure ratio for the upstream dam slope.

Unraveling the roles of fault asperities over earthquake cycles

Fault asperities can produce concentrated slips during large earthquakes and intensify damage. However, how asperities control fault behavior during other phases of the earthquake cycle remains poorly known. Here, we conduct friction experiment on a laboratory fault featuring two prominent geometric asperities to directly image their influences on long-term and short-term seismicity, and the nucleation and propagation of the mainshock. The laboratory observations and supporting numerical simulations reveal that one asperity located in the fault interior behaves firstly as a mechanical attractor to long-term seismicity, then a barrier to preseismic slow slip, and ultimately a source of large coseismic slip. In contrast, another asperity located near the fault margin primarily undergoes persistent aseismic slip throughout most phases of the earthquake cycle. These results provide new insights into how asperities partition strain across a broad spatiotemporal domain, establishing a physical link between long-term and short-term fault behaviors and the occurrence of large earthquakes.

Detection of a recent large Hyuga-nada long-term slow slip event and estimation of its spatiotemporal slip distributions

Using horizontal and vertical GNSS time series data from the GSI in sotheastern Kyushu from January 1, 2017 to June 30, 2022, we detected a recent long-term slow slip event (L-SSE) that occurred in the Hyuga-nada region, southwest Japan, and estimated its spatiotemporal slip distribution. We performed such analysis considering the piecewise evolution of slip over time with time windows of 0.8 years to ensure a good signal-to-noise ratio in the horizontal displacements in each time window. The results showed that a slip of more than 2 cm occurred in the central part of Miyazaki Prefecture from 2018.5 to 2019.3 (the units were 0.1 years = 36.5 days). Then, the amount of slip increased and expanded slightly to the southern part of Miyazaki Prefecture from 2019.3 to 2020.1, and the amount of slip reached a maximum of 3.9 cm in the subsequent period (2020.1–2020.9). A smaller slip occurred at almost the same location in the following 0.8 years. Therefore, the duration of the L-SSE was approximately 3.2 years from approximately 2018.5 to 2021.7. The annual average maximum slip rate was approximately 4.9 cm/yr (3.9 cm/0.8 yr) during the period from 2020.1 to 2020.9. The maximum total slip was estimated to be approximately 12.9 cm, and the equivalent release moment was 4.9 × 1019 Nm, corresponding to Mw7.1. Compared to previous L-SSEs that occurred in the Hyuga-nada region, the annual average maximum slip rate was relatively low, and the main total slip was estimated at almost the same location on the plate interface. On the other hand, the slip duration of this L-SSE was the longest, and the release moment and moment magnitude were the greatest. The total slip area of more than approximately 10 cm in the Hyuga-nada L-SSE estimated in this study almost overlapped with the afterslip area of the December 3, 1996 Hyuga-nada earthquake, with a depth range of approximately 30–40 km. Short-term SSEs also occurred in the same depth range. These coincidental same depth ranges of interplate seismic events stem from large thermal gradients. The large thermal gradients play an important role in narrowing the depth range of frictional parameters (a–b) at the plate interface. In addition, the pore pressure and normal stress are closely related to the critical stiffness. In particular, contrasting values of pore pressure and normal stress at the shallow and deep sides of the mantle wedge corner may also contribute to transient aseismic slip within the same depth range.