Lunar–solar tide effects in the Earth’s crust and atmosphere

The gravitational interaction in the Earth–Moon–Sun system is considered from the standpoint of influencing the formation of time variations in the geophysical fields and some natural processes. The analysis of the results of instrumental observations revealed the main periodicities and cycles in the time variations of subsoil radon volumetric activity with the same periods as the vertical component of the variations of the tidal force. The amplitude modulation of seismic noise by the lunar-solar tide is demonstrated. It is shown that the intensity of relaxation processes in the Earth’s crust has a near-diurnal periodicity, whereas the spectrum of groundwater level fluctuations includes clearly expressed tidal waves. Based on the data on the tilts of the Earth’s surface, the role of tidal deformation in the formation of the block motions in the Earth’s crust is analyzed. A new approach is suggested for identifying tidal waves in the atmosphere by analyzing micropulsations of the atmospheric pressure with the use of adaptive rejection filters.     

Cyclic variations in the Earth’s flattening and questions of seismotectonics

For more than a decade, the global network of GPS stations whose measurements are part of the International GPS Service (IGS) have been recording cyclic variations in the radius vector of the geodetic ellipsoid with a period of one year and amplitude of ~10 mm. The analysis of the figure of the Earth carried out by us shows that the observed variations in the vertical component of the Earth’s surface displacements can induce small changes in the flattening of the Earth’s figure which are, in turn, caused by the instability of the Earth’s rotation. The variations in the angular velocity and flattening of the Earth change the kinetic energy of the Earth’s rotation. The additional energy is ~10²¹ J. The emerging variations in the flattening of the Earth’s ellipsoid lead to changes in the surface area of the Earth’s figure, cause the development of deformations in rocks, accumulation of damage, activation of seismotectonic processes, and preparation of earthquakes. It is shown that earthquakes can be caused by the instability of the Earth’s rotation which induces pulsations in the shape of the Earth and leads to the development of alternating-sign deformations in the Earth’s solid shell.     

Annual variations in the electromagnetic field of the earth and electrical conductivity of the geological environment

Synchronous annual variations in the geoelectric and geomagnetic field are studied on the basis of long-term electromagnetic monitoring. It is shown that the annual geoelectric variations have intraterrestrial origin and are not related to the annual geomagnetic variations. Temporal variations in the magnetotelluric impedance and magnetic tipper, which characterize the electrical conductivity of the geological environment, are analyzed. It is established that annual variations in the magnetotelluric impedance mainly describe the variations in the electrical conductivity of surface crustal layers and are less sensitive to the deep electrical conductivity of the Earth. The annual variations in the imaginary magnetic tipper at the periods of 1000–3000 s probably reflect the changes in conductivity of a deep transversal low-resistive zone (the fault). It is suggested that annual variations in the geoelectrical and geomagnetic fields, as well as in the electrical conductivity of the geological environment, arise as a response to the changes in the geodynamical processes caused by the revolution of the Earth around the Sun.

     

Manifestation of the Lunar–Solar Tide and Free Oscillations of the Earth in the Variations of the Magnetic Field

The results of processing and analyzing the instrumental observations of the Earth’s magnetic field at the Geophysical Observatory Mikhnevo of the Institute of Geosphere Dynamics of the Russian Academy of Sciences (IGD RAS) for 2010–2015 are presented. Quasi-harmonic components with the periods close to the lunar–solar tidal waves are revealed in the spectra of geomagnetic variations over a period of 0.4 to 30 days. The elliptical S₁ tidal wave which is detected in the geomagnetic variations has modulations with periods of 1/3, 1/2, and 1 year. The spectra of the geomagnetic variations contain peaks corresponding to the free oscillations of the Earth. The analysis of the time series of the magnetic field for the period of the strong earthquakes in the absence of geomagnetic disturbances revealed the fine structure of the Earth’s fundamental spheroidal mode ₀S₂, which splits into five singlets. The established features of the spectrum of geomagnetic variations are helping the development of the new method for studying the deep structure of the Earth and the properties of the inner geospheres for estimating the viscosity of the Earth’s outer core and dynamics of the current systems in the outer (liquid) core, as well as for exploring, with the use of empirical data, the general regularities governing the regimes of energy exchange processes in the geospheres.     

Tidal components of geomagnetic variations

The natural geomagnetic field is constantly disturbed. The total registered effect of geomagnetic variations depends on both planetary and local processes. Planetary sources and sources in the Earth’s core respond to tidal effects. In the accepted model, the complex MHD processes in the Earth’s outer core are approximated by the assumed ring current in the equatorial plain of the liquid core. The geomagnetic variation originating as a result of tidal deformations of ring currents are ～10^?4 and 0.10–1 nT in the liquid core and magnetosphere, respectively. The calculated values coincide in order of magnitude with the processed geomagnetic measurements at Paratunka observatory (Kamchatka region).     

