Nonlinear Whole Seismology, Topological Seismology, Magni-tude-Period Formula of Earthquakes and Their Predictions

: First, we propose the nonlinear whole seismology and its three basic laws. Next, based on the nonlinear equations of fluid dynamics in Earth ’ s crust, we obtain a chaos equation, in which chaos corresponds to the earthquake, and shows complexity on seismology. But, combining the Carlson-Langer model and the Gutenberg-Richter relation, a simplified nonlinear solution and corresponding magnitude-period formula of earthquakes may be derived approximately. Further, we research the topological seismology. From these theories some predictions can be calculated quantitatively and are already tested. Combining the Lorenz nonlinear model, we may discuss the earthquake migration to and fro. Finally, if various modern scientific instruments, different scientific theories and some paranormal ways for earthquake are combined each other, the accuracy of multilevel prediction will be increased.


Introduction
It is recognized that the earthquakes are very complex nonlinear phenomena, and many theories and some phenomenological descriptions have been proposed, such as fractals, and propagation and interaction on the seismic wave in the nonlinear media [1][2][3]. Carlson, Langer, et al. [4][5][6][7][8], presented the Burridge-Knopoff block-and-spring model of an earthquake fault, and discussed basic properties, predictability and so on.
Based on the Gutenberg-Richter (GR) relation, we proposed the magnitude-period formula of the earthquake [9,10]: Topology is the mathematical study on properties of space preserved under continuous deformations including stretching and bending, but not tearing or gluing. It includes connectedness, continuity and boundary. Topology is widely applied to various aspects [11,12]. In this paper, we propose the nonlinear whole seismology, topological seismology and their applications.

Nonlinear Whole Seismology and Its Basic Laws
Based on the most basic features of modern seismology and combined geological structure and seismic mechanism, we propose the nonlinear whole seismology and its three basic laws: Yi-Fang Chang 2 of 10 First law: Seismology as whole is a complex huge system of violent crustal movement. Second law: Various interactions in seismology all are nonlinear. Corresponding equations all are nonlinear, and possess chaos and fractals, etc. Different magnitudes M possess fractals with self-similarity.
Third law: Earth and crust as environment are closely related to earthquakes. They are regarded as boundary conditions of the system, which and the seismological systems have often various nonlinear relations.
Two main characters of seismology are nonlinearity and whole. Three laws are a whole: Because of complexity, the description of seismological systems must apply the nonlinear theory with the interaction terms. Its future development is the nonlinear whole earth science.
We researched possible relations on earthquakes and various cycles of crustal movement, terrestrial magnetic field, Earthrotation, change in climate, and sunspot [9] which is related to nonlinear Earth-Sun physics. They show also whole and complexity.

Nonlinear Geodynamics and Seismology
Based on general geodynamics [13], and combines the Burridge-Knopoff block-andspring model of an earthquake fault presented by Carlson, Langer, et al. [4][5][6][7][8], we proposed that the fundamental nonlinear equations of fluid mechanics in Earth's crust with momentum conservation are: Here the frictional force for Newtonian fluid is: Usual solutions of the nonlinear fluid mechanics are very difficult. In order to simplify the above case, assume that 0 F =0, the pressure gradient gradp=C" is independent of the velocity V, so Eq. (3) turns into This is a nonlinear partial differential equation of the velocity V. Using a similar way for the soliton solution, let This solution is: whose solution is: Therefore, the kinetic energy density is: When r=0 and t=0,   Let the integral constant C>0, the larger the region is, the smaller the energy density  is at the same time t; the larger  be- It is a well-known nonlinear chaos equation. For momentum the control parameter is: This determines the branch-chaos. If our consideration is in more detail,  will be more complex. Various similar results may be approximately derived from the energy conservation equations. This fluid layer may be a general asthenosphere. As an example, the matter density is   4 3 g cm / and   10 22 t is enough long,   1 , so the usual Yi-Fang Chang 4 of 10 crust is stable. But, if the asthenosphere consists of solid, fluid, gas and plasma,  will decrease very quickly. For a special asthenosphere, when   0 75 . , the bifurcation will appear, perhaps it corresponds to the momentum migration between a couple of places.
...., it infers a chaos will appear, and corresponding earthquake occurs. In fact, earthquake will occur as the shorter the series of calm time, which is consistent with a branch-chaos process. This could be a basis of the renormalization schemes [9], which research some seismic activity prior to the main earthquake.
Some factors may cause earthquake, for example, volcanism, reservoir, underground nuclear explosion and so on [15][16][17], which are merely boundary or initial conditions. But, because the seismic system obeys the nonlinear equation, they are extremely sensitive to the initial conditions.
It is known that earthquakes are very complex nonlinear phenomena. At present the nonlinear seismology is an important direction of the development. A main character of nonlinearity is chaos. In Webster's Third Dictionary "chaos" possesses these meanings of chasm, gulf and abyss. These are namely results of earthquake. Turcotte discussed generally fractals and chaos in geology and geophysics [2]. Levin researched nonlinear oscillating structures in the earthquake and seaquake dynamics [3].
The earthquake corresponds to chaos, which shows complexity on seismology, and is useful for the explanation of seismic essence, and casts a fatal shadow in an earthquake exact prediction. However, the possible periodicity in earthquakes may be determined according to the approximate results from a simplified model.

