Natural Fractals as Irreversible Disorder: Entropy Approach from Cracks in the Semi Brittle-Ductile Lithosphere and Generalization

The seismo-electromagnetic theory describes the growth of fractally distributed cracks within the lithosphere that generate the emission of magnetic anomalies prior to large earthquakes. One of the main physical properties of this theory is their consistency regarding the second law of thermodynamics. That is, the crack generation of the lithosphere corresponds to the manifestation of an irreversible process evolving from one steady state to another. Nevertheless, there is still not a proper thermodynamic description of lithospheric crack generation. That is why this work presents the derivation of the entropy changes generated by the lithospheric cracking. It is found that the growth of the fractal cracks increases the entropy prior impending earthquakes. As fractality is observed across different topics, our results are generalized by using the Onsager’s coefficient for any system characterized by fractal volumes. It is found that the growth of fractality in nature corresponds to an irreversible process.


Introduction
Earthquakes are complex natural phenomena that can be studied by using different approaches. Some of them are focused on fault mechanics. That is, the dislocation and rupture of preexisting faults within the lithosphere which generate seismic radiation [1]. Others are focused on the geodetic deformation of the Earth's surface due the tectonic plates' drift which generate long-term stress accumulation [2]. Finally, a recent focus corresponds to the electromagnetic signals that can be linked to crack and microcrack generation within the lithosphere which could be considered as the manifestation of impending earthquakes [3].
The first approach makes impossible the forecast of major earthquakes due to two main reasons: the vanishing seismic source information due to the attenuation of seismic waves [4] and, the interplay among different heterogeneous physical processes that make the rupture chaotic [5]. Despite this, it has been recently shown that any prediction should be done using multidisciplinary precursors [6]. In that sense, the seismo-electromagnetic approach could be considered as multidisciplinary because is based on experimental studies on brittle lithospheric rock samples, magnetic data and physical analysis. These studies can be summarized with the following five points: (1) the volume of cracks increases before the macroscopic failure [7,8], (2) the external mechanical loads generate detectable electromagnetic signals known as pressure stimulated currents (PSC) [9], (3) electromagnetic signals generated from cracks are well described by the dislocation process known as motion of charged edge dislocations (MCD) [10][11][12][13], (4) the number of cumulative magnetic anomalies near and prior impending megathrust earthquakes also increases [14][15][16], and (5) fractally distributed cracks within the lithosphere could explain and link the above mentioned points in the so-called seismo-electromagnetic theory [17].
The physical basis of the seismo-electromagnetic theory is focused on electromagnetic signals that rise due to pre-failure states such as the material cracking or complex multiparametric statistical approach [18,19]. In addition, the theory also takes into account the concept of "earthquake entropy" which relates the lithospheric stress changes to the b-value of the Gutenberg-Richter's law within faults [20]. That is, a fault's properties give information about stress states within the lithosphere. In other words, stress states within the lithosphere can be linked to frictional properties which could be related to other seismic parameters such as seismic magnitude or seismic moment [21]. As lithospheric stress generates the fractally distributed cracks and earthquakes correspond to an irreversible process [22], it is expected that the physical process that generates fractal cracks also increases the entropy of the lithosphere before impending earthquakes. The entropy change for one single crack is well understood in terms of electrical current, friction, and temperature [23]. Nevertheless, realistic cracks are fractally embedded within materials.
On the other hand, it has been shown that other systems similar to fractal cracks evolving in time are also present in several natural phenomena. For instance, some of those systems are (including seismology studies): • Earthquakes spatial distributions [24], • Earthquake slip patterns [25], • Cracks in rocks and lithosphere [17,26], • Structural geology [27], • Galaxies clustering [28], • Self-organized criticality (SOC) systems [29], • Quantum scale properties [30,31] [43] allows scientists to propose that these systems are governed by fractal laws such as the "Constructal law" [44] or principle of least action [45,46]. Nevertheless, fractals are geometrically well described but, a general description of thermodynamics fluxes, such as Onsager's coefficients [47,48], which generate fractality is still missing. In order to obtain a clearer physical meaning of one specific fractal system, Section 2.1 demonstrates the extent to which the entropy of one single crack can lead to fractally distributed cracks. Section 2.2 describes the relationship between entropy and seismic moment. Section 2.3 generalizes the entropy change for any system characterized by fractal geometry by using the linear nonequilibrium thermodynamics framework as the Onsager's coefficients. Then, Section 2.4 discusses the extent of the results while the conclusion is in Section 3.