Dynamic processes in the system of interacting geospheres at the Earth’s crust-atmosphere boundary

Results of various instrumental observations of geodynamic processes in the Earth’s crust and geophysical fields at the crust-atmosphere boundary in the Oka area of the Nelidovo-Ryazan tectonic structure and adjacent fractures and in the zone of the Gornyi Altai earthquake of September 27, 2003, are presented. The correlations between the geophysical fields are determined from the results of processing and analysis of microseismic vibrations, the emanation field of natural radon, and variations in the electric field in ground and the magnetic field in the surface atmospheric layer. Tidal deformations of the crust were considered as external effects enhancing the interactions between the geospheres. It is noted that tectonic faults determining the block structure of the crust are recognizable as anomalous variations in the geophysical fields and their high cross-correlation.     

Peculiarities of geomagnetic variations in the central part of the Russian platform

The instrumental observations of geomagnetic variations in the middle latitudes are analyzed. The main periodicities are revealed in the background local variations of the magnetic field of the Earth. Besides the 27–29-day variations, which are associated with the rotation of the Sun around its axis, also the harmonic components with periods of ～1, 6–8, 13–14, 57–60 days, and about a year are identified. It is shown that the periodicities in the background variations are both regular and sporadic. The peculiarities in the time behavior of some spectral components of geomagnetic variations are established.     

Modeling of changes in the tidal response of the environment to an impending earthquake

It has become possible to apply efficiently Beaumont and Berger’s idea (1974) about the possibility to predict earthquakes using continuous monitoring of temporal variations of the tidal response only in recent years, with the appearance of GPS data on the field of the Earth’s surface displacements in seismic regions, as well as with the appearance of the Global Seismographic Network (GSN) containing the tidal tiltmetric and gravimetric data in the vicinities of the epicenters of strong earthquakes in Peru (2001), Sumatra (2004), and Chile (2010) before and after the strongest seismic events. Below, we present the results of the modeling of temporal variations of the tidal response of the environment in the vicinities of the focus of catastrophic earthquakes. The model of the focus includes data on the expansion and orientation of the fault surface, as well as on the discontinuity of the tangential component of the displacement vector on the opposite shores. The model is constructed from the GPS data on horizontal and vertical displacements of the Earth’s surface.     

Study of the recent geodynamic processes in the Kopet-Dag region

The time series of the uniquely long geodetic observations of recent geodynamic processes in the Kopet-Dag region are analyzed. The regional observations of contemporary vertical movements cover a 75-year period; the zonal and local systems of leveling measurements, which provide an increased degree of spatiotemporal detail (the distance between the benchmarks is less than 1 km and the measurements are repeated with a frequency of once per month to two times per annum), have been functioning for 50 years. It is shown that during the last 40–50 years, the regional stress field in the forefront of the Main Kopet-Dag thrust and collision zone of the Turanian and Iranian plates is quasi static. The annual average trend rate of strains estimated from a set of the time series of uniquely long geodetic observations is (3–5) × 10^?8 yr^?1, which is only one to two amplitudes of tidal deformations of the solid Earth. The local deformations in the fault zone reach the rates that are by 1.5–2 orders of magnitude higher than in the block part of the region. It is found that the segments of the Earth’s surface within the axial part of the depression experience persistent uplifting, which indicates that they do not follow the scheme of inherited evolution characteristic of the fault-block structures of the region. It is demonstrated that these anomalous uplifts can be caused by the variations in the weak seismicity in the zone of the North Ashgabat Fault.     

3D models of slow motions in the Earth’s crust and upper mantle in the source zones of seismically active regions and their comparison with highly accurate observational data: I. Main relationships

