Impedances for induction soundings of the Earth’s mantle

Determination of impedances is necessary in order to eliminate some shortcoming of our knowledge about structures of the exciting source fields and their fickleness. The experimental impedances for induction soundings result from the impedance boundary conditions or heuristic models. The simplified models give just a rough idea of their domain of applicability. Impedances can depend on many factors, including the exciting field structures of several source types which are present in the period range of the mantle soundings (10⁴−4×10⁸ s). The problem in the mantle investigations arises if impedances measured by different methods have to be jointly inverted in order to essentially prolongate the analyzed period range and hence to increase the reliability and depth of induction soundings on land. The subject of our work is an analysis of the known magnetotelluric and magnetovariation impedances to suggest a physically substantiated approach for their joint inversions.     

Seismic Studies of the Earth’s Core

Izvestiya, Physics of the Solid Earth - Abstract—The paper presents the review of the conceptually most important results of seismological studies of the Earth’s core and their...     

Models of the Earth’s plasmapause position

正The plasmasphere is a region of relatively dense(~10–10000 cm~(–3))plasma,surrounding the Earth and extending to distances of about five Earth radii(R_E).It is filled with large amount of cold(~1 e V)plasma originated from the Earth’s ionosphere and co-rotating with the Earth due to the large scale co-rotation electric field.The outermost     

Using the Digital Models of Gravity Anomalies for Zoning of the Earth’s Gravity Field

This paper addresses the questions associated with using the models in regioning the Earth’s gravity field depending on its roughness and type of anomalies. A formalized technique for zoning based on the analysis of the intensity and variability of gravity anomalies is suggested. The suggested technique includes preprocessing the initial gravity anomalies aimed at eliminating the noise component; partitioning the studied region into the elementary segments; calculating the primary characteristics of the complexity and roughness (the intensity and variability of gravity anomalies) for each segment; and, based on these characteristics, classifying each segment into a particular category of complexity. The proposed system of classification of segments relies on the use of three classes of intensity and three classes of variability of gravity anomalies and four categories of complexity of the regions, which are derived from these classes. As a result of applying the technique, the mapped territory is subdivided into the regions that are uniform in terms of their geophysical properties. The developed technique is used for comparing the quality of the different digital models of gravity anomalies, including the Russian RGM model (A.P. Karpinsky Russian Geological Research Institute (VSEGEI)) and the WGM model (Bureau Gravimétrique International (BGI)) in the region of the Mongol–Okhotsk orogenic belt and Sikhote–Alin fold system. The results of the study can be used for zoning of the other geophysical fields and for planning the locations of the new survey networks for increasing the accuracy of the initial data used for compiling the maps.     

Use of Green’s function for determining the disturbing potential of an ellipsoidal Earth

A spherical approximation makes the basis for a majority of formulas in physical geodesy. However, the present-day accuracy in determining the disturbing potential requires an ellipsoidal approximation. The paper deals with constructing Green’s function for an ellipsoidal Earth by an ellipsoidal harmonic expansion and using it for determining the disturbing potential. From the result obtained the part that corresponds to the spherical approximation has been extracted. Green’s function is known to depend just on the geometry of the surface where boundary values are given. Thus, it can be calculated irrespective of the gravity data completeness. No changes of gravity data have an effect on Green’s function and they can be easily taken into account if the function has already been constructed. Such a method, therefore, can be useful in determining the disturbing potential of an ellipsoidal Earth.     

Analysis of Methods for Estimating the Energy of Sources of Acoustic-Gravity Waves in the Earth’s Atmosphere

Izvestiya, Physics of the Solid Earth - Abstract—A comparative analysis is presented of the approaches and methods for estimating the energy of the sources of acoustic-gravity waves (AGW) in...     

Spectral expressions for modelling the gravitational field of the Earth’s crust density structure

We derive expressions for computing the gravitational field (potential and its radial derivative) generated by an arbitrary homogeneous or laterally varying density contrast layer with a variable depth and thickness based on methods for a spherical harmonic analysis and synthesis of gravity field. The newly derived expressions are utilised in the gravimetric forward modelling of major known density structures within the Earth’s crust (excluding the ocean density contrast) beneath the geoid surface. The gravitational field quantities due to the sediments and crust components density contrasts, shown in numerical examples, are computed using the 2 × 2 arc-deg discrete data from the global crustal model CRUST2.0. These density contrasts are defined relative to the adopted value of the reference crustal density of 2670 kgm⁻³. All computations are realised globally on a 1 × 1 arc-deg geographical grid at the Earth’s surface. The maxima of the gravitational signal due to the sediments density contrast are mainly along continental shelf regions with the largest sedimentary deposits. The corresponding maxima due to the consolidated crust components density contrast are over areas of the largest continental crustal thickness with variable geological structure.     

Character of turbulence in the boundary regions of the Earth’s magnetosphere