Entropy of Fractals Cracks Distribution
The change of entropy ds 0 /dt required for the generation of a single crack is given by [23]: where µ is the coefficient of friction generated by the inner relative displacement of the crack boundaries, N the normal force, v the relative velocity, T the temperature, and X e and J e are the voltage and electrical current, respectively. Here, the first term (right side) corresponds to the frictional heating while the second one corresponds to the electrical production due the electrical imbalance in the semi-brittle plastic regime [10,11,49,50]. The generation of electromagnetic signals prior to main failure has been widely reported in laboratory experiments [3,[51][52][53][54][55] and on a geodynamic scale [16,[56][57][58]. By contrast, temperature changes prior to earthquakes has been poorly supported in [59]. Then, it is expected that the entropy changes are manly driven by electrical charge generation within microcracks. This implies that: Here and after, the temperature is considered as constant. In addition, it has been shown that the electrical current generated by external uniaxial stress change is given by [60]: where Y e f f is the effective Young's modulus and σ u the applied uniaxial stress. The voltage can be obtained by using the Ohm's law for continuum medium. That is where σ e is the electrical conductivity and → E the electric field [61]. The voltage definition 62] allows us to write the voltage as: where ρ e is the electrical resistivity and a is the distance where the electrical current flows within the crack. Using Equations (2)-(4), the entropy change for a single crack is: where C 0 = ρ e cos θ/T 0 . Let us consider that a and J e are parallel, then C 0 is maximum. The parameter a relates the length where the electrical current flows within a single crack, it can also be considered as proportional to the volume occupied by the cracks. That is, a = V/A 0 , where A 0 corresponds to an area of reference. This implies that Equation (5) is ds 0 /dt ∼ V. This equation stands for the volume of a single crack. However, there are several cracks within a macroscopic material under uniaxial stress prior to the failure [7,8].
Then, each contribution should be considered for macroscopic material. On the other hand, it has been shown that the total entropy change of dS corresponds to the sum of all entropy change at different length scale [63]. For example, dS = d nano S + d meso S + d macro S means the nano, meso, and macro scale contribute to the total entropy change [63]. Here, cracks are also observed at different length scale. Specifically, it has been observed that the microcracks are fractally distributed [64,65]. This implies that entropy is a sum of all the volume contributions with different length scale: continuum fractal distribution, the entropy change turns into: dS ∝ dV = V. Then, the total change of entropy dS/dt depends on a fractal volume: where σ is the macroscopic uniaxial stress and V f r is the fractal volume defined as [66]: where D is the fractal dimension of rocks which mainly lies between 2 and 3, l min is the smallest rupture radius considered and l max is the largest crack within the fractal distribution. For simplicity, the largest crack is considered as circular crack [17]. Experimentally, it has been shown in X-ray tomography studies on rock samples that the generation of cracks is dominated by the growth of preexisting fractures before the main failure [7]. This implies that the volume is increasing with time. In other words, l max = l max (t). Then, Equation (6) in (7) implies that the entropy change is proportional to the growth of the largest crack l max (t) within the fractal volume, growth that is generated by the external stress change dσ/dt: where Here, dS ≥ 0 because l max is a distance which is positive by definition. The only manner in which dS = 0 is when dσ/dt ≈ 0 or when l min ≈ 0. The entropy change imposes that dσ/dt = 0 and l min = 0. The former means that there is an input or external force required in order to generate the cracks. In other words, no stress changes, no cracks growth, and no entropy changes in the system. The latter implies that there is a lower boundary for fractality in nature. The relationship between the external forces applied by F that generates the fractal volume V f r and the entropy change dS can be summarized as: Thus, Equations (8) and (9) indicate that the increase of entropy is obtained by the generation of tridimensional fractals in a certain domain due an increasing external force.

Entropy Change in Terms of Spatial Parameters
Let us consider now a brief description of Section 2.1. Equation (8) shows the relationship between entropy and spatial properties as fractal geometry. Let us consider now how the entropy is changing in terms of specific spatial parameters as the maximum length l max and the fractal dimension D while dσ/dt is constant. This can be seen in Figure 1a. The entropy changes dS increases when l max increases regardless of the value of D. This can be seen as the change from blue to red colors in Figure 1a. Nevertheless, the dependency of dS on D is different and also dependent on l max . For example, dS increases when D increases for small values of l max (l max < 1). That can be seen as the increase dS from ∼ −7.3 a.u. to −4.8 a.u. in log scale. Note that a.u. means arbitrary units and the colormap of Figure 1a shows ln(dS) which allows negative values. Contrary, dS increases when D decreases for large values of l max . The latter can be seen in the marked grey rectangle in Figure 1a. That area is characterized by a dark red color and represents large values of dS at values of l max close to 10 3 a.u and a fractal dimension lower than 2.5. The dS values decrease for the same l max and large D values (D > 2.5).
shows that the domain is mainly governed by small cracks while large cracks are rare ( Figure 1c). The opposite is found in (Figure 1b).

Entropy Change in Terms of External Stress Change
Now it is relevant to consider if all the fractally distributed cracks increase the entropy. Let us consider two cases, one where the stress increases linearly and sigmoidshaped in time. The first case is when / and are constant (equal to 1 arbitrary unit or a.u.). This case can be seen in the black curve in Figure 2a. Here, the entropy changes and increases linearly (in a log-log plot) when the fractal dimension for granite ( = 2.6 ) is considered [66]. The second case is when / could be considered as proportional to the sigmoidal shape which has been found for real earthquakes [14][15][16]. That is, ( ) ∼ ln(1 + ( − ) ) and / ∼ (1 + − ( − ) ) −1 [17], where is constant (here equal to 1) and corresponds to the time where the macroscopic failure occurs. Figure 2b shows the evolution of (black curve) and / (red curve) for the sigmoidal case. This input generates the entropy changes shown in Figure 2a (red curve). It is clear that the incorporation of / reduce the entropy change for small to medium values. That is, the red curve is lower than the black one in Figure 2a. Nevertheless, in both cases, the black and red curves reach the same values for large . Note that the vertical purple segmented line in Figure 2b shows when = . This can be seen in Figure 1a as a purple dot in red curve ( Figure 2a). This shows that the incorporation of an increasing stress ( / > 0) increases the entropy change . In other words, the external forces increase the entropy from one state to another. Figure 1b,c shows spheres that represent the randomly distributed fractal cracks for two values of fractal dimension 2.99 and 2.01 respectively when l max is large. Small D shows that the domain is mainly governed by small cracks while large cracks are rare ( Figure 1c). The opposite is found in (Figure 1b).

Entropy Change in Terms of External Stress Change
Now it is relevant to consider if all the fractally distributed cracks increase the entropy. Let us consider two cases, one where the stress increases linearly and sigmoid-shaped in time. The first case is when dσ/dt and Λ are constant (equal to 1 arbitrary unit or a.u.). This case can be seen in the black curve in Figure 2a. Here, the entropy changes dS and l max increases linearly (in a log-log plot) when the fractal dimension for granite (D = 2.6 ) is considered [66]. The second case is when dσ/dt could be considered as proportional to the sigmoidal shape which has been found for real earthquakes [14][15][16].
That is, σ(t) ∼ ln 1 + e a s (t−t C ) and dσ/dt ∼ 1 + e −a s (t−t C ) −1 [17], where a s is constant (here equal to 1) and t C corresponds to the time where the macroscopic failure occurs. Figure 2b shows the evolution of σ (black curve) and dσ/dt (red curve) for the sigmoidal case. This input generates the entropy changes shown in Figure 2a (red curve). It is clear that the incorporation of dσ/dt reduce the entropy change for small to medium l max values. That is, the red curve is lower than the black one in Figure 2a. Nevertheless, in both cases, the black and red curves reach the same dS values for large l max . Note that the vertical purple segmented line in Figure 2b shows when t = t C . This can be seen in Figure 1a as a purple dot in red curve ( Figure 2a). This shows that the incorporation of an increasing / . The purple point shows the and where the earthquake occurs for the sigmoid function. (b) Stress evolution (black curve) and shear stress change / (red curve) of the lithosphere prior to and after the main failure. Here, the earthquake time is = = 5 a. u. Note that a. u. means arbitrary units. (c) Entropy increases by using Equation (10). It is shown that there are two main behaviors: the initial slow increase that lasts up to ∼ 3.5 a. u. (this trend is represented as a black dotted line) and the fast increases between 3.5 a. u. and = 5 a.u. (d) Magnitude expected in terms of the stress change.
On the other hand, the length of correlation describes the length at which the stress perturbation in cracks affect the surrounding volume [67]. It has been considered similar to the length of the largest crack during the load cycle: ∝ [68,69]. This is related to a second order power law for heterogeneous materials given by [69]: where is the stress required for the macroscopic failure, and and are parameters that describe the stress evolution. Note that Equation (10) corresponds to the best fit from experiments on rock samples under compression stress. By replacing Equations (7) and (10) into Equation (6), the entropy changes of the whole system prior to the macroscopic failure is: where 2 = 5− . Figure 2c shows how the entropy increases by using Equation (11)  and shear stress change dσ/dt (red curve) of the lithosphere prior to and after the main failure. Here, the earthquake time is t = t C = 5 a. u. Note that a. u. means arbitrary units. (c) Entropy increases by using Equation (10). It is shown that there are two main behaviors: the initial slow increase that lasts up to t ∼ 3.5 a. u. (this trend is represented as a black dotted line) and the fast increases between 3.5 a. u. and t C = 5 a.u. (d) Magnitude expected in terms of the stress change.
On the other hand, the length of correlation ξ describes the length at which the stress perturbation in cracks affect the surrounding volume [67]. It has been considered similar to the length of the largest crack during the load cycle: ξ ∝ l max [68,69]. This ξ is related to a second order power law for heterogeneous materials given by [69]: where σ P is the stress required for the macroscopic failure, and k and p are parameters that describe the stress evolution. Note that Equation (10) corresponds to the best fit from experiments on rock samples under compression stress. By replacing Equations (7) and (10) into Equation (6), the entropy changes of the whole system prior to the macroscopic failure is: where C 2 = Λk 5−D . Figure 2c shows how the entropy increases by using Equation (11) and p = 0.64 [69], σ P = 5 (a.u.) and C 2 = 1 a.u. Overall entropy increases prior to the earthquake. Nevertheless, it increases particularly fast from t = 3.5 a.u. up to the main failure at t = 5 a.u. (purple curve in Figure 2c). This means that the generation of cracks prior the macroscopic failure is part of an irreversible process which maximizes the entropy in time (see reference [70] for the relationship between maximum entropy principle and irreversibility). Figure 3 shows the schematic representation of stress evolution, crack generation, entropy increases, and final rupture. For example, Figure 3a shows the onset of the system when no stress changes are applied. That is, no cracks are generated, and entropy is at lowest value (blue circles). An increase in the stress will create a small number of cracks while the entropy slowly increases (Figure 3b). Before the main earthquake, the stress increases even more while the entropy rises faster (blue circles in Figure 3b). This state is characterized by a large number of cracks within the lithosphere. Right before the earthquake, the entropy, and crack numbers are at maximum values while the stress change is at its maximum acceleration (Figure 3d). time (see reference [70] for the relationship between maximum entropy principle and irreversibility). Figure 3 shows the schematic representation of stress evolution, crack generation, entropy increases, and final rupture. For example, Figure 3a shows the onset of the system when no stress changes are applied. That is, no cracks are generated, and entropy is at lowest value (blue circles). An increase in the stress will create a small number of cracks while the entropy slowly increases (Figure 3b). Before the main earthquake, the stress increases even more while the entropy rises faster (blue circles in Figure 3b). This state is characterized by a large number of cracks within the lithosphere. Right before the earthquake, the entropy, and crack numbers are at maximum values while the stress change is at its maximum acceleration (Figure 3d). Finally, (d) shows that the maximum entropy change is found right before the impending main failure (marked as green area).

Seismic Moment and Entropy
As the entropy is rapidly increasing before the main earthquake, let us consider the effect of Equations (8) and (11) into other seismic parameters as the earthquake's magnitude. The ruptured area for earthquakes can be approximate by = 2 [17]. Then the area in terms of the entropy from Equation (8)   Finally, (d) shows that the maximum entropy change is found right before the impending main failure (marked as green area).

Seismic Moment and Entropy
As the entropy is rapidly increasing before the main earthquake, let us consider the effect of Equations (8) and (11) into other seismic parameters as the earthquake's magnitude. The ruptured area A for earthquakes can be approximate by A = πl 2 max [17]. Then the area in terms of the entropy from Equation (8) is: On the other hand, the experimental relation between seismic moment M 0 and the rupture area is given by A = a A M 2/3 0 , where a A = 1.34 × 10 −10 m 5 s kg 2 1/3 [71]. By introducing Equation (12), the seismic moment is: where Γ = a −3/2 A Γ 3/2 0 . Note that Equation (8) states that dS ≥ 0. Here, the seismic moment is only possible if dS = 0 and dσ/dt = 0. That is: dS > 0. In other words, earthquake's occurrence requires a change of stress which generates an irreversible process that increases entropy. The Equation (13) into the seismic magnitude equation [1] gives: Equation (14) corresponds to the most probable expected magnitude of impending earthquake by regarding the entropy of the lithosphere. That is, if entropy change or stress change are known, an estimation of the expected magnitude at that moment could be obtained by using Equation (14). This can be seen in Figure 2d for the sigmoidal stress change. Note that if no earthquake occurs, it is implied that entropy or stress will keep increasing. This means that the expected magnitude will also increase.

Entropy and Fractal Geometry Generalization for Linear Nonequilibrium Thermodynamics
Let us consider now the generalization for Equations (8) and (9). The first step is to consider the entropy change in time which is defined as [72]: where X k and J k are the N thermodynamics forces and flows in the system. These forces and flows are related by the phenomenological equations which are defined as: where K kl and L kl are the phenomenological coefficient. Specifically, the resistance and conductance coefficients respectively. The phenomenological coefficient obeys the Onsager's relations when no external magnetic field or Coriolis force are present [47,48]: The condition from Equation (17) is satisfied because the magnetic anomalies are the results of the crack generation. Magnetism is not generating or affecting the lithospheric stress states as shown in Equation (3). Then, Equation (17) for different indices becomes: Additionally, the matrix L kl and K kl are related by: Let us consider the entropy change dS µ /dt for different subdomain µ (µ = 1, 2, . . . , m). Then, Equation (15) becomes: and the thermodynamic forces and flows are now expressed as: Replacing the forces from Equation (21) into (20) gives: Here, the phenomenological coefficient cannot be negative while flows (or forces) are quadratic in forms which implies that dS µ /dt ≥ 0 [73]. Equations (16) and (21) show that all the forces X i can be generated by all the flows J i (and vice versa) by a linear combination of the resistance (or conductance) coefficients. For example, Equation (4) shows that the proportionality between electrical currents and volts holds when the distance at which the volt is considered is not zero. In other words, the phenomenological coefficients are related to the domain's length for this case. Similarly, thermal conductivity depends on the Knudsen number which depends on the characteristic length [74][75][76] while the hydraulic conductivity is proportional to the hydraulic permeability which is a measure of the pore geometry of the pore structure [77]. In addition, the resistivity K can be related to the local metric tensor which represents the geometric measurement of the distance traveled [78].
Other examples of the relation among thermodynamics forces, flows, and phenomenological coefficients that require geometrical properties as volume or characteristic length can be seen in Table 3.1 in reference [79]. The above-mentioned examples suggest that phenomenological coefficients K µ kl describe the geometrical properties as the characteristic length r µ kl which defines those domains where X µ k and J µ k are valid. This means that the phenomenological coefficients can be written as: where η 0 is dimensionless constants, r 0 a constant of units of length, K µ kl correspond to a constant with the phenomenological coefficient units, and α is a constant that determine the dimension of the characteristic length r. That is, α represents Euclidian dimensions or fractal dimensions. Then, the total entropy change is: In other words, the total entropy depends on the entropy of each subdomain (dS/dt = dS/dt(dS µ /dt), (µ = 1, 2, . . . , m)). On the first hand, a general self-affine trans- And ξ 11 , ξ 22 , . . . , ξ µµ are constants that represent dilatancy (ξ µµ > 1) or contractions (ξ µµ < 1) among different subdomains M µ [80]. Here, Equation (24)  Mathematically, this scale invariance can be written as [80]: where λ is a factor and α 0 is a constant. To consider our system as self-affine, each length ξ µµ must be constant and characterize the length of each subdomain µ. This means that if the self-affine property is applied to Equation (23), it implies that r µ kl does not depend on indices k and l. In other words, the N net forces and flows of each subdomain µ are restricted within a specific length scale (r µ ) α . Thus, the latter turn Equation (23) into: whereX k ,Ĵ k , andŜ µ are defined as the thermodynamics forces density, thermodynamics flows density, and entropy change density and include the quadratic forces. Then, it is possible to observe from 26 that f (λ, r) = dS dt (λ, r) = λ α dS dt (r) = λ α f (r) which satisfy the scale invariance definition from Equation (25) and where α 0 = α. In the continuum self-affine limit, the total entropy is: On the second hand, a general fractal geometrical volume V E is obtained when V E ∼ r D E dN, where N = k 1 r −D , k 1 is a constant D E is the Euclidian dimension, and D is the fractal dimension [66]. This gives the relation: V E ∼ −Dk 1 r D E −D−1 dr. Then, the entropy is fractally distributed when α = D E − D − 1. In other words, Equation (27) is a generalization of the relation dS/dt = V dS v /dtdV used in reference [73], where S v is the entropy density (chapter 3.9.1). If the entropy change density is fractally distributed, the contributions are only a dilated or contracted version of a constant dŜ 0 /dt (Equation (26). Thus, the contribution of entropy change of Equation (27) is independent of r. Thus, the total entropy change is: The constant k 1 can be found by using the topological dimension D T that obeys A T = −Dk 1 r D T −D−1 dr, where A T is a known positive constant. Then, the total entropy change is: As Equation (22) is defined always as nonnegative, Equations (29) and (30) can be written in a more general manner as: . Here a fixed r min is used. Note that Equation (6) is recovered if we consider the volumetric contribution V by using the Euclidean dimension D E = 3, D T = D E − 1, η 0 = 4π/3 and A T = πr 2 max . While the quadratic forces are represented by dŜ 0 /dt as shown in Equation (26).
Finally, it is possible to observe that the volume part V in Equation (31) grows when the larger fractal structure length r max = r max (t) is larger. In other words the larger the fractal volume, the larger the entropy increases.
Multifractal Entropy for Linear Nonequilibrium Thermodynamics Equation (31) depends on the fractal dimension D. Nevertheless, D is not a fixed value because there could be infinite D values between D T and D E . This raises the possibility to take into account a system composed by subdomains characterized by different fractal dimensions. If each subdomain is restricted to a specific volume, it is possible to consider the volume growth of several non-interacting fractal structures by adding k different entropy contribution. That is: Equivalently: Then, the contribution of several non-interacting fractals volumes' also increase the total entropy of the system. The interacting systems case is different.

Discussion
Cracks and fractures within the lithosphere are well described by scaling laws or fractals distributions [81]. This implies that experiments on rock samples could also give information regarding the geodynamic scale. One relevant property of rock samples is the increasing number of cracks before the main failure [7,8,54,69]. This means that the role of cracking and its electromagnetic signals correspond to a pre-failure feature that can be used as a forecast for major earthquakes [4]. On the other hand, it is known that the lithospherical surface deformation is a feature of the interseismic cycle which is observed in the middle of two large earthquakes [82]. Then, the link between interseismic and the pre-failure process could be stated as it follows: the surface deformation is the first reaction to the stress increases while the cracking generation rises when the lithosphere cannot hold more strain. That is, the cracking is generally generated when no deformation is clearly observed. Up to some point, the strain and cracking will not be sufficient to hold the still-increasing stress. It is at this moment when the main earthquake is expected to occur.
Here, as Equation (8) shows that the crack generation rapidly increases the entropy, it could be considered that the fractally distributed cracking process corresponds to the main manifestation of lithospheric dynamics despite the not clearly observed deformation. Note that large scale strain-rate perturbation still could exist [83]. As the main earthquake is also a manifestation of entropy increases, the entropy change can be used in order to link the prefailure dynamic to the main failure (or crack) in order to estimate an approximate magnitude of the impending earthquake. This can be done by considering the Equations (13) and (14). These Equations show that the larger the entropy (or stress) change, the larger the expected earthquake magnitude. In addition, Equations (13) and (14) also state that no earthquake can be expected if there is no entropy (stress) change. In other words, both the main rupture and crack generation are linked by the entropy increases. This means that the earthquake process must be preceded by the crack generation within the lithosphere that also could generate electromagnetic signals. This contradicts those claims that consider the space weather and solar activity as the main source of pre-failure electromagnetic signals [84]. Despite this, there is one manner in which space weather could influence the lithospherical cracking. That is, by means of external forces. Specifically, part of the stress changes might be generated by solar activity. This also implies that a reliable space weather mechanism must be presented. Otherwise, if no relation between solar activity and lithospherical stress changes is confirmed, implies that the cracking process required to increase the lithospherical entropy is driven only by tectonic forces.
As those tectonic forces increase the number of cracks, it is implied that they also increase the interaction between fractures and pores which leads to the increase of the fracture permeability [85]. In addition, the Darcy's law states that the larger the permeability, the larger the fluid flow within the media [86]. This means that the increase of fractal cracks imply the increase of the permeability, which also allows fluid migration within the lithosphere. This fluid migration disturbs the lithospheric effective stresses [87] which could destabilize stress states within faults [88]. Additionally, the fluid migration could also flow upward the Earth's surface while carrying high temperatures, different gases, and electric charges that could ionize the lower atmosphere [89][90][91][92]. The latter physical description can be added to the cascade of physical processes that rise due the entropy increases. The whole schematic representation can be seen in Figure 4. Here, the main physical property that triggers the cascade of other physical phenomena is the stress increase and the increase of the entropy. The green arrow represents the different process that can be explained by the seismo-electromagnetic theory when non-equilibrium thermodynamics is considered. This includes the fractal cracking, electromagnetic signals, fluids migration, frictional changes, b-value changes, main earthquake generation, and other less direct processes such as gases liberation [93], electrical charge movement [94], or ionospheric anomalies (see Figure 10 in reference [89] for the cascade of physical process considered as not direct effect). The black arrows show the classical seismological relations such as tsunamis, aftershocks, and gravitational signals. Figure 4 shows that the seismoelectromagnetic theory in the context of non-equilibrium thermodynamics complements classical seismology and gives a multidisciplinary physical explanation of the earthquake's generation. Note that this description is regardless of the geological context because all of these cascades of phenomena have been reported in different seismic events and rises because of the cracking of intact rock [95][96][97].
Regarding the geometrical distribution of cracks, this is a key feature that allows us to obtain Equations (8), (13), and (14). This is the same fractal distribution which is also observed in other fields under other names as "Constructal law" [44,98]. That is why Equations (29)-(31) were obtained in Section 2.3 and correspond to the generalization of Equation (8). Specifically, it is shown how geometrical properties of the phenomenological coefficients (Onsager's relations), which relates thermodynamics forces and flows, can also be valid for fractal distribution (Equations (29)-(31) in Section 2.3). Note that Equations (31) and (33)   It is possible to observe that the cascade of physical phenomena start with stress changes that increase the entropy of the lithosphere. Then, the green arrows represent the seismo-electromagnetic branch that explain the observed seismic and non-seismic pre-earthquake measurements. The black arrows represent the classical seismic domain. Note that the physics that explain earthquake occurrences come from the seismo-electromagnetic domain. That is, the change of b-value, the main earthquake, and secondary effects and fault's frictional changes are due the entropy increases. Those green-black arrows represent relation that can be stablished by classical seismology and seismo-electromagnetic phenomena.
It is also important to note the link between the Constructal law, the Onsager's relations, the metric tensor, and multiscale thermodynamics. As the Constructal Law describes the energy flux of natural systems that are characterized by fractal geometry [99], and the Onsager's coefficients describe the thermodynamic flux in non-equilibriums systems (Equation 16), it is possible to state that the Constructal law is equivalent to the phenomenological coefficients when time-dependent fractal geometry is considered. Furthermore, as the metric tensor is related to the Onsager's resistivity [78], and the fractal volume in Equations (31) and (33) scales (Equation (26)), which implies that the metric tensor could have fractal or multiscale properties. This also implies that future works should relate this fractal entropy (Equation (32)) to the multiscale thermodynamics [100] or even cosmological evolution [101] and quantum gravity [102,103]. This is because these kinds of links would allow other deepest question to be asked, such as: are fractals the results of random fluxes and forces or the results of space time properties? Are those forces and fluxes being shaped by the fractal metric tensor? Are fractals the milestone required to link quantum and gravity realms? More work must be done in this direction. It is possible to observe that the cascade of physical phenomena start with stress changes that increase the entropy of the lithosphere. Then, the green arrows represent the seismo-electromagnetic branch that explain the observed seismic and non-seismic pre-earthquake measurements. The black arrows represent the classical seismic domain. Note that the physics that explain earthquake occurrences come from the seismo-electromagnetic domain. That is, the change of b-value, the main earthquake, and secondary effects and fault's frictional changes are due the entropy increases. Those green-black arrows represent relation that can be stablished by classical seismology and seismo-electromagnetic phenomena.
Equations (29)-(31) also show the general relation between the growth of a fractal distribution and the increase of entropy in time. By considering this, it is possible to claim that the growth of fractals in nature correspond to the rise of one kind of irreversible 'disorder' governed by Equations (29)- (31).
It is also important to note the link between the Constructal law, the Onsager's relations, the metric tensor, and multiscale thermodynamics. As the Constructal Law describes the energy flux of natural systems that are characterized by fractal geometry [99], and the Onsager's coefficients describe the thermodynamic flux in non-equilibriums systems (Equation 16), it is possible to state that the Constructal law is equivalent to the phenomenological coefficients when time-dependent fractal geometry is considered. Furthermore, as the metric tensor is related to the Onsager's resistivity K [78], and the fractal volume in Equations (31) and (33) scales K (Equation (26)), which implies that the metric tensor could have fractal or multiscale properties. This also implies that future works should relate this fractal entropy (Equation (32)) to the multiscale thermodynamics [100] or even cosmological evolution [101] and quantum gravity [102,103]. This is because these kinds of links would allow other deepest question to be asked, such as: are fractals the results of random fluxes and forces or the results of space time properties? Are those forces and fluxes being shaped by the fractal metric tensor? Are fractals the milestone required to link quantum and gravity realms? More work must be done in this direction.

Conclusions
This work has described the thermodynamics of fractals cracks presented in the seismoelectromagnetic theory prior to main earthquakes. This example has been useful in order to generalize the thermodynamics of systems characterized by fractal geometries. Regarding the fractal cracking process, it is possible to conclude that:

•
As Equation (8) is always positive, it is implied that the generation of cracks are the manifestation of irreversible process.

•
The pre-failure and failure process can be linked by means of the entropy changes.

•
The seismic moment and magnitude exist if external stress, that increases of the entropy of the lithosphere, and increases in the number of cracks and electromagnetic signals also exist. • It is possible to estimate an expected seismic magnitude in terms of the entropy change/stress change. • Entropy rapidly increases before earthquakes.

•
No entropy increase, no earthquake.

•
The seismo-electromagnetic theory explains the non-seismic pre-earthquakes signals and gives physical foundations to the generation of earthquakes.
Regarding the generalization of other non-equilibrium system characterized by fractal properties, it can be concluded:

•
The tendency in which nature creates fractals corresponds to a geometrical manifestation of that tendency in which the universe increases the entropy. • Fractals rising in several fields and topics reveals the increase of 'disorder' of those systems.

•
The phenomenological coefficients can describe geometrical properties of forces and fluxes.

•
The Constructal law is one geometrical application of Onsager's relations.

•
The entropy density is defined as dŝ, which represents the quadratic time derivative of those forces ( . F 2 ) that generate the fractal geometry V f r . No changing force F implies no fractality. • More work must be done in order to link metric tensor, fractal entropy and multiscale thermodynamics.
Author Contributions: P.V.-A. proposed the core idea, mathematical development, and initial draft of the project. E.G.C. and D.L. contributed to the scientific discussions of the work. All authors have read and agreed to the published version of the manuscript.
Institutional Review Board Statement: Not applicable.

Informed Consent Statement: Not applicable.
Data Availability Statement: Not applicable.