Constructing detailed models for postseismic and coseismic deformations of the Earth’s surface has become particularly important because of the recently established possibility to continuously monitor the tectonic stresses in the source zones based on the data on the time variations in the tidal tilt amplitudes. Below, a new method is suggested for solving the inverse problem about the coseismic and postseismic deformations in the real non-ideally elastic, radially and horizontally heterogeneous, self-gravitating Earth with a hydrostatic distribution of the initial stresses from the satellite data on the ground surface displacements. The solution of this problem is based on decomposing the parameters determining the geometry of the fault surface and the distribution of the dislocation vector on this surface and elastic modules in the source in the orthogonal bases. The suggested approach includes four steps: 1. Calculating (by the perturbation method) the variations in Green’s function for the radial and tangential ground surface displacements with small 3D variations in the mechanical parameters and geometry of the source area (i.e., calculating the functional derivatives of the three components of Green’s function on the surface from the distributions of the elastic moduli and creep function within the volume of the source area and Burgers’ vector on the surface of the dislocations); 2. Successive orthogonalization of the functional derivatives; 3. Passing from the decompositions of the residuals between the observed and modeled surface displacements in the system of nonorthogonalized functional derivatives to their decomposition in the system of orthogonalized derivatives; finding the corrections to the distributions of the sought parameters from the coefficients of their decompositions in the orthogonalized basis; and 4. Analyzing the ambiguity of the inverse problem solution by constructing the orthogonal complement to the obtained basis. The described approach has the following advantages over the method of steepest descent which was used in our previous works: 1. Application of the perturbation method significantly reduces the volume of the computations in the real problems of coseismic and postseismic deformations (by three to four orders of magnitude when the data from a few dozens of observation points are used); 2. In contrast to the method of steepest descent, the suggested method always provides stable results. This means that adding the new satellite data does not alter the previously calculated coefficients in the low-order harmonics of the distributions of the sought parameters in the orthogonalized basis; this only changes the coefficients of the increasingly higher harmonics which determine the smallscale details in the sought distributions. 3. In contrast to the method of steepest descent, the suggested method is not only capable of constructing stable partial solutions of the inverse problem but also estimating the ambiguity of these solutions. The ambiguity is represented in terms of the superposition of the known functions contained in the orthogonal complement and, hence, with the growth of the amount of the analyzed data it is determined by the linear combination of the increasingly higher harmonics. In the second part of the paper, we present the results of the model numerical computations of Green’s function for the elastic displacements of the ground surface, which correspond to the case of the arbitrary geometry of the dislocation surface and arbitrary orientation of the dislocation vector for the real model of the radially heterogeneous gravitating Earth with the hydrostatic distribution of the initial stresses. The numerical calculations of the creep function in the upper mantle for the coseismic deformations and the ambiguity of the models of postseismic deformations in the vicinity of the source of the Great Tohoku earthquake (Japan) of March 11, 2011 are illustrated by the examples.     

Temporal variations in the tidal response of the medium in the vicinities of the sources of catastrophic earthquakes

The idea of predicting earthquakes by continuously monitoring temporal variations in the tidal response of the medium was suggested by Beaumont and Berger in 1974. However, it became possible to implement the idea only recently. This possibility has arisen due to the deployment of Global Seismic Network (GSN), which collects the data on tidal tilts and gravity in the epicenters of strong earthquakes before and after the strongest events. In this paper, we present the results of model analytical and numerical calculations of the elastic displacements of the Earth’s surface caused by the earthquakes and their preparatory processes. The analytical calculations are limited to the model of a uniform elastic halfspace; the numerical calculations, in addition to this model, also cover the models with radially heterogeneous distributions of elastic moduli in the crust and in the upper mantle, which are determined by the PREM model. We describe the results of modeling temporal variations in the tidal response of the medium in the vicinity of the source of a catastrophic earthquake. The model of seismic source is specified by the length and the orientations of the fault plane and by the value of the discontinuity in the tangential component of the displacement vector on the opposite sides of the fault. The model is based on the GPS data on the horizontal and vertical displacements of the Earth’s surface. We suggest the method for determining temporal changes in the tidal response of the medium in the seismically active regions. This method improves the sensitivity and time resolution of the standard techniques of sliding-window analysis by more than an order of magnitude. The comparative analysis of temporal variations in the tidal response of the medium in the zones of the magnitude 9 earthquake in Japan (March 9, 2011) illustrates the described approach.

     

Temporal variations in the tidal response of the medium in the vicinities of the sources of catastrophic earthquakes

The idea of predicting earthquakes by continuously monitoring temporal variations in the tidal response of the medium was suggested by Beaumont and Berger in 1974. However, it became possible to implement the idea only recently. This possibility has arisen due to the deployment of Global Seismic Network (GSN), which collects the data on tidal tilts and gravity in the epicenters of strong earthquakes before and after the strongest events. In this paper, we present the results of model analytical and numerical calculations of the elastic displacements of the Earth??s surface caused by the earthquakes and their preparatory processes. The analytical calculations are limited to the model of a uniform elastic halfspace; the numerical calculations, in addition to this model, also cover the models with radially heterogeneous distributions of elastic moduli in the crust and in the upper mantle, which are determined by the PREM model. We describe the results of modeling temporal variations in the tidal response of the medium in the vicinity of the source of a catastrophic earthquake. The model of seismic source is specified by the length and the orientations of the fault plane and by the value of the discontinuity in the tangential component of the displacement vector on the opposite sides of the fault. The model is based on the GPS data on the horizontal and vertical displacements of the Earth??s surface. We suggest the method for determining temporal changes in the tidal response of the medium in the seismically active regions. This method improves the sensitivity and time resolution of the standard techniques of sliding-window analysis by more than an order of magnitude. The comparative analysis of temporal variations in the tidal response of the medium in the zones of the magnitude 9 earthquake in Japan (March 9, 2011) illustrates the described approach.     