Magnitude-Period Formula of Earthquake
For the forecasts of earthquakes, Kiremidjian, et al. [18] presented a stochastic slippredictable model based on Markov renewal theory for earthquake occurrences. Harris, et al., discussed changes and influence of static stress on earthquake in southern California [19][20][21]. Borodich [14] described some renormalization schemes for earthquake prediction, which can be reduced to a power-law or the log-periodic approximation of the regional seismic-activity data. Marzocchi, et al. [22] provided insights that might contribute to better formally defining the earthquake-forecasting problem, both in setting up and in testing the validity of the forecasting model, and found that the forecasting capability of these algorithms is very likely significantly overestimated. Helmstetter, et al. [23] developed a time-independent forecast for southern California by smoothing the locations of magnitude 2 and larger earthquakes, and using small m 2 earthquakes gives a reasonably good prediction of m 5 earthquakes.
We combine the Carlson-Langer model [4][5][6][7][8], the movable asthenosphere corresponds to the surface, and the faulting rock layer corresponds to the block. From this the Gutenberg-Richter relation and the same results may be obtained, while in such case the masses of blocks are usually different, and the velocities may be changeable, especially at the fissure. Further, from some formulas of the Carlson-Langer model, the magnitudeperiod formula [24] may be derived. In the model, the magnitude M is be defined to be the natural logarithm of the earthquake moment H, Let the moment between two times T 0 and T T 0 + . In the model, the frictional force is which is also a result in a simple mechanical model for earthquake dynamics [25]. In present paper the frictional force is If both forces F are equal, Only the time variable is related, so    − + ut c , which is replaced to Eq.
It is namely my formula (2). Its general formula is (1). In this case, b is a tectonic parameter in the Gutenberg-Richter relation. The fractal characteristics of earthquakes show the fractal dimension D=2b [26]. If k=1, b=0.43429 agrees with the earthquake belt of the longer period of earthquakes, whose b=0.4-0.6. If k=0.5-0.333, b=0.86858-1.3029 [27]. The longer the period T of earthquakes is, the less the seismic number N is. The formula (23) and the Tsubokawa-Whitcomb-Scholz relation [28,29] log , T aM b = − (24) between duration of crustal movement and magnitude of earthquake, have the same form and but different meanings and signs of parameters. The solutions of the one dimensional partial differential equations for Carlson-Langer model have the form [21]: which implies a period. Almost all formulas of the forecasts on earthquakes are based on the Gutenberg-Richter relation, which can estimate the average rate of earthquakes. We assume that T and N are inversely proportional in a first-order approximation, so the magnitude-period formula (2) of earthquakes may be derived from the Gutenberg-Richter relation [9]. But, the formula is built on the nonlinear dynamical theory, so it may be researched widely.

