The total energy of geomagnetic ﬁeld

Geomagneticﬁeldof Earthhas been existed during the lastthree billionyears.The geomagnetic ﬁeld isaformoflosingthegeothermalenergyofEarth.Inthis paper, the formula of total energy of geomagnetic ﬁeld, which is used in the volume from the surface of Earth to unlimited space, is obtained by the function of scalar potential of IGRF (1900–2010). Generally, the total energy of geomagnetic ﬁeld is decreasing by the linear form in these years. The general power of geomagnetic ﬁeld about surface of Earth was also steady about 949.7 kW. However, a damping oscillation form with about 60 year’s period is appeared in the total energy of geomagnetic ﬁeld. Maybe the damping oscillation form is sourced by effects of ﬂuid motions on the core mantle boundary. Maybe the linear decreasing form in the total energy is the variation of the strong ﬁeld sourced at the center of the Earth’s inner core.


Introduction
Strong field is sourced by the dynamic processes in sub-shell of georeactor at the center of the Earth's inner core (Herndon 2009). The strong field spreads via the inner core, the liquid core, the mantle and the crust as well as it is observed as geomagnetic field.
The chain reaction's energy of nuclear fission of heavy nuclear isotopes (U-238, U-235 and Th-232) in georeactor becomes a main source of geothermal energy of Earth. The geothermal energy is lost by all geological activity phenomena such as volcanic activity, Earthquakes, continental drift and geomagnetic field. It is difficult that the total energy of volcanic activity, Earthquakes and the continental drift is estimated for whole Earth. The geomagnetic field is a form of losing the geothermal energy of Earth. However, the total energy of geomagnetic field solved for whole Earth is easier than others. The geomagnetic field of Earth has been also existed during the last three billion years. Therefore any phenomenon could be observed in the total energy of geomagnetic field.
Thus our purpose is to study the total energy of geomagnetic field for whole Earth. We are interesting the power of geomagnetic field about Earth's surface as well.
in 1839. This method obtains the scalar potential for geomagnetic field form the Laplace's equation. International Geomagnetic Reference Field (IGRF) is the negative spatial gradient of the scalar potential as well. It can be written as where V-is the scalar potential for geomagnetic field solved form the Laplace's equation. And the V-is rewritten in the following form as here a-is reference radius of Earth (a = 6371.2 km). The r, h and u are the geocentric coordinates. r-is the radius in kilometers, h-is the co-latitude (h = 90°-latitude), and u-is the longitude. The coefficients g m n and h m n À Á are Gaussian coefficients put forth by the IAGA for the IGRF, and P m n cos h ð Þ ð Þ represents the Schmidt quasi-normalized associated Legendre functions of degree n and order m.
If the Eq.
(2) is substituted to the Eq. (1), the components of the Eq. (1) could be written as the following forms B r ¼ ÀoV and where B r , B h and B u represent the field strength in local tangential coordinates. P m n cos h ð Þ ð Þis rewritten in the following form as here d 0 m -is Kronecker delta and l equals with cos(h) (l = cos(h)) (Winch et al. 2005).
The Eq. (6) is called the quasi-normalized function of associated Legendre function. Now, the density of energy of geomagnetic field could be solved from the Eqs. (3-5) as here l 0 equals 4p 9 10 -7 H/m (l 0 = 4p 9 10 -7 H/ m). In other word, the density of energy of geomagnetic field is solved how the Eqs. (3-5) are substituted to the Eq. (7). And the total energy of geomagnetic field could be solved from the Eq. (7) as To solve the integral in Eq. (8) The integral on the Eq. (8) could be divided two terms with g m n À Á 2 and h m n À Á 2 as The first term (W 1 ) in Eq. (13) can be rewritten as If the orthogonal properties on the Eqs. (9-12) are used, the expression in the Eq. (14) could be written as And the second term (W 2 ) of the Eq. (13) could be written as Now, the Eq. (13) could be rewritten in the following form The total energy of geomagnetic field is expressed by the Eq. (17). The values of Gaussian coefficients g m n ; h m n À Á in the Eq. (17) are taken from the website (doi: 10.1111/j.1365-246X.2010.04804.x).

Results and discussion
The Eq. (17) could be called ''the formula of total energy of geomagnetic field''. The total energy of geomagnetic field can be solved by the Eq. (17) in the interesting volume from the Earth's surface to unlimited space. The values of Gaussian coefficients in the Eq. (17) are determined by measured values on the Earth's the surface and satellites in orbits about Earth. Therefore we think that the formula in the Eq. (17) fits clearly in the volume from the Earth's surface to the ionosphere. We are also very interesting the power of geomagnetic field about surface of Earth. Thus the total energy of geomagnetic field is solved in the volume from the Earth's surface (a = 6371.2 km) to the ionosphere (a ? 80 km) in this paper. Values of the total energy solved by the Eq. (17) are written on the Table 1.
The dependence on time of total energy of geomagnetic field is depicted on the Fig. 1. There are two types of dependence on time in the Fig. 1. The first type is linear and it is expressed as The power of geomagnetic field could be solved from the Eq. (18) as where P-is power of geomagnetic field about surface of Earth. Now, the unit of the power must be expressed by the unit of SI system as   If the total energy of geomagnetic field is changed clearly by the linear decreasing law as the Eq. (18), it must be equaled zero after 1126 years or in 3026. But geomagnetic field of Earth has been existed during the last three billion years. Thus, the linear decreasing form on the Fig. 1 could be a small part, which is appeared in the linear form, of any oscillation with the very long period such as longer than 600 years.
If the linear dependence on time is subtracted from the Fig. 1, the clear periodical dependence on time in the energy could be differed. Its analytic form could be obtained by the least square method as where t is time or variable, and a 1 = 0.025219512, b = 0.00892494537, a 2 = 0.098482076, a 3 = -1.62740646 and a 4 = -0.04020871 are parameters. It is depicted on the Fig. 2. The coefficient of determination of Eq. (21) is R 2 = 0.9037. The function in Eq. (21) is a damping oscillation. In other word, the amplitude of oscillation on the Fig. 2 is decreasing. The periodical dependence on time in the Fig. 2 could be a part with 110 years of a large damping oscillation.

Conclusions
Generally, the total energy of geomagnetic field is decreasing along to the line expressed in the Eq. (18) as the Fig. 1. The general power of geomagnetic field was also steady about 949.7 kW in the years . Therefore, the variation of the power of geomagnetic field about Earth's surface could be equaled zero in the few years. This condition, which the variation of the power equals zero in the few years, is mathematical main condition to solve the drift of geomagnetic field about Earth's surface.
Two types of dependencies on time are observed in the total energy of geomagnetic field. The first type is the linear decreasing form as the Eq. (18). The second type is the damping oscillation form with about 60 year's period on the Fig. 2.
We think that the linear decreasing form could be the variation of the strong field sourced at the center of the Earth's inner core. Maybe the total energy of geomagnetic field is an opportunity to study the properties of the strong field.
However, the damping oscillation with about 60 year's period could be sourced by the fluid motions on the core mantle boundary (CMB). In other words, the effects of the fluid motions on CMB could be appeared as the damping oscillation form in the total energy of geomagnetic field about Earth's surface.