The slow self-arresting nature of low-frequency earthquakes

Low-frequency earthquakes are a series of recurring small earthquakes that are thought to compose tectonic tremors. Compared with regular earthquakes of the same magnitude, low-frequency earthquakes have longer source durations and smaller stress drops and slip rates. The mechanism that drives their unusual type of stress accumulation and release processes is unknown. Here, we use phase diagrams of rupture dynamics to explore the connection between low-frequency earthquakes and regular earthquakes. By comparing the source parameters of low-frequency earthquakes from 2001 to 2016 in Parkfield, on the San Andreas Fault, with those from numerical simulations, we conclude that low-frequency earthquakes are earthquakes that self-arrest within the rupture patch without any introduced interference. We also explain the scaling property of low-frequency earthquakes. Our findings suggest a framework for fault deformation in which nucleation asperities can release stress through slow self-arrest processes.

Anatomy of strike-slip fault tsunami genesis

Tsunami generation from earthquake-induced seafloor deformations has long been recognized as a major hazard to coastal areas. Strike-slip faulting has generally been considered insufficient for triggering large tsunamis, except through the generation of submarine landslides. Herein, we demonstrate that ground motions due to strike-slip earthquakes can contribute to the generation of large tsunamis (>1 m), under rather generic conditions. To this end, we developed a computational framework that integrates models for earthquake rupture dynamics with models of tsunami generation and propagation. The three-dimensional time-dependent vertical and horizontal ground motions from spontaneous dynamic rupture models are used to drive boundary motions in the tsunami model. Our results suggest that supershear ruptures propagating along strike-slip faults, traversing narrow and shallow bays, are prime candidates for tsunami generation. We show that dynamic focusing and the large horizontal displacements, characteristic of strike-slip earthquakes on long faults, are critical drivers for the tsunami hazard. These findings point to intrinsic mechanisms for sizable tsunami generation by strike-slip faulting, which do not require complex seismic sources, landslides, or complicated bathymetry. Furthermore, our model identifies three distinct phases in the tsunamic motion, an instantaneous dynamic phase, a lagging coseismic phase, and a postseismic phase, each of which may affect coastal areas differently. We conclude that near-source tsunami hazards and risk from strike-slip faulting need to be re-evaluated.

A Spectral Boundary-Integral Method for Quasi-Dynamic Ruptures of Multiple Parallel Faults

ABSTRACT Numerical models of rupture dynamics provide great insights into the physics of fault failure. However, resolving stress interactions among multiple faults remains challenging numerically. Here, we derive the elastostatic Green’s functions for stress and displacement caused by arbitrary slip distributions along multiple parallel faults. The equations are derived in the Fourier domain, providing an efficient means to calculate stress interactions with the fast Fourier transform. We demonstrate the relevance of the method for a wide range of applications, by simulating the rupture dynamics of single and multiple parallel faults controlled by a rate- and state-dependent frictional contact, using the spectral boundary integral method and the radiation-damping approximation. Within the antiplane strain approximation, we show seismic cycle simulations with a power-law distribution of rupture sizes and, in a different parameter regime, sequences of seismogenic slow-slip events. Using the in-plane strain approximation, we simulate the rupture dynamics of a restraining stepover. Finally, we describe cycles of large earthquakes along several parallel strike-slip faults in three dimensions. The approach is useful to explore the dynamics of interacting or isolated faults with many degrees of freedom.

Preseismic atmospheric radon anomaly associated with 2018 Northern Osaka earthquake

Despite the challenges in identifying earthquake precursors in intraplate (inland) earthquakes, various hydrological and geochemical measurements have been conducted to establish a possible link to seismic activities. Anomalous increases in radon (222Rn) concentration in soil, groundwater, and atmosphere have been reported prior to large earthquakes. Although the radon concentration in the atmosphere is lower than that in groundwater and soils, a recent statistical analysis has suggested that the average atmospheric concentration over a relatively wide area reflects crustal deformation. However, no study has sought to determine the underlying physico-chemical relationships between crustal deformation and anomalous atmospheric radon concentrations. Here, we show a significant decrease in the atmospheric radon concentration temporally linked to the seismic quiescence before the 2018 Northern Osaka earthquake occurring at a hidden fault with complex rupture dynamics. During seismic quiescence, deep-seated sedimentary layers in Osaka Basin, which might be the main sources of radon, become less damaged and fractured. The reduction in damage leads to a decrease in radon exhalation to the atmosphere near the fault, causing the preseismic radon decrease in the atmosphere. Herein, we highlight the necessity of continuous monitoring of the atmospheric radon concentration, combined with statistical anomaly detection method, to evaluate future seismic risks.

Seismogenic potential of the Main Himalayan Thrust constrained by coupling segmentation and earthquake scaling

<p>Recent studies have shown that the Himalayan region is under the threat of earthquakes of magnitude 9 or larger. These estimates are based on comparisons of the geodetically inferred moment deficit rate with the seismicity of the region. However, these studies do not account for the physics of fault slip, specifically the influence of frictional barriers on earthquake rupture dynamics, which controls the extent and therefore the magnitude of large earthquakes. Here, we propose a methodology for incorporating outcomes of physics-based earthquake cycle models into hazard estimates. The methodology takes also into account the moment deficit rate, the magnitude-frequency of the current and historical catalogs, and the moment-area earthquake scaling law.</p><p>For the Himalaya setting, we estimate an improved probabilistic estimate moment deficit rate using coupling estimates inferred using a Bayesian framework. The locking distribution of the fault suggests an along-strike segmentation of the MHT with three segments that may act as aseismic barriers. The effect of the barriers on rupture propagation is assessed using results from dynamic models of the earthquake cycle. We show that, accounting for measurement and methodological uncertainties, the MHT is prone to rupturing in M8.7 earthquakes every T>200 yr, with M>9.5 events being greatly improbable. The methodology also allows to estimate the probability of the position of earthquakes on the fault based on the effect of the seismic barriers and their magnitude. This study provides a straightforward and computationally efficient method for estimating regional seismic hazard accounting for the full physics of fault slip.</p>

Online measurement of floc size, viscosity, and consistency of cellulose microfibril suspensions with optical coherence tomography

In this study, cellulose microfibril (CMF) suspensions were imaged during pipe flow at consistencies of 0.4%, 1.0%, and 1.6% with optical coherence tomography (OCT) to obtain images of the structure and the local velocity of the suspension. The viscosities obtained by combining pressure loss measurement with the OCT velocity data showed typical shear thinning behavior and were in excellent agreement with viscosities obtained with ultrasound velocity profiling. The structural OCT images were used to calculate the radial and the axial floc sizes of the suspension. A fit of power law to the geometrical floc size–shear stress data gave the same power law index for all consistencies, suggesting that floc rupture dynamics is independent of consistency. The dependence of viscosity and floc size on shear stress was similar, indicating that the shear thinning behavior of CMF suspensions is closely related to the rupture dynamics of flocs. The results also showed that an apparent attenuation coefficient of the OCT signal can be used to determine the consistency of CMF suspensions.