- About
- Editorial board
- Articles
- Special issues
- Highlight articles
- Manuscript tracking
- Book reviews
- Subscribe to alerts
- Peer review
- For authors
- For reviewers
- EGU publications
Research article 27 May 2020
Research article | 27 May 2020
Correspondence: Patricio Venegas-Aravena (plvenegas@uc.cl)
HideCorrespondence: Patricio Venegas-Aravena (plvenegas@uc.cl)
Recently, it has been shown theoretically how the lithospheric stress changes could be linked with magnetic anomalies, frequencies, spatial distribution and the magnetic-moment magnitude relation using the electrification of microfractures in the semibrittle–plastic rock regime (Venegas-Aravena et al., 2019). However, this seismo-electromagnetic theory has not been connected with the fault's properties in order to be linked with the onset of the seismic rupture process itself. In this work we provide a simple theoretical approach to two of the key parameters for seismic ruptures which are the friction coefficient and the stress drop. We use sigmoidal functions to model the stress changes in the nonelastic regime within the lithosphere. We determine the temporal changes in frictional properties of faults. We also use a long-term friction coefficient approximation that depends on the fault dip angle and four additional parameters that weigh the first and second stress derivative, the spatial distribution of the nonconstant stress changes, and the stress drop. We found that the friction coefficient is not constant in time and evolves prior to and after the earthquake occurrence regardless of the (nonzero) weight used. When we use a dip angle close to 30^{∘} and the contribution of the second derivative is more significant than that of the first derivative, the friction coefficient increases prior to the earthquake. During the earthquake event the friction drops. Finally, the friction coefficient increases and decreases again after the earthquake occurrence. It is important to mention that, when there is no contribution of stress changes in the semibrittle–plastic regime, no changes are expected in the friction coefficient.
The electromagnetic phenomena that could be linked with earthquake occurrences are usually considered within the lithosphere–atmosphere–ionosphere-coupling effect (or LAIC effect; e.g., De Santis et al., 2019a). Some of these electromagnetic phenomena have been recorded both prior to and after earthquakes using different methodologies, data, earthquakes and instrumentation. For example, some researchers have shown coseismic magnetic variations during some earthquakes (e.g., Utada et al., 2011, during the 2011 Tōhoku earthquake), while other researchers have focused on the oscillation frequency (microhertz–kilohertz range) of the magnetic field prior to the occurrence of some earthquakes (Schekotov and Hayakawa, 2015; Cordaro et al., 2018; Potirakis et al., 2018a, b, among others). One of the most recent methodologies corresponds to the measurement of magnetic anomalies, which detects the number of magnetic peaks (or events) that exceed a given threshold that represents the normal conditions of the magnetic field over long periods of time (De Santis et al., 2017). There is some research that shows that the increases in the number of anomalies could be linked with impending earthquakes. For instance, De Santis et al. (2019b) have recently found an increase in the number of daily magnetic anomalies prior to 12 major earthquakes between 2014 and 2016. These magnetic anomalies were confirmed by other researchers (e.g., Marchetti and Akhoondzadeh, 2018).
It is well known from laboratory experiments that rock samples undergoing fast changes in the semibrittle–plastic rock regime generate microfractures, displacement of dislocations and electrification (e.g. Anastasiadis et al., 2004; Zhang et al., 2020 and references therein). This physical mechanism of rock electrification is described mathematically by the motion of charged edge dislocations (MCD) model (e.g., see Vallianatos and Tzanis, 1998, 2003, for a comprehensive derivation of the MCD model). This model is a plausible electromechanical mechanism that could explain the magnetic measurements. Starting from the experimental evidence Venegas-Aravena et al. (2019) developed a seismo-electromagnetic theory based on microcracks and stress changes. This theory showed how the fractal nature of the cracks explains the observed magnetic frequency range, the coseismic magnetic field and the conditions for generating magnetic anomalies. However, this theory (like the others, e.g., Freund, 2003; De Santis et al., 2019a) does not explain any change in the parameters that control the generation of seismic ruptures using magnetic measurements. This prevents the linking of seismo-electromagnetism with classical seismology. In this work, we address this link using one of the key points that controls seismic rupture and slip on the fault: friction force and stress drop. Our model uses the tectonic geometry and stress drop of the M_{w} 8.8 2010 Maule (Chile) earthquake in order to provide real data and validate our analysis. The paper is divided as follows: in Sect. 2 we consider the friction coefficient, adding the brittle–plastic stress changes contribution to the usual elastic stress. In Sect. 3 we discuss the temporal variations in the brittle–plastic friction and their connections with the spatial distribution on the fault and in the lithosphere. In Sect. 4 we provide the relation between the stress drop and the coseismic magnetic field. The spatial–temporal friction coefficient along the fault is calculated by adding the elastic stress drop. In Sect. 5, the Gutenberg–Richter law is written in terms of the semibrittle–plastic shear stress. The rupture time is discussed in Sect. 6. Finally, the discussion takes place and conclusions are drawn in Sect. 7.
The standard friction force can be understood as the complex dissipation of mechanical energy in the form of plastic or elastic deformation of asperities (mechanical interaction), thermal dissipation (heat) and the adhesion (interatomic interaction) of two sliding surfaces (e.g., Sun and Mosleh, 1994). When we consider a particular contact area between two dry surfaces, the static friction coefficient μ that describes this interaction can be written approximately as the ratio of the shear τ and the normal N stress (load) as shown by Eq. (1) (e.g., Byerlee, 1978; Chen, 2014, and references therein).
The static friction coefficient μ gives some information about the contact behavior. For instance, μ tends to be large when the contact area is increased due to the surface's plastic deformation (e.g., Chen, 2014). If we also consider the pure plastic regime, we can add a small plastic shear stress (δτ_{plastic}) and a small plastic normal stress (δN_{plastic}) contribution into Eq. (1) leading to the following expressions:
If we do not consider the pure plastic effects, the plastic contribution vanishes and the ratio of τ_{elastic} and N_{elastic} describes the usual (nonlinear) friction behaviors that occur during the complete frictional cycle: presliding (increase in friction coefficient when there is no apparent or residual displacement), gross sliding (observable displacement and decrease in friction coefficient) and healing (friction coefficient recovery; e.g., Parlitz et al., 2004; Marone and Saffer, 2015; Papangelo et al., 2015; and references therein).
Venegas-Aravena et al. (2019) state that the electrification within rocks is mainly due to a nonconstant stress change during the semibrittle–plastic transition. This means that the temporal changes in the semibrittle–plastic stress (δσ_{sbp}) rules the total plastic stress (δσ_{plastic}). Thus, it implies that the plastic shear and normal stress can be written in terms of a linear combination of the temporal changes in δσ_{sbp} as shown in Eqs. (4) and (5).
where k_{1}, k_{2}, k_{3} and k_{4} are dimensionless constants to be determined later, δt and (δt)^{2} are the temporal variations from the first- and second-order time contribution, respectively. The temporal variations in the shear and normal semibrittle–plastic stress contributions are ${\dot{\mathit{\tau}}}_{\mathrm{sbp}}$, ${\ddot{\mathit{\tau}}}_{\mathrm{sbp}}$, ${\dot{N}}_{\mathrm{sbp}}$ and ${\ddot{N}}_{\mathrm{sbp}}$, respectively. The series expansions in Eqs. (4) and (5) are convenient and allow us to write the plastic contribution as a sum of stresses that depend on time. This is relevant because the seismo-electromagnetic theory seeks to relate the temporal variable of the earthquake and the stress within the lithosphere. From now on, we refer to the semibrittle–plastic stress as the uniaxial stress σ.
It is possible to relate this shear, and the normal stresses from Eqs. (4)–(5), with the uniaxial stress using the geometry shown in Fig. 1. In this figure, we have a simple schematic representation of the lithosphere under uniaxial stress change dσ∕dt in the presence of a fault with a static friction coefficient μ and a dip angle of 30^{∘}. We can write this uniaxial stress change in terms of the dip angle θ in the normal and tangential direction (in red on Fig. 1) of the fault, as shown in the following expression:
where dσ∕dt corresponds to the magnitude of uniaxial temporal stress change. Using this expression, we write the brittle–plastic stress contributions in terms of uniaxial stress change as
Substituting Eqs. (2)–(8) into Eq. (1) and considering a second-order linear combination, we obtain the static friction coefficient as a function of time, the fault angle and the semibrittle–plastic changes within the lithosphere given by
where the dots above σ are the notation for the first and second temporal derivative of the uniaxial stress. According to Venegas-Aravena et al. (2019), the temporal stress change $\dot{\mathit{\sigma}}$ has approximately a sigmoidal shape, which can be written as $\dot{\mathit{\sigma}}\left(t\right)=a/(b+{e}^{({t}_{\mathrm{0}}-t)\times w})$, where a, b, w and t_{0} are constants. In Fig. 2 we represent $\dot{\mathit{\sigma}}$ and $\ddot{\mathit{\sigma}}$ as a function of time where we have used $a=b=w=\mathrm{1}$ and t_{0}=10 a.u. (arbitrary units). According to De Santis et al. (2019b), most of the earthquakes recorded occurred close to the center of Fig. 2, which is when t_{EQ}=t_{0}. For instance, Marchetti and Akhoondzadeh (2018) have shown that the M_{w} 8.2 2017 Chiapas (Mexico) earthquake occurred after this time (t_{EQ}>t_{0}). They also use daily values of magnetic anomalies ($B\propto \mathrm{d}\mathit{\sigma}/\mathrm{d}t$, which comes directly from the experimental equation $I={\mathit{\alpha}}_{\mathrm{0}}\mathrm{d}\mathit{\sigma}/\mathrm{d}t$, where I is the electric current and α_{0} is a constant of proportionality; see Vallianatos and Triantis, 2008, and references therein), thus δt=1 d =86 400 s.
Let us now estimate the values of the constants of Eq. (9). First, we consider the dip or subduction angle θ. According to Maksymowicz (2015), this angle is close to 20^{∘} at the depth (∼30 km) and location (35^{∘}54^{′}32^{′′} S, 72^{∘}43^{′}59^{′′} W) of the 2010 Maule earthquake. Maksymowicz (2015) claimed that the static friction coefficient in the Chilean convergent margin is close to μ_{ch}≈0.5. Lamb (2006) calculated that the initial value of τ_{0} is 15.4 MPa (constant) in southern Chile. Using μ_{ch}, it is expected that N_{0} should be approximately equal to 30.8 MPa (constant).
In addition, rock experiments show that the values of $\dot{\mathit{\sigma}}$ are close to 1 MPa s^{−1} (e.g., Saltas et al., 2018). This implies that $\left|{k}_{\mathrm{1}}\right|$ and $\left|{k}_{\mathrm{3}}\right|$ must be close to $\sim {\mathrm{10}}^{-\mathrm{4}}$, in order to balance the δt factor. The values of k_{2} and k_{4} must be equal to or less than $\sim {\mathrm{10}}^{-\mathrm{9}}$, otherwise the values of the friction coefficient would be greater than 1. If we consider an initial increase in the normal stress, the sign of the constants should be negative for k_{3}, k_{4} and positive for k_{1}, k_{2}. In Fig. 3 we can see how the friction coefficient changes in time, when using values of k_{1} and k_{3} described above and different values of k_{2} and k_{4} (second-order contribution). When we use values on k_{2} and k_{4} of the order of $\sim {\mathrm{10}}^{-\mathrm{9}}$, it is possible to observe how the friction decreases after the earthquake. Furthermore, it is important to note that the earthquake does not occur when the friction has its maximum value but when it is close to its maximum. When we use values of k_{2} and k_{4} that are similar or less than $\sim {\mathrm{10}}^{-\mathrm{10}}$, the contribution of $\ddot{\mathit{\sigma}}$ in Eq. (9) vanishes.
Another critical point is related to the differential time. For instance, when we consider δt≤1 s, the semibrittle–plastic stress term vanishes and the usual friction is recovered. This fact is especially relevant because the semibrittle–plastic contribution to the friction coefficient seems to be relevant only on large timescales. To summarize, the friction coefficient of Eq. (9) could be seen as a generalization of the standard friction coefficient when large timescales are considered.
In the previous section we linked the friction coefficient to the semibrittle–plastic regime that generates microcracks and electrification within the rocks. An important caveat is that this phenomenon does not occur everywhere. For instance, Dobrovolsky et al. (1979) described a specific “preparation zone” required close to the future hypocenter in order to accumulate sufficient stress to trigger the earthquake. This criterion has been widely used by researchers to establish a threshold for determining where the magnetic measurements can be associated with earthquakes (e.g., De Santis et al., 2019a, b). In other words, the phenomenon is local. However, if this is applied, a variation in the friction coefficient related to the rock electrification phenomenon close to the fault would be expected, while a friction variation outside the zones of semibrittle–plastic influence would not.
In order to consider the local variations, we add to our model a separable spatial function of the uniaxial stress as given by $\stackrel{\mathrm{\u203e}}{\mathit{\sigma}}(x,t)=\mathit{\gamma}\left(x\right)\mathit{\sigma}\left(t\right)$ (for simplicity we choose a uniaxial dependence), where σ(t) corresponds to the same time function as considered in Sect. 2 and γ(x) is the dimensionless spatial distribution parallel to the fault (see the coordinate system in Fig. 4). Let us mention that the values of γ(x) are different when constant and nonconstant stress changes are considered. With this in mind, we can rewrite Eq. (9) as
After straightforward derivations, we get that the gradient of the friction coefficient along the fault using Eq. (10) is given by
where $A=({k}_{\mathrm{1}}\dot{\mathit{\sigma}}+{k}_{\mathrm{2}}\ddot{\mathit{\sigma}}\mathit{\delta}t)\mathrm{cos}\mathit{\theta}\mathit{\delta}t$, $B=({k}_{\mathrm{3}}\dot{\mathit{\sigma}}+{k}_{\mathrm{4}}\ddot{\mathit{\sigma}}\mathit{\delta}t)\mathrm{sin}\mathit{\theta}\mathit{\delta}t$, $\mathit{\alpha}={N}_{\mathrm{0}}{k}_{\mathrm{1}}\mathrm{cos}\mathit{\theta}+{\mathit{\tau}}_{\mathrm{0}}{k}_{\mathrm{3}}\mathrm{sin}\mathit{\theta}$ and $\mathit{\beta}={N}_{\mathrm{0}}{k}_{\mathrm{2}}\mathrm{cos}\mathit{\theta}+{\mathit{\tau}}_{\mathrm{0}}{k}_{\mathrm{4}}\mathrm{sin}\mathit{\theta}$. Equations (10) and (11) imply that, if we distribute the constant and nonconstant temporal stress change in a nonuniform manner along the fault, then it could be said that this phenomenon is local. This can be understood using the example of Fig. 4: Fig. 4a shows a blue area of length L where we have a nonconstant temporal stress change. There one observes a change in the friction coefficient inside the projected gray area (Fig. 4b). Conversely, outside the area of length L (Fig. 4c), we assume a constant stress change, and there is no change in the friction coefficient (Fig. 4d). Let us remind ourselves that the entire fault suffered from stress accumulation during the entire process. However, only the gray area is directly affected by the friction change. This shows that the temporal friction changes are restricted only to a specific area (gray area) on the fault. Hence, one observes a nonzero friction coefficient gradient on the fault and we have a local phenomenon ($\mathrm{\nabla}\stackrel{\mathrm{\u203e}}{\mathit{\mu}}\ne \mathrm{0}\iff \mathrm{Local}$).
The process described in Fig. 4 reveals why the magnetic measurements are not global. It also validates the locality criteria in terms of fault properties. Furthermore, the spatial distribution of magnetic field measurements is expected to have a length scale comparable to the friction coefficient change due to the dependency of γ(x) in Eq. (10). That is, the larger the detection area of magnetic anomalies, the larger the area where fault friction is changing. In addition, Venegas-Aravena et al. (2019) described how the uniaxial stress change implies a change in the b-value of the Gutenberg–Richter law. It also implies that a significant earthquake is needed in order to satisfy this change in the Gutenberg–Richter law. Therefore, a greater expected magnitude could be related to changes in the coefficient of friction within an area within the fault. However, changes in the friction coefficient do not directly imply the earthquakes generation itself.
Up to this point, no changes in the elastic stresses have been considered in Eq. (10). This section focuses on one of the elastic parameter that is involved in the seismic rupture process: the stress drop Δτ. This parameter is one of the most relevant because it shows the shear stress differences prior to and after the onset of the earthquake within the fault rupture area (e.g., Aki, 1966). Furthermore, this parameter is also related to the seismic waves radiated (through the corner frequency of waves) and the seismic moment M_{0} (e.g., Eshelby, 1957; Brune, 1970; Baltay et al., 2011, and references therein). If we consider a circular rupture area with radius d_{crack}, the stress drop Δτ is linked with the seismic moment M_{0} through the following equation (Eshelby, 1957).
On the other hand, the seismic moment M_{0} and the moment magnitude M_{w} in terms of the coseismic magnetic field are also related. Hence, the seismic moment is given by
where μ_{sm} is the shear modulus, d the average slip, D the fractal dimension of rock and B_{cs} the coseismic magnetic field; J corresponds to the total electric current density; μ_{m} is the magnetic permeability of the medium; r is the distance to the fault; l_{max} and l_{min} are the radius of the circular rupture area and the smallest microcrack length, respectively. The circular rupture is calculated using ${l}_{max}=\sqrt{S/\mathit{\pi}}$, where S corresponds to the total rupture area.
In this case, the rupture geometry is circular in both formulations, thus, ${d}_{\mathrm{crack}}={l}_{max}$. Substituting this into Eqs. (13) and (12) we obtain
Equation (14) relates the stress drop with the coseismic magnetic field, seismic rupture, and the electrical and mechanical properties of rocks (lithosphere). We can use the data from the Maule earthquake in order to compare the result of Eq. (14) with those found by other researchers. If we use the fault values ${\mathit{\mu}}_{\mathrm{sm}}=\mathrm{3.3}\times {\mathrm{10}}^{\mathrm{10}}$ Pa, d=4 m and $S=\mathrm{450}\times \mathrm{120}$ km^{2} (Vigny et al., 2011; Yue et al., 2014); the granite rock and brittle properties ${\mathit{\mu}}_{\mathrm{m}}=\mathrm{13.5}\times {\mathrm{10}}^{-\mathrm{7}}$ N A^{−2} (Scott, 1983), $J=\mathrm{5}\times {\mathrm{10}}^{-\mathrm{6}}$ A m^{−2} (Tzanis and Vallianatos, 2002), ${l}_{min}={\mathrm{10}}^{-\mathrm{3}}$ m (Shah, 2011) and D=2.6 (Turcotte, 1997); and the magnetic data B_{cs}≈0.1 nT at r≈250 km (Fig. 5 in Venegas-Aravena et al., 2019), we obtain a stress drop Δτ≈3.4 MPa. This result is in close agreement with the result of Luttrell et al. (2011; 4 MPa). Using this value, we can calculate the elastic shear stress as
where H(t−t_{0}) corresponds to the step function centered at t_{0} (the moment of the earthquake occurrence) and γ_{2} is a dimensionless second step function that represents the fault area where the stress drop exists. This means that γ_{2}=0 outside the rupture area and γ_{2}=1 if the point is within the rupture area. Combining this result into Eq. (10), we calculate the total friction coefficient ${\stackrel{\mathrm{\u203e}}{\mathit{\mu}}}_{\mathrm{T}}$ of the fault as
This total friction coefficient ${\stackrel{\mathrm{\u203e}}{\mathit{\mu}}}_{\mathrm{T}}$ is especially relevant because of the dependence of the coseismic magnetic field B_{cs} (through the stress drop Δτ) and the magnetic anomalies (through the relation $B\propto \dot{\mathit{\sigma}}$). Furthermore, Eq. (16) explains the spatial distribution of friction along the fault in addition to its time variations. In Eq. (16), the spatial changes in friction (represented by γ) are not necessarily related to the seismic rupture area (represented by γ_{2}). However, in the case where they are related, one expects γ_{2} to be a function of γ. That is γ_{2}=γ_{2}(γ(x)). In the general case the total friction coefficient gradient can be written as
where $\mathrm{\Gamma}(\mathit{\gamma},{\mathit{\gamma}}_{\mathrm{2}})=\mathrm{1}/({N}_{\mathrm{0}}-\mathit{\gamma}B{)}^{\mathrm{2}}[{\mathit{\gamma}}_{\mathrm{2}}^{\prime}(\mathit{\gamma}B-{N}_{\mathrm{0}})-{\mathit{\gamma}}^{\prime}{\mathit{\gamma}}_{\mathrm{2}}B]$, $\mathrm{\nabla}\stackrel{\mathrm{\u203e}}{\mathit{\mu}}$ is the friction gradient already defined in Eq. (11), and the same definitions for A and B are used. The second term of Eq. (17) implies that a more complex spatial friction distribution on the fault is expected after the earthquake. When no brittle–plastic contribution is considered (γ=0), the friction is only proportional to the gradient of fault rupture distribution ($\mathrm{\nabla}{\stackrel{\mathrm{\u203e}}{\mathit{\mu}}}_{\mathrm{T}}=-\mathrm{\Delta}\mathit{\tau}H(t-{t}_{\mathrm{0}}){\mathit{\gamma}}_{\mathrm{2}}^{\prime}/{N}_{\mathrm{0}}$). When γ_{2}≈γ, the two terms of the total friction coefficient will be proportional to the spatial distribution gradient ($\mathrm{\nabla}{\stackrel{\mathrm{\u203e}}{\mathit{\mu}}}_{\mathrm{T}}\propto {\mathit{\gamma}}^{\prime}$). If the earthquake does not occur, the second term vanishes and Eq. (11) is recovered.
As an illustration, if we consider one point affected by the rupture area and the same values as in Fig. 3, we can calculate the shape of the total friction coefficient as shown in Fig. 5. The three cases show an increase in the friction coefficient prior to the earthquake, and in the three cases the friction reaches its maximum values before the earthquake (at t=10 a.u.). The subsequent decrease is due to the stress drop influence calculated using the coseismic magnetic field B_{cs}. After the earthquake, in none of the cases is the friction completely recovered to its values before the earthquake.
As a second illustration, it is also possible to compare the temporal behavior of the friction coefficient at different points within the fault. For instance, in Fig. 6 (top) we show the different friction behaviors expected if we consider the total friction coefficient. The red area indicates the area where the friction drop occurs. A friction drop is not expected to occur outside the red zone. However, it is possible to observe differences in the coefficient of friction close to the rupture zone where γ_{2}≠γ (yellow area). In both the rupture zone and its surroundings (marked as red and yellow areas, respectively, in Fig. 6) there are changes in the coefficient of friction prior to the rupture. In stark contrast, the standard friction drop is expected only in the rupture zone (red area) and a nonmeasurable change in friction is expected prior to the earthquake (Fig. 6, bottom).
We can summarize our analysis by using two different distributions of γ and γ_{2}. For instance, in Fig. 7 we show the double-sigmoidal distribution for γ and γ_{2} (black and red curves, respectively) along the fault x direction (of total length 2x_{half}). The dimensionless distribution is defined as a combination of the sigmoidal function used in Sect. 2 as
If we consider $a=b=\mathrm{1}$ and x_{half}=10 a.u. in both distributions; x_{0}=5 a.u., L=10 a.u. and w=1 a.u. for γ; and x_{0}=8 a.u., L=4 a.u. and w=10 a.u. for γ_{2}, we get the two distributions of Fig. 7. These values indicate the same situation as discussed in Fig. 6, that is a rupture length (L=4 a.u. in γ_{2} represents the x direction of the red area in Fig. 6) influenced less by the friction coefficient than by semibrittle–plastic stress (here L=10 a.u. in γ represents the x direction of the yellow area in Fig. 6). Using both sets of values, those used in the stress drop (Eq. 14) and also the same k parameters used in Eq. (9), we calculate the total friction coefficient (Eq. 16) as shown in Fig. 8 (case ${k}_{\mathrm{2}}={k}_{\mathrm{4}}={\mathrm{10}}^{-\mathrm{9}}$). At time t=0 a.u., no friction changes occur (${\stackrel{\mathrm{\u203e}}{\mathit{\mu}}}_{\mathrm{T}}=\mathrm{0.5}$). However, the friction increases at t=5 a.u., where the spatial distribution γ is initially defined as nonzero ($x\in [\mathrm{5},\mathrm{15}]$ a.u.). The friction increases up to t=10 a.u., the time corresponding to the earthquake occurrence. The earthquake rupture length is shown as a sudden friction decrease ($x\in [\mathrm{8},\mathrm{12}]$ a.u.) from 0.76 to 0.67. In the zone immediately next to the rupture, the friction increases even more up to the maximum values (0.77), while in the rupture section the friction decreases. After this time (t∼12 a.u.), the rupture and the surrounding section experience a friction decrease.
The case of ${k}_{\mathrm{2}}={k}_{\mathrm{4}}={\mathrm{10}}^{-\mathrm{10}}$ is shown in Fig. 9. The rupture is shown as a blue area at t=10 a.u. in section $x\in [\mathrm{8},\mathrm{12}]$ a.u. At this time and location, we observe that the friction decreases to similar initial values (∼0.5). The friction of this rupture section increases after the earthquake up to ∼0.6. On the other hand, the rupture's surrounding section increases up to the maximum values (∼0.7). Despite those facts, the initial (t=0 a.u.) and final (t=20 a.u.) friction values are almost identical for both cases. For instance, the rupture area has values close to 0.6 in Figs. 8 and 9. The surrounding rupture section has values close to 0.69 in both cases, and the section away from the rupture (close to x=0 a.u. and x=20 a.u.) always has the same initial value (0.5) in both cases. However, both cases exhibit a complex time behavior after the earthquake occurrence, as Eq. (17) indicates.
Finally when the semibrittle–plastic contribution is not considered, we get that γ=0. The result is shown in Fig. 10 and is only observed with the sudden friction decrease at t=10 a.u. (and $x\in [\mathrm{8},\mathrm{12}]$ a.u.). In this later case none of the other complex friction behaviors is observed.
A general expression for friction has been obtained in Eq. (16). Following this we can adapt the Gutenberg–Richter law in this semibrittle–plastic context. Let us remind ourselves that this law establishes a relation between the number of earthquakes in a region and their magnitude (Gutenberg and Richter, 1944). Mathematically, this link is written as $\mathrm{log}N=a-b\times m$, where N is the number of earthquakes with magnitude equal to or greater than m and a and b are constants. However, the b-value of this law might evolve in time. Indeed, De Santis et al. (2011) determined an equation that relates the temporal evolution of the parameter b (b=b(t)) and the measure of the entropy (stress) H(t) of that region as $b\left(t\right)={b}_{max}{\mathrm{10}}^{-H\left(t\right)}$, where ${b}_{max}=e{\mathrm{log}}_{\mathrm{10}}e$. The entropy H(t) is assumed to be proportional to the real stress in the lithosphere (Venegas-Aravena et al., 2019). If we consider the total shear stress τ_{T} in the fault described by Eq. (16), we can write
where ${\mathit{\tau}}_{\mathrm{T}}={\mathit{\tau}}_{\mathrm{0}}-{\mathit{\gamma}}_{\mathrm{2}}\mathrm{\Delta}\mathit{\tau}H(t-{t}_{\mathrm{0}})+\mathit{\gamma}({k}_{\mathrm{1}}\dot{\mathit{\sigma}}+{k}_{\mathrm{2}}\ddot{\mathit{\sigma}}\mathit{\delta}t)\mathrm{cos}\mathit{\theta}\mathit{\delta}t$ and k_{0} is a constant with dimensions of an inverse stress. If we use the same values as in Fig. 5 (${k}_{\mathrm{2}}={k}_{\mathrm{4}}={\mathrm{10}}^{-\mathrm{9}}$), we get the temporal evolution of the b-value (Fig. 11a). Figure 11a shows a decrease in the b-value until the earthquake occurs (at t=10 a.u.). Figure 11b shows the Gutenberg–Richter law at three instants of time (and k_{0}=0.01, a=1): initial (t=0 a.u.), prior to the earthquake (t=9 a.u.) and final (t=20 a.u.). This figure shows how large earthquakes (M_{w}∼6) are not expected at the initial moment (blue line). However, just before the earthquake, an M8-class earthquake should be expected (green line). After the earthquake one would only expect earthquakes no greater than M_{w}∼7 to exist (red line). Figure 12a and b show the same previous case but now considering ${k}_{\mathrm{2}}={k}_{\mathrm{4}}={\mathrm{10}}^{-\mathrm{10}}$. In this case there is no difference between time immediately before the earthquake (t=9 a.u.) and at the final time (t=20 a.u.). In addition, using these parameters (${k}_{\mathrm{2}}={k}_{\mathrm{4}}={\mathrm{10}}^{-\mathrm{10}}$), smaller magnitudes are reached than when using ${k}_{\mathrm{2}}={k}_{\mathrm{4}}={\mathrm{10}}^{-\mathrm{9}}$ (M_{w}∼7 and M_{w}∼8, respectively). This is why it is more likely that earthquakes of greater magnitude will be found when considering the contribution of the $\ddot{\mathit{\sigma}}$ term in the analysis.
In Sect. 5 we determine that ${k}_{\mathrm{2}}={k}_{\mathrm{4}}={\mathrm{10}}^{-\mathrm{9}}$ were the most adequate values. Despite this fact, the Gutenberg–Richter law does not provide an approximate time t_{0} for the earthquake occurrences. If we look at Eq. (16), the term t_{0} appears explicitly (in the step function). However, t_{0} appears only after the earthquake, so it is not possible to find t_{0} analytically before the rupture (using the rupture itself). This means that we must find an estimate using the other parameters. One way to do it is to consider the differential total friction coefficient $\mathrm{d}{\stackrel{\mathrm{\u203e}}{\mathit{\mu}}}_{\mathrm{T}}$ to find an approximate rupture time t_{0}. Figure 13 shows $\mathrm{d}{\stackrel{\mathrm{\u203e}}{\mathit{\mu}}}_{\mathrm{T}}$ considering ${\mathit{\gamma}}_{\mathrm{2}}=\mathit{\gamma}=\mathrm{1}$, ${k}_{\mathrm{2}}={k}_{\mathrm{4}}={\mathrm{10}}^{-\mathrm{9}}$ and ${k}_{\mathrm{2}}={k}_{\mathrm{4}}={\mathrm{10}}^{-\mathrm{10}}$. When major earthquakes (${k}_{\mathrm{2}}={k}_{\mathrm{4}}={\mathrm{10}}^{-\mathrm{9}}$) are considered, the rupture occurs after the maximum value ($(\mathrm{d}{\stackrel{\mathrm{\u203e}}{\mathit{\mu}}}_{\mathrm{T}}{)}_{\mathrm{MAX}}$), when $\mathrm{d}{\stackrel{\mathrm{\u203e}}{\mathit{\mu}}}_{\mathrm{T}}\approx \mathrm{1}/\mathrm{2}(\mathrm{d}{\stackrel{\mathrm{\u203e}}{\mathit{\mu}}}_{\mathrm{T}}{)}_{\mathrm{MAX}}$ (Fig. 13, top). When k_{2}, ${k}_{\mathrm{4}}={\mathrm{10}}^{-\mathrm{10}}$, the rupture also occurs after the maximum value, in this case $\mathrm{d}{\stackrel{\mathrm{\u203e}}{\mathit{\mu}}}_{\mathrm{T}}\approx \mathrm{0.9}(\mathrm{d}{\stackrel{\mathrm{\u203e}}{\mathit{\mu}}}_{\mathrm{T}}{)}_{\mathrm{MAX}}$ (Fig. 13, bottom). Considering these two limiting cases in this framework we conclude that earthquakes occur after ($\mathrm{d}{\stackrel{\mathrm{\u203e}}{\mathit{\mu}}}_{\mathrm{T}}{)}_{\mathrm{MAX}}$ when the differential total friction coefficient decreases to $\mathrm{d}{\stackrel{\mathrm{\u203e}}{\mathit{\mu}}}_{\mathrm{T}}\approx C(\mathrm{d}{\stackrel{\mathrm{\u203e}}{\mathit{\mu}}}_{\mathrm{T}}{)}_{\mathrm{MAX}}$, where $C\in [\mathrm{0.5},\mathrm{0.9}]$.
The time between $\mathrm{d}{\stackrel{\mathrm{\u203e}}{\mathit{\mu}}}_{\mathrm{T}}\approx (\mathrm{d}{\stackrel{\mathrm{\u203e}}{\mathit{\mu}}}_{\mathrm{T}}{)}_{\mathrm{MAX}}$ and $\mathrm{d}{\stackrel{\mathrm{\u203e}}{\mathit{\mu}}}_{\mathrm{T}}\approx C(\mathrm{d}{\stackrel{\mathrm{\u203e}}{\mathit{\mu}}}_{\mathrm{T}}{)}_{\mathrm{MAX}}$ can be denoted by δ. This δ parameter increases when C decreases and vice versa, so δ is inversely proportional to C ($\mathit{\delta}\propto {C}^{-\mathrm{1}}$). Then, we arrive at an expression for the general rupture time t_{0} as follows:
where ${t}_{(\mathrm{d}{\stackrel{\mathrm{\u203e}}{\mathit{\mu}}}_{\mathrm{T}}{)}_{\mathrm{MAX}}}$ is the time when $\mathrm{d}{\stackrel{\mathrm{\u203e}}{\mathit{\mu}}}_{\mathrm{T}}\approx (\mathrm{d}{\stackrel{\mathrm{\u203e}}{\mathit{\mu}}}_{\mathrm{T}}{)}_{\mathrm{MAX}}$. Note that Eq. (20) is only valid after the maximum value ($\mathrm{d}{\stackrel{\mathrm{\u203e}}{\mathit{\mu}}}_{\mathrm{T}}{)}_{\mathrm{MAX}}$ is reached. Equation (20) is general; however, considering the Gutenberg–Richter law, we expect that C values close to 0.5 are necessary to represent earthquakes of greater magnitude in this theoretical framework.
In this work we have linked key parameters associated with the geological fault with magnetic measurements. Both stress drops and the semibrittle–plastic stress were linked to a friction coefficient (on the fault) equation in terms of magnetic measurements (Eq. 16). One of the main points of Eq. (16) corresponds to the fact that the total friction coefficient ${\stackrel{\mathrm{\u203e}}{\mathit{\mu}}}_{\mathrm{T}}(x,t)$ is entirely determined by the spatial distribution of the nonconstant stress changes within the lithosphere. If the seismo-electromagnetic theory is applied, it means that the rupture process is controlled by the nonconstant stress changes that surround the fault and not by the fault itself. In this scenario, the earthquakes might occur at places on the fault that are being affected by a continuous friction increase prior to rupture (this friction increase occurs regardless of the values of k_{2} and k_{4} used in the analysis). In our work we have shown that the total friction coefficient depends on two different spatial distributions. The first one is associated with the uniaxial stress changes γ(x) and the second one with the rupture area γ_{2}(x). These two distributions are not necessarily correlated. In the case where they are comparable (that is $\left|\mathit{\gamma}\right(x)-{\mathit{\gamma}}_{\mathrm{2}}(x\left)\right|\approx \mathrm{0}$), it means that the lithosphere area affected by nonconstant uniaxial stress changes is a determining factor for predicting the earthquake magnitude and location before it occurs. This is $\mathit{\gamma}\left(x\right)\approx {\mathit{\gamma}}_{\mathrm{2}}\left(x\right)\approx \mathrm{0}$ if x belong to sections where $\dot{\mathit{\sigma}}=\mathrm{cte}$ and γ(x) and γ_{2}(x)≠0 if x∈L, where L is the rupture length (at places where $\dot{\mathit{\sigma}}\ne \mathrm{cte}$). The seismic moment M_{0} is proportional to this length L (Aki, 1966), and the seismic moment magnitude M_{w} depends on the seismic moment (Hanks and Kanamori, 1979), implying that the seismic moment magnitude depends on the spatial distribution of the total friction coefficient variations. As the nonconstant uniaxial stress changes could also create magnetic signals due the microcracks of rocks (Venegas-Aravena et al., 2019, and references therein), it is reasonable to assume that a larger area of magnetic anomalies could imply a larger earthquake. This is the locality (or Dobrovolsky) criterion used by some researchers when they relate some electromagnetic measurements to earthquakes (e.g., De Santis et al., 2019a, and references therein). If in this case we also consider the initial time of friction increase (impending rupture time), the approximate magnitude, then the approximate location and the approximate imminent time could be theoretically determined through Eq. (20).
For the limiting case where $\left|\mathit{\gamma}\right(x)-{\mathit{\gamma}}_{\mathrm{2}}(x\left)\right|\gg \mathrm{0}$ we have seen that the locality criterion does not hold anymore. Hence, the earthquakes occurrences cannot be related to the nonconstant stress changes and the magnetic field. This limiting case prevents the possibility of a real earthquake prediction using this simple theoretical base. This later case may also imply that the cumulative stress on the fault is not enough to generate a seismic rupture at any point on the fault (that is γ_{2}(x)=0). This means that the accumulated elastic energy that is injected into the fault is not sufficient to spread the rupture (this last energy is called fracture energy G_{C}; see Ohnaka, 2013; Nielsen et al., 2016). This scenario could indicates that the stress changes in the semibrittle–plastic regime are not a sufficient condition for the earthquake generation; however, they could be a necessary condition.
With respect to the size of earthquakes in the present model, Sect. 5 has shown that earthquakes have greater magnitudes when ${k}_{\mathrm{2}}={k}_{\mathrm{4}}\approx {\mathrm{10}}^{-\mathrm{9}}$. If we consider those values in Fig. 5, we see that the total friction coefficient is also higher. This is ∼0.75 when ${k}_{\mathrm{2}}={k}_{\mathrm{4}}\approx {\mathrm{10}}^{-\mathrm{9}}$ and sensibly larger than ∼0.6 when ${k}_{\mathrm{2}}={k}_{\mathrm{4}}\approx {\mathrm{10}}^{-\mathrm{10}}$. This indicates that there is a correlation between the size of the earthquake and the total friction coefficient. Hence, the earthquake has a greater magnitude when there is a higher total friction coefficient (or shear stress τ). This means that ${\mathit{\gamma}}_{\mathrm{2}}={\mathit{\gamma}}_{\mathrm{2}}\left({\stackrel{\mathrm{\u203e}}{\mathit{\mu}}}_{\mathrm{T}}\right)$; therefore, the rupture length L of Eq. (18) is proportional to ${\stackrel{\mathrm{\u203e}}{\mathit{\mu}}}_{\mathrm{T}}$ (that is $L=\mathit{\alpha}{\stackrel{\mathrm{\u203e}}{\mathit{\mu}}}_{\mathrm{T}}\left({t}_{\mathrm{0}}\right)$). Note that this is independent of the value of $\left|\mathit{\gamma}\right(x)-{\mathit{\gamma}}_{\mathrm{2}}(x\left)\right|$, since it comes directly from the Gutenberg–Richter law. However, in this case ${\stackrel{\mathrm{\u203e}}{\mathit{\mu}}}_{\mathrm{T}}$ is homogeneous, so more studies will be needed in order to determine L.
Finally, this theoretical work has provided a possible mechanism that explains several magnetic measurements performed recently. Equally this work provides some necessary conditions of the fault in order to trigger earthquakes in terms of their magnetic properties. We think that future research of the LAIC effect community should focus on lithospheric-fault dynamics as one of the promising topics in the field. When the lithosphere part of this effect is better understood, the other effects will have a stronger theoretical basis that will help to perform measurements and/or make accurate predictions.
All the data are open source and can be found using the references that are listed in the text. The numerical data can be easily generated by everyone using the equations and indications of the text. If you have any problems, do not hesitate to write and ask the authors.
PVA conceived the theoretical derivation and numerical implementation. PVA also performed the initial paper writing and created the figures. EGC and DL carried out the writing corrections and generated the scientific view, background and support.
The authors declare that they have no conflict of interest.
This paper is in memory of Marcela Larenas Clerc (1960–2019). Patricio Venegas-Aravena acknowledges Patricia Aravena, Alejandro Venegas, Patricia Venegas and Richard Sandoval for their outstanding support in carrying out this work and Valeria Becerra-Carreño for her scientific support. David Laroze acknowledges partial financial support from Centers of Excellence with CONICYT Basal financing, grant AFB180001, CEDENNA.
This research has been partially supported by Centers of Excellence with CONICYT Basal (grant no. AFB180001, CEDENNA).
This paper was edited by Filippos Vallianatos and reviewed by Angelo De Santis and Michael E. Contadakis.
Aki, K.: Generation and propagation of G waves from the Niigata earthquake of June 14, 1964. Part 2. Estimation of earthquake moment, released energy and stress-strain drop from G wave spectrum, B. Earthq. Res. I., 44, 73–88, 1966.
Anastasiadis, C., Triantis, D., Stavrakas, I. and Vallianatos, F.: Pressure Stimulated Currents (PSC) in marble samples, Ann. Geophys., 47, 21–28, 2004.
Baltay, A., Ide, S., Prieto, G., and Beroza, G.: Variability in earthquake stress drop and apparent stress, Geophys. Res. Lett., 38, L06303, https://doi.org/10.1029/2011GL046698, 2011.
Brune, J. N.: Tectonic stress and the spectra of seismic shear waves from earthquakes, J. Geophys. Res., 75, 4997–5009, https://doi.org/10.1029/JB075i026p04997, 1970.
Byerlee, J. D.: Friction of Rocks, Pure Appl. Geophys., 116, 615–626, https://doi.org/10.1007/BF00876528, 1978.
Chen, G. S.: Fundamentals of contact mechanics and friction, Handbook of Friction-Vibration Interactions, 1st edn., Woodhead Publishing, 71–152, https://doi.org/10.1533/9780857094599.71, 2014.
Cordaro, E. G., Venegas, P., and Laroze, D.: Latitudinal variation rate of geomagnetic cutoff rigidity in the active Chilean convergent margin, Ann. Geophys., 36, 275–285, https://doi.org/10.5194/angeo-36-275-2018, 2018.
De Santis, A., Cianchini, G., Favali, P., Beranzoli, L., and Boschi, E.: The Gutenberg–Richter Law and Entropy of Earthquakes: Two Case Studies in Central Italy, B. Seismol. Soc. Am., 101, 1386–1395, https://doi.org/10.1785/0120090390, 2011.
De Santis, A., Balasis, G., Pavón-Carrasco, F. J., Cianchini, G., and Mandea, M.; Potential earthquake precursory pattern from space: The 2015 Nepal event as seen by magnetic Swarm satellites, Earth Planet. Sc. Lett., 461, 119–126, https://doi.org/10.1016/j.epsl.2016.12.037, 2017.
De Santis, A., Abbattista, C., Alfonsi, L., Amoruso, L., Campuzano, S. A., Carbone, M., Cesaroni, C., Cianchini, G., De Franceschi, G., De Santis, A., Di Giovambattista, R., Marchetti, D., Martino, L., Perrone, L., Piscini, A., Rainone, M. L., Soldani, M., Spogli, L., and Santoro, F.: Geosystemics View of Earthquakes, Entropy, 21, 412, https://doi.org/10.3390/e21040412, 2019a.
De Santis, A., Marchetti, D., Spogli, L., Cianchini, G., Pavón-Carrasco, F. J., Franceschi, G. D., Di Giovambattista, R., Perrone, L., Qamili, E., Cesaroni, C., De Santis, A., Ippolito, A., Piscini, A., Campuzano, S. A., Sabbagh, D., Amoruso, L., Carbone, M., Santoro, F., Abbattista, C., and Drimaco, D.: Magnetic Field and Electron Density Data Analysis from Swarm Satellites Searching for Ionospheric Effects by Great Earthquakes: 12 Case Studies from 2014 to 2016, Atmosphere, 10, 371, https://doi.org/10.3390/atmos10070371, 2019b.
Dobrovolsky, I. R., Zubkov, S. I., and Myachkin, V. I.: Estimation of the size of earthquake preparation zones, Pure Appl. Geophys., 117, 1025–1044, https://doi.org/10.1007/BF00876083, 1979.
Eshelby, J. D.: The determination of the elastic field of an ellipsoidal inclusion, and related problems, P. Roy. Soc. Lond. A Mat., A241, 376–396, https://doi.org/10.1098/rspa.1983.0054, 1957.
Freund, F.: Rocks That Crackle and Sparkle and Glow: Strange Pre-Earthquake Phenomena, Journal of Scientic Exploration, 17, 37–71, 2003.
Gutenberg, B. and Richter, C. F.: Frequency of earthquakes in California, B. Seismol. Soc. Am., 34, 185–188, 1944.
Hanks, T. C. and Kanamori, H.: A moment magnitude scale, J. Geophys. Res., 84, 2348–2350, https://doi.org/10.1029/JB084iB05p02348, 1979.
Kyriazopoulos, A., Anastasiadis, C., Triantis, D., and Brown, C. J.: Non-destructive evaluation of cement-based materials from pressure-stimulated electrical emission – Preliminary results, Constr. Build. Mater., 25, 1980–1990, https://doi.org/10.1016/j.conbuildmat.2010.11.053, 2011.
Lamb, S.: Shear stresses on megathrusts: Implications for mountain building behind subduction zones, J. Geophys. Res., 111, B07401, https://doi.org/10.1029/2005JB003916, 2006.
Luttrell, K. M., Tong, X., Sandwell, D. T., Brooks, B. A., and Bevis, M. G.: Estimates of stress drop and crustal tectonic stress from the 27 February 2010 Maule, Chile, earthquake: Implications for fault strength, J. Geophys. Res., 116, B11401, https://doi.org/10.1029/2011JB008509, 2011.
Maksymowicz, A.: The geometry of the Chilean continental wedge: Tectonic segmentation of subduction processes off Chile, Tectonophysics, 659, 183–196, https://doi.org/10.1016/j.tecto.2015.08.007, 2015.
Marchetti, D. and Akhoondzadeh, M.: Analysis of Swarm satellites data showing seismo-ionospheric anomalies around the time of the strong Mexico (M_{w}=8.2) earthquake of 08 September 2017, Adv. Space Res., 62, 614–623, https://doi.org/10.1016/j.asr.2018.04.043, 2018.
Marone, C. J. and Saffer, D. M.: The Mechanics of Frictional Healing and Slip Instability During the Seismic Cycle, Treatise on Geophysics: Second Edition, 4, 111–138, https://doi.org/10.1016/B978-0-444-53802-4.00092-0, 2015.
Nielsen, S., Spagnuolo, E., Violay, M., Smith, S., Di Toro, G., and Bistacchi, A.: G: Fracture energy, friction and dissipation in earthquakes, J. Seismol., 20, 1187–1205, https://doi.org/10.1007/s10950-016-9560-1, 2016.
Ohnaka, M.: The Physics of Rock Failure and Earthquakes, Cambridge University Press, Cambridge, UK, https://doi.org/10.1017/CBO9781139342865.007, 2013.
Papangelo, A., Ciavarella, M., and Barber, J. R.: Fracture mechanics implications for apparent static friction coefficient in contact problems involving slip-weakening laws471, P. R. Soc. A, 471, 20150271, https://doi.org/10.1098/rspa.2015.0271, 2015.
Parlitz, U., Hornstein, A., Engster, D., Al-Bender, F., Lampaert, V., Tjahjowidodo, T., Fassois, S. D., Rizos, D., Wong, C. X., Worden, K., and Manson, G.: Identification of pre-sliding friction dynamics, Chaos, 14, 420–430, https://doi.org/10.1063/1.1737818, 2004.
Potirakis, S. M., Contoyiannis, Y., Asano, T., and Hayakawa, M.: Intermittency-induced criticality in the lower ionosphere prior to the 2016 Kumamoto earthquakes as embedded in the VLF propagation data observed at multiple stations, Tectonophysics, 722, 422–431, 2018a.
Potirakis, S. M., Asano, T., and Hayakawa, M.: Criticality Analysis of the Lower Ionosphere Perturbations Prior to the 2016 Kumamoto (Japan) Earthquakes as Based on VLF Electromagnetic Wave Propagation Data Observed at Multiple Stations, Entropy, 20, 199, https://doi.org/10.3390/e20030199, 2018b.
Saltas, V., Vallianatos, F., Triantis, D., and Stavrakas, I.: Complexity in laboratory seismology: From electrical and acoustic emissions to fracture, in: Complexity of seismic time series; measurement and Application, edited by: Chelidze, T., Telesca, L., and Vallianatos, F., Elsevier, Amsterdam, the Netherlands, 239–273, https://doi.org/10.1016/B978-0-12-813138-1.00008-0, 2018.
Schekotov, A. and Hayakawa, M.: Seismo-meteo-electromagnetic phenomena observed during a 5-year interval around the 2011 Tohoku earthquake, Phys. Chem. Earth, 85–86, 167–173, https://doi.org/10.1016/j.pce.2015.01.010, 2015.
Scott, J. H.: Electrical and Magnetic properties of rock and soil,USGS Publications Warehouse, USGS Open-File Report 83-915, https://doi.org/10.3133/ofr83915, 1983.
Shah, K. P.: The Hand Book on Mechanical Maintenance, Practical Maintenance, compiled by: Shah, K. P., available at: http://practicalmaintenance.net/?p=1135 (last access: 19 May 2020), 2011.
Sun, N. P. and Mosleh, M.: The Minimum Coefficient of Friction: What Is It?, CIRP Annals, 43, 491–495, https://doi.org/10.1016/s0007-8506(07)62260-4, 1994.
Turcotte, D. L.: Fractals and Chaos in Geology and Geophysics, 2nd edn., Cambridge University Press, Cambridge, UK, 397 pp., 1997.
Tzanis, A. and Vallianatos, F.: A physical model of electrical earthquake precursors due to crack propagation and the motion of charged edge dislocations, in: Seismo Electromagnetics: Lithosphere–Atmosphere–Ionosphere Coupling, edited by: Hayakawa, M. and Molchanov, O. A., TERRAPUB, Tokyo, Japan, 117–130, 2002.
Utada, H., Shimizu, H., Ogawa, T., Maeda, T., Furumura, T., Yamamoto, T., Yamazaki, N., Yoshitake, Y., and Nagamachi, S.: Geomagnetic field changes in response to the 2011 off the Pacific Coast of Tohoku earthquake and tsunami, Earth Planet. Sc. Lett., 311, 11–27, https://doi.org/10.1016/j.epsl.2011.09.036, 2011.
Vallianatos, F. and Triantis, D.: Scaling in Pressure Stimulated Currents related with rock fracture, Physica A, 387, 4940–4946, https://doi.org/10.1016/j.physa.2008.03.028, 2008.
Vallianatos, F. and Tzanis, A.: Electric Current Generation Associated with the Deformation Rate of a Solid: Preseismic and Coseismic Signals, Phys. Chem. Earth, 23, 933–938, 1998.
Vallianatos, F. and Tzanis, A.: On the nature, scaling and spectral properties of pre-seismic ULF signals, Nat. Hazards Earth Syst. Sci., 3, 237–242, https://doi.org/10.5194/nhess-3-237-2003, 2003.
Venegas-Aravena, P., Cordaro, E. G., and Laroze, D.: A review and upgrade of the lithospheric dynamics in context of the seismo-electromagnetic theory, Nat. Hazards Earth Syst. Sci., 19, 1639–1651, https://doi.org/10.5194/nhess-19-1639-2019, 2019.
Vigny, C., Socquet, A., Peyrat, S., Ruegg, J.-C., Metois, M., Madariaga, R., Morvan, S., Lancieri, M., Lacassin, R., Campos, J., Carrizo, D., Bejar-Pizarro, M., Barrientos, S., Armijo, R., Aranda, C., Valderas-Bermejo, M.-C., Ortega, I., Bondoux, F., Baize, S., Lyon-Caen, H., Pavez, A., Vilotte, J. P., Bevis, M., Brooks, B., Smalley, R., Parra, H., Baez, J.-C., Blanco, M., Cimbaro, S., and Kendrick, E.: The 2010 M_{w} 8.8 Maule Megathrust Earthquake of Central Chile, monitored by GPS, Science, 332, 1417–1421, 2011.
Yue, H., Lay, T., Rivera, L., An, C., Vigny, C., Tong, X., and Báez Soto, J. C.: Localized fault slip to the trench in the 2010 Maule, Chile M_{w}=8.8 earthquake from joint inversion of high-rate GPS, teleseismic body waves, InSAR, campaign GPS, and tsunami observations, J. Geophys. Res.-Sol. Ea., 119, 7786–7804, https://doi.org/10.1002/2014JB011340, 2014.
Zhang, H., Sun, Q., and Ge, Z: Analysis of the characteristics of magnetic properties change in the rock failure process, Acta Geophys. 68, 289–302, https://doi.org/10.1007/s11600-020-00406-3, 2020.
An interactive open-access journal of the European Geosciences Union