The statistical features of the magnetic field and ion flux fluctuations in the boundary regions of the Earth’s magnetosphere have been studied on different timescales based on the Interball satellite measurements. Changes in the form and parameters of the probability density function have been studied for the periods when the satellite was in the solar wind plasma, different magnetosheath regions, and the turbulent boundary layer (TBL) at the polar cusp outer boundary. Variations in the probability density function maximum (P ₀) and the kurtosis value as characteristics of the turbulence property evolution on different timescales have been studied. Two asymptotic regimes of P ₀, which are characterized by different power laws, have been found. The structural functions of different orders and the types of diffusion processes in different regions, depending on time variations in the generalized diffusion coefficient, have been studied in order to analyze the character of diffusion processes. For the magnetosheath regions, TBL, and polar cusp, it has been found that the diffusion coefficient increases in the course of time (i.e., the regime of superdiffusion has been obtained). In the foreshock region before the main shock, turbulent processes are described by the Kolmogorov model of classical diffusion.     

Model of muon generation in the Earth’s atmosphere

The muon fluxes on the Earth’s surface and at depths of 7, 20, and 40 m of water equivalent are calculated based on a simple model of pion generation by primary particles with different energies. This generation model is based on the known concepts of multiple pion production. The model parameters are compared with the data obtained using accelerating machines.     

High Resolution Constituents of the Earth’s Gravitational Field

During the General Assembly of the European Geosciences Union in April 2008, the new Earth Gravitational Model 2008 (EGM08) was released with fully-normalized coefficients in the spherical harmonic expansion of the Earth’s gravitational potential complete to degree and order 2159 (for selected degrees up to 2190). EGM08 was derived through combination of a satellite-based geopotential model and 5 arcmin mean ground gravity data. Spherical harmonic coefficients of the global height function, that describes the surface of the solid Earth with the same angular resolution as EGM08, became available at the same time. This global topographical model can be used for estimation of selected constituents of EGM08, namely the gravitational potentials of the Earth’s atmosphere, ocean water (fluid masses below the geoid) and topographical masses (solid masses above the geoid), which can be evaluated numerically through spherical harmonic expansions. The spectral properties of the respective potential coefficients are studied in terms of power spectra and their relation to the EGM08 potential coefficients is analyzed by using correlation coefficients. The power spectra of the topographical and sea water potentials exceed the power of the EGM08 potential over substantial parts of the considered spectrum indicating large effects of global isostasy. The correlation analysis reveals significant correlations of all three potentials with the EGM08 potential. The potential constituents (namely their functionals such as directional derivatives) can be used for a step known in geodesy and geophysics as the gravity field reduction or stripping. Removing from EGM08 known constituents will help to analyze the internal structure of the Earth (geophysics) as well as to derive the Earth’s gravitational field harmonic outside the geoid (geodesy).     

Ranges of AGW propagation in the Earth’s atmosphere

The propagation of atmospheric gravity waves (AGWs) is studied in the context of geometrical optics in the nonisothermal, viscous, and thermal-conductive atmosphere of Earth in the presence of wind shifts. Parametric diagrams are plotted, determining the regions of allowed frequencies and horizontal phase velocities of AGWs depending on the altitude. It is shown that a part of the spectrum of AGWs propagates in stationary air in an altitude range from the Earth’s surface through the ionospheric F1 layer. AGW from nearearth sources attenuate below 250 km, while waves generated at altitudes of about 300 km and higher do not reach the Earth’s surface because of the inner reflection from the thermosphere base. The pattern changes under strong thermospheric winds. AGW dissipation decreases with an adverse wind shift and, hence, a part of the wave spectrum penetrated from the lower atmosphere to the altitudes of F2 layer.     

Impact of Ponderomotive Forces on the Earth’s Magnetosphere

The theory of plasma density redistribution and polar wind acceleration acceleration under the affect on the magnetosphere of the ponderomotive forces induced by the ultra-low frequency electromagnetic waves is presented. Our attention is focused mainly on the important question about the necessity of experimental verification of fairly certain theoretical predictions. It is pointed out that experimental validation is not only necessary for the development of the theory but also for replenishing the knowledge about the structure and dynamics of the near-Earth space. An original method for indirect verification is presented. The idea of this method is based on the dependence of the foreshock locations on the orientation of the field lines of the interplanetary magnetic field (IMF) in front of the magnetosphere relative to the plane of the geomagnetic equator.     

Motion of the Earth’s magnetic poles in the last decade

Based on vector magnetic data from the CHAMP German satellite, average daily spherical harmonic models of the main geomagnetic field to n = m = 10 have been constructed for the period from May 2001 to the end of 2009 at an interval of 4 days. The obtained 16 models, which were averaged over half a year, have been used to calculate the coordinates of the north and south magnetic poles (the points where magnetic field lines are vertical). The changes in these coordinates during these eight and a half years have been traced. Both poles continue moving northward and westward. The north magnetic pole has traveled 400 km during this period. The velocity of its motion has increased up to the year 2003, reaching 62.5 km yr⁻¹, and then started decreasing and reached 45 km yr⁻¹ by the end of 2009. In addition, the direction of motion changed from north-northwestward to northwestward; i.e., the pole started turning slightly towards Canada. The south magnetic pole moved slower by an order of magnitude and has traveled 42 km during this period. The coordinates of the geomagnetic (dipole) poles and the eccentric dipole parameters have also been calculated. The dynamics of these poles has been traced.     