Seismological Topology
Based on the chaos equation (13) of earthquakes, they may have periods two, four, eight and so on, for example, period is two 0.5130 and 0.7995, or four 0.38, 0.83, 0.88 and 0.50, etc. This type of earthquakes may use the Fourier series: It is the superposition of various periods, and whose formula It may replace into the magnitude-period formula (2). According to the Landau equation of the turbulence: Its two simple solutions are: (1). A=0(r=0) is zero-earthquake; (2). The limit cycle of the radius l r / 2 = is the periodic earthquakes. An earthquake field may be proposed phenomenally. Generally, the nonlinear differential equation has the stable multi-solutions, so that the nonlinear earthquakes have the multi-periods. Probably, many periods reach finally chaos, which obtains the strange attractor and determines the central point of earthquake, earthquake region and fractal. And b in Eq.(2) and a self-similar period may be evaluated from the fractal dimension D. The three dimensional qualitative analysis theory of nonlinear equation may determine the longitude-latitude-depth, or the longitude-latitude-time. While a stable singularity shows that this point will be necessarily scene of earthquake. Attractor may be applied to a formation of earthquake.
The continuous fluid may form a topological space, and it is a flow topology. The dynamical system of chaos may be defined by topology [12], such a chaotic meaning of earthquake contacts naturally the topology of earthquake. Further, we research the topological seismology, in which 1). Earthquake is topology torn by tension. 2). The fixed point is an earthquake region. 3). For 2-, 3-, 4-dimensional topology, earthquake point corresponds to a cut point (a sensitive point), and may develop to a small region, which contacts network, for which the self-similarity and fractal exist possibly.
which is a topological invariant. For the 2-and 3-dimensional spaces, the simplest cases are  =1 and 2. If V(t) is a function of time, E(t) and F(t) will be relative to time t, but  should be invariant. 5). A curved surface with k-holes may be projected to a plane with hole (≤k).
The equations of fluid dynamics on geodynamics combine the phase transition point  and Lorenz model, etc., we can research the identification scales on stability of the earth's crust and on occurred earthquake. These may introduce the Lyapunov exponent.
The topological seismology shows seismology as a whole. Its development is the topological earth science [30].
Moreover, we propose that research of earthquake may be: (1). The space-time point of earthquake is described by computer, and further computer simulation and programs of earthquakes. (2). Usual earthquake is stochastic, and corresponds to chaos. But, there is the self-similarity (for magnitude, etc.) and fractal. (3). Carlson-Langer model [5][6][7][8] and many nonlinear equations derive various soliton solutions, which may describe the earthquake wave. (4). When b in Eq.(2) develops to a function of time b(t), the period T is related to time t, and N=C/T (C is a function). (5). The sink and source by a qualitative analysis theory of nonlinear earthquake equations should be the center earthquake. (6). For a region, we may describe the magnitude-time, the space-time-magnitude, and the general 4dimensional longitude-latitude-time-magnitude figure.

Predictions of Earthquakes
According to the magnitude-period formula (2), so long as one supposes a possible period of earthquakes, for example, the period 34 year in California, other periods can be calculated approximately. Bi and Yuan proposed a periodic scale 250 year for change of weather in China from various aspects. We extended the periodic scale T 0 =250 year to a period of earthquakes for M 0 =7, and results are consistent with many known earthquakes in China [9]. From this the quantitative calculations can be given. The period T that corresponds to the magnitude M will depend on b only, if the scale is determined. Either evaluation of T from b or evaluation of b from T agrees approximately with the facts. It is a simplified and calculable model.
On earthquakes in California Smith, et al., searched strain adjustments associated with earthquakes in southern California [31]. Kovach researched source mechanisms for Wilmington oil field, California, subsidence earthquakes [32]. Mavko, et al., discussed effects of the 1983 Coalinga, California, earthquake on creep along the San Andreas fault [33]. Oppenheimer, et al., obtained the fault plane solutions for the 1984 Morgan Hill, California, earthquake sequence as evidence for the state of stress on the Calaveras fault [34]. Reasenberg, et al., discussed earthquake hazard after a main shock in California [35]. Jaumé, et al., researched change in the state of stress on the southern San Andreas fault resulting from the California earthquake sequence of April to June 1992 [36]. Stein, et al., searched stress triggering of the 1994 M = 6.7 Northridge, California, earthquake by its predecessors [37]. Dodge, et al., observed detail on California foreshock sequences as implications for the earthquake initiation process [38]. Deng and Sykes studied evolution of the stress field in southern California and triggering of moderate-size earthquakes for a 200-year perspective [39], and discussed stress evolution in southern California and triggering of moderate-, small-, and micro-size earthquakes [40]. Rhoades applied the EEPAS model to forecasting earthquakes of moderate magnitude in Southern California [41]. Based on the basic principle of general model of forecasts, Gerstenberger, et al., discussed real-time forecasts of tomorrow's earthquakes in California [42]. Ebel, et al., searched the non-Poissonian earthquake clustering and the hidden Markov model as bases for earthquake forecasting in California [43].