No Arabic abstract
Earthquakes at seismogenic plate boundaries are a response to the differential motions of tectonic blocks embedded within a geometrically complex network of branching and coalescing faults. Elastic strain is accumulated at a slow strain rate of the order of $10^{-15}$ s$^{-1}$, and released intermittently at intervals $>100$ years, in the form of rapid (seconds to minutes) coseismic ruptures. The development of macroscopic models of quasi-static planar tectonic dynamics at these plate boundaries has remained challenging due to uncertainty with regard to the spatial and kinematic complexity of fault system behaviors. In particular, the characteristic length scale of kinematically distinct tectonic structures is poorly constrained. Here we analyze fluctuations in GPS recordings of interseismic velocities from the southern California plate boundary, identifying heavy-tailed scaling behavior. This suggests that the plate boundary can be understood as a densely packed granular medium near the jamming transition, with a characteristic length scale of $91 pm 20$ km. In this picture fault and block systems may rapidly rearrange the distribution of forces within them, driving a mixture of transient and intermittent fault slip behaviors over tectonic time scales.
Seismogenic plate boundaries are presumed to behave in a similar manner to a densely packed granular medium, where fault and blocks systems rapidly rearrange the distribution of forces within themselves, as particles do in slowly sheared granular systems. We use machine learning and show that statistical features of velocity signals from individual particles in a simulated sheared granular fault contain information regarding the instantaneous global state of intermittent frictional stick-slip dynamics. We demonstrate that combining features built from the signals of more particles can improve the accuracy of the global model, and discuss the physical basis behind decrease in error. We show that the statistical features such as median and higher moments of the signals that represent the particle displacement in the direction of shearing are among the best predictive features. Our work provides novel insights into the applications of machine learning in studying frictional processes that take place in geophysical systems.
The particle discrete element simulation of the instability and failure process of the granular slope accumulator model when the metal plate continues downward is obtained, and the two-dimensional total velocity vector of soil particle velocity and slope slip during the instability and failure of the slope accumulator are obtained. Macro-response processes such as removing the angle of the crack surface and the average velocity in the y-direction of the slope top of the slope accumulation body. Construct a normal force chain undirected network model of the slope accumulation body particles under natural accumulation, and study the location of its slip surface, and The results are compared with the experimental results. Finally, the complex network method is used to analyze the topological characteristics of the contact force chain network of the particles on the slope top of the slope accumulation body, and the average degree, clustering coefficient and average shortest path are obtained during the slope instability of the slope accumulation body. The evolutionary rule of the method is used to verify its accuracy in combination with the strength reduction method. The research results show that the average shortest path can provide a more effective early warning of the instability and failure of slope deposits. A complex network theory is used to study the macro response of the slope deposits and its force chain. The interrelationship between the macroscopic structure of the network provides a new mathematical analysis method for the study of slope instability.
We report results of 3D Discrete Element Method (DEM) simulations aiming at investigating the role of the boundary vibration in inducing frictional weakening in sheared granular layers. We study the role of different vibration amplitudes applied at various shear stress levels, for a granular layer in the stick-slip regime and in the steady-sliding regime. Results are reported in terms of friction drops and kinetic energy release associated with frictional weakening events. We find that larger vibration amplitude induces larger frictional weakening events. The results show evidence of a threshold below which no induced frictional weakening takes place. Friction drop size is found to be dependent on the shear stress at the time of vibration. A significant increase in the ratio between the number of slipping contacts to the number of sticking contacts in the granular layer is observed for large vibration amplitudes. These vibration-induced contact rearrangements enhance particle mobilization and induces a friction drop and kinetic energy release. This observation provides some insight into the grain-scale mechanisms of frictional weakening by boundary vibration in a dense sheared granular layer. In addition to characterizing the basic physics of vibration induced shear weakening, we are attempting to understand how a fault fails in the earth under seismic wave forcing. This is the well know phenomenon of dynamic earthquake triggering. We believe that the granular physics are key to this understanding.
The coupled mechanics of fluid-filled granular media controls the behavior of many natural systems such as saturated soils, fault gouge, and landslides. The grain motion and the fluid pressure influence each other: It is well established that when the fluid pressure rises, the shear resistance of fluid-filled granular systems decreases, and as a result catastrophic events such as soil liquefaction, earthquakes, and accelerating landslides may be triggered. Alternatively, when the pore pressure drops, the shear resistance of these systems increases. Despite the great importance of the coupled mechanics of grains-fluid systems, the basic physics that controls this coupling is far from understood. We developed a new multi-scaled model based on the discrete element method, coupled with a continuum model of fluid pressure, to explore this dynamical system. The model was shown recently to capture essential feedbacks between porosity changes arising from rearrangement of grains, and local pressure variations due to changing pore configurations. We report here new results from numerical experiments of a continuously shearing layer of circular two-dimensional grains, trapped between two parallel rough boundaries. The experiments use a fixed confining stress on the boundary walls, and a constant velocity applied to one of the boundaries, as if this system was the interior of a sliding geological fault filled with fault gouge. In addition, we control the layer permeability and the drainage boundary conditions. This paper presents modeling results showing that the localization of shear (into a narrow shear band within the shearing layer) is strongly affected by the presence of fluids. While in dry granular layers there is no preferred position for the onset of localization, drained systems tend to localize shear on their boundary. We propose a scaling argument to describe the pressure deviations in a shear band, and use that to predict the allowable positions of shear localizations as a function of the fault and gouge properties.
We explore interactions of elastic waves propagating in plates (with soil parameters) structured with concrete pillars buried in the soil. Pillars are 2 m in diameter, 30 m in depth and the plate is 50 m in thickness. We study the frequency range 5 to 10 Hz, for which Rayleigh wave wavelengths are smaller than the plate thickness. This frequency range is compatible with frequency ranges of particular interest in earthquake engineering. It is demonstrated in this paper that two seismic cloaks configurations allow for an unprecedented flow of elastodynamic energy associated with Rayleigh surface waves. The first cloak design is inspired by some approximation of ideal cloaks parameters within the framework of thin plate theory. The second, more accomplished but more involved, cloak design is deduced from a geometric transform in the full Navier equations that preserves the symmetry of the elasticity tensor but leads to Willis equations, well approximated by a homogenization procedure, as corroborated by numerical simulations. The two cloakss designs are strikingly different, and the superior efficiency of the second type of cloak emphasizes the necessity for rigor in transposition of existing cloakss designs in thin plates to the geophysics setting. Importantly, we focus our attention on geometric transforms applied to thick plates, which is an intermediate case between thin plates and semi-infinite media, not studied previously. Cloaking efficiency (reduction of the disturbance of the wave wavefront and its amplitude behind an obstacle) and protection (reduction of the wave amplitude within the center of the cloak) are studied for ideal and approximated cloaks parameters. These results represent a preliminary step towards designs of seismic cloaks for surface Rayleigh waves propagating in sedimentary soils structured with concrete pillars.