Comparative Assessment of Global Models of the Earth’s Gravity Field

Izvestiya, Physics of the Solid Earth - Abstract—Empirical comparative study of the modern global models of the Earth’s gravity field (EGF) in the form of geopotential spherical...     

Influence of Total Solar Irradiance on the Earth’s Climate

Geomagnetism and Aeronomy - The results of analysis of variations in the total solar irradiance in the 17–24th solar activity cycles and their relation to the climate global warming are...     

Convection in anelastic models of the Earth’s liquid core

The effects of adiabatic cooling on the convection in the anelastic model of the liquid core of the Earth are considered. It is shown that even minor adiabatic cooling causes significant changes in the pattern of the convection, shifting the peak in the convection intensity to the inner part of the core. Just as in the Boussinesq model, both direct and inverse kinetic energy cascades are simultaneously present, and the direct cascade of entropy is observed.     

Seismic response to electromagnetic sounding of the Earth’s lithosphere

The seismic catalogues of 1967?C2008 for the Bishkek geodynamical test site are analyzed for the purpose of studying the response of seismic activity to the electromagnetic sounding of the Earth??s crust during two series of field experiments with high-power controlled sources. The first series of the experiments, which were carried out in 1982?C1990, utilized the pulses provided by a magnetohydrodynamic (MHD) generator. The sounding signals in the second series of the experiments (2000?C2005) were generated by the capacitor-thyristor source ERGU 600-2. In these experiments, temporal variations of the set of statistical parameters characterizing the seismicity, which are typically used in the studies of the background and transient modes of seismicity, were investigated in a selected spatial domain within 150 km from the current electrodes. In terms of time, the analysis was conducted on two levels of detail. The study on a temporal scale of a few years was focused on the variations that preceded, accompanied, and followed the series of the experiments, while the day-scale analysis considered variations that were observed within 10 days after each sounding event. The day-scale analysis yields the following results. The slope of the frequency-magnitude diagram of the earthquakes (b value) during the sounding events is substantially larger than its background value. The slope of the graph gradually becomes gentler within about a day and a half after termination of sounding. The seismic activity slightly enhances during the interval of sounding and abates after its termination to a minimum, which corresponds to the interval of decreasing b value. This character of variations in seismicity differs from the scenario previously established for other transitional seismic regimes. The analysis on a temporal scale of a few years revealed variations in the studied parameters of the seismicity, some of which fall in both sounding intervals of 1983?C1990 and 2000?C2005. However, these variations are not unique; their character and durations suggest their being associated with the processes of preparation and after effects of the strong earthquakes that occurred in the vicinity of the sounding dipole.     

Three distinct types of hotspots in the Earth’s mantle

The origin of mantle hotspots is a controversial topic. Only seven (‘primary’) out of 49 hotspots meet criteria aimed at detecting a very deep origin (three in the Pacific, four in the Indo-Atlantic hemisphere). In each hemisphere these move slowly, whereas there has been up to 50 mm/a motion between the two hemispheres prior to 50 Ma ago. This correlates with latitudinal shifts in the Hawaiian and Reunion hotspots, and with a change in true polar wander. We propose that hotspots may come from distinct mantle boundary layers, and that the primary ones trace shifts in quadrupolar convection in the lower mantle.     

Plume Mode of Thermal Convection in the Earth’s Mantle

In our previous works, based on numerical models, it was shown that under certain conditions a hot material can rise in portions in the tails of thermal mantle plumes. The spectrum of these pulsations can correspond to the observed spectra of catastrophic hotspot eruptions. Since most of the existing numerical models of thermal convection for the mantle of the present Earth do not reveal these pulsations, in this work, we analyze the physical cause and initiation conditions of pulsations of thermal plumes. The results of a numerical solution of the thermal convection equations for a material with varying parameters in the extended Boussinesq approximation are presented. It is shown how the structure of the convection is transformed with the increase of convection intensity. At the Rayleigh numbers Ra > 10⁶, convection becomes unsteady, and the configuration of the ascending and descending flows changes. The new flow emerging at the mantle bottom acquires a mushroom shape with a head and a tail. After the rise of the plume’s head to the surface, the tail remains in the mantle in the form of a quasi-stationary hot steam. It turns out that at Ra ~ 5 × 10⁷, the thermal mantle plume becomes pulsating and its tail is in fact a heated channel through which the hot material rises in successive portions. At the Rayleigh numbers Ra > 5 × 10⁸, the tail of the thermal plume breaks and the plume becomes a regular conveyor of separate ascending portions of the hot material, which are referred to as thermals. Thus, thermal convection with pulsating plumes takes place at the transitional stage from the regime of quasi-stationary plumes to the regime of thermals.