Tidal and low-frequency currents off the Jaguaribe River estuary (4° S, 37° 4′ W), northeastern Brazil

This article characterizes the spatial and temporal current variations, in the subtidal and tidal ranges, during the rainy and dry seasons, at the continental shelf off the Jaguaribe River, through measurements of continuous current field data from an acoustic Doppler current profiler (ADCP) mooring during 124 days, from June 12 to October 14, 2009. To support this dataset, we collected corresponding data from a meteorological station located at the estuary. The spatial variation showed that highest current speeds occur near the coast, with an offset of a NNW coastal jet, decreasing intensity, monotonically, towards offshore up to 0.1 ms^?1. In the rainy season, small inversions of the wind field were observed, lasting 2 to 3 days on average and were accompanied by the direction of surface currents only. In the dry season, the period of reversal of wind fields and currents lasted 14 and 35 h, respectively. The analysis of empirical orthogonal functions in rainy and dry seasons showed that the continental shelf is predominantly barotropic, where the second and third modes explained only 7% of the total variance, during the dry season. The tidal currents are more intense in the direction normal to the coast, showing a semidiurnal tidal regime. Energy distribution between tidal currents and currents of longer periods showed that for the component parallel to bathymetry, subtidal frequency currents are dominant, contributing to more than 70% of the variance. For the normal component to the coastline bathymetry, there is a significant increase of power concerning tidal currents, at all depths, so they contribute with about 55% of the total variance.     

中国大陆精密重力潮汐改正模型

利用理论和实验重力固体潮模型,充分考虑全球海潮和中国近海潮汐的负荷效应,建立了中国大陆的精密重力潮汐改正模型.结果表明,采用不同的固体潮模型会对重力潮汐结果产生相对变化幅度小于0.06％的差异;在沿海地区海潮负荷的影响约为整个潮汐的4％,而中部地区约为1％,其中中国近海潮汐模型的影响约占整个海潮负荷的10%,内插或外推潮波的负荷约占海潮负荷的3％.通过比较实测的重力数据表明,本文给出的重力潮汐改正模型的精度远远优于0.5×10^-8 m·s^-2,说明了本文构建的模型的实用性,可为中国大陆高精度重力测量提供有效参考和精密的改正模型.     

Deep electrical conductivity features in the transition zone from the Pacific to Eurasia

Understanding the processes that occur in the transition from the Pacific Ocean to Eurasia is key to constructing the tectonic models of the Earth’s shells and the convection models of the upper mantle. The electromagnetic methods permit estimating the temperature and fluid content (and/or carbon (graphite) content) in the Earth’s interior. These estimates are independent of the traditionally used estimates based on seismic methods because the dependence of electrical conductivity on the physical properties of the rock is based on different principles than the behavior of the elastic waves. The region is characterized by a complicated geological structure with intense three-dimensional (3D) surface heterogeneities, which significantly aggravate the retrieval of the information about the deep horizons in the structure of the Earth’s mantle from the observed electromagnetic (EM) fields. The detailed analysis of the nature of the deep electrical conductivity and structural features of the transition from the Pacific to Eurasia included numerical modeling of the typical two- and three-dimensional models has been carried out. Based on this analysis, the approaches that increase the reliability of the interpretation of the results of the EM studies are suggested.     

Correlation of variations in the thermal neutron flux from the Earth’s crust with the Moon’s phases and with seismic activity

The results of the long-term recording of thermal neutron flux near the Earth’s surface with the use of an unshielded scintillation thermal-neutron detector are presented. The data obtained indicate the presence of periodic variations in the thermal neutron flux with the lunar diurnal and the lunar monthly periods. A hypothesis about the existence in the Earth’s crust of radon-neutron tidal variations in the concentration of thermal neutrons, correlated with the Moon’s phases and which have the gravitational origin, is formulated and confirmed experimentally. A simple mathematical model is proposed, which satisfactorily describes the observed variations. The case of the anomalous behavior of thermal neutrons is presented, which correlates with the high local seismic activity.     

Simulation of tides in hydrocarbon lakes on Saturn’s moon Titan

Numerous hydrocarbon lakes have recently been detected on Saturn’s largest moon Titan, representing the only known large bodies of liquids on a planetary surface outside the Earth. In the context of comparative oceanography, tides and tidal currents in two representative lakes on Titan (Kraken Mare and Ontario Lacus) are simulated by a three-dimensional baroclinic ocean circulation model. Since the tide-generating force on Titan is an order of magnitude larger than on Earth and the gravitational acceleration is small, tides and currents are substantially larger than in Earth’s lakes and are more comparable with those in Earth’s oceans. The predicted maximum tidal range in Kraken Mare is 4 m. The tidal wave propagates around the basin of Kraken Mare, while a nearly standing tidal wave is excited in Ontario Lacus. Titan’s rotation is too slow to affect the tidal flow in any Titan’s lake. The tidal current velocity in Kraken Mare amounts to a few centimeters per second except in the vicinity of a narrow strait, where it is enhanced by an order of magnitude. In summer, when the lake is stratified, internal tides can develop. Seiches cannot be caused by tide. In the largest lakes, atmospheric tide may cause additional lake surface displacements.     

Cosmic ray flux variations, modulated by the solar and earth’s magnetic fields, and climate changes. 1. Time interval from the present to 10–12 ka ago (the Holocene Epoch)

Direct and indirect data on variations in cosmic rays, solar activity, geomagnetic dipole moment, and climate from the present to 10–12ka ago (the Holocene Epoch), registered in different natural archives (tree rings, ice layers, etc.), have been analyzed. The concentration of cosmogenic isotopes, generated in the Earth’s atmosphere under the action of cosmic ray fluxes and coming into the Earth archives, makes it possible to obtain valuable information about variations in a number of natural processes. The cosmogenic isotopes ¹⁴C in tree rings and ¹⁰Be in ice layers, as well as cosmic rays, are modulated by solar activity and geomagnetic field variations, and time variations in these concentrations gives information about past solar and geomagnetic activities. Since the characteristics of natural reservoirs with cosmogenic ¹⁴C and ¹⁰Be vary with climate changes, the concentrations of these isotopes also inform about climate changes in the past. A performed analysis indicates that cosmic ray flux variations are apparently the most effective natural factor of climate changes on a large time scale.     

Radon emanation as an indicator of the intensity of Intergeospheric interactions at the Earth’s crust-atmosphere interface

The results of measurements of the volumetric activity of subsurface radon Rare given for several sections of the Earth’s crust located within the central part of the East European Platform. The influence of tectonic structures on the spatial zonality of radon emanation and the intensity of its temporal variations is shown. The temporal variations in radon emanation are characterized by periodicity with periods of approximately a day, as well as of 4, 14, and 29 days, and, as correlation analysis shows, they are determined by the cyclic recurrence of tidal deformations and by cyclonic processes in the atmosphere. The high correlation of the mean values of the volumetric activity of subsurface radon with the value of atmospheric pressure is established. The dependence of the volumetric activity of subsurface radon and its temporal variations on the depth within the range of depths 0.8–11.5 m is determined. As a result of laboratory experiments, the nonlinear nature of the dependence of radon emanation on the frequency of vibration action is established. For granites, the intensity of radon emanation is maximal at frequencies of 16 and 32 Hz. The results of complex measurements of the microseismic background and the volumetric activity of subsurface radon carried out in the region of the Nelidovo-Ryazan tectonic structure demonstrate the high correlation between the volumetric activity of subsurface radon and the relative amplitude of the quasi-monochromatic component of the microseismic background with a peak frequency of 16.5 Hz.     

中国大陆井水温度潮汐动态的统计与调和分析

用收集到的全国356个井水温度测点的数据, 分析了水温对地球固体潮汐的响应, 统计出 35个存在水温潮汐现象的测点。 利用Baytap-G调和分析方法, 计算了水温潮汐分波的振幅、 振幅比和相位差。 结果表明: 水温潮汐现象是一类较普遍的地球物理现象, 其机制与水位潮汐相关, 可用水动力学模式解释; 水温潮汐变化特征还受太阳辐射热、 含水层和地温的影响, 自流井水温记录潮汐现象的能力高于非自流井、 东部地区水温测点记录潮汐现象的能力高于西部, 与太阳辐射热的影响有关, 在含水层附近的水温测点, 其潮汐动态比其他井段显著, 在受地温影响较大的井段, 水温的潮汐变化幅度与水温梯度成正比; 水温的应力-应变灵敏度量级为0.01～10℃/10^-6m·s^-2。