


编号 5 (2023)
Articles
Theory of Anhysteretic Remanent Magnetization of Single-Domain Grains
摘要
Abstract—A new approach to the solution of kinetic equations describing the process of formation of anhysteretic remanent magnetization (ARM) is proposed, which made it possible to accelerate the numerical calculation of the process of formation of ARM by two orders of magnitude for uniaxial oriented non-interacting single-domain particles, while practically not yielding in accuracy to a strict numerical solution. It follows from the calculation results that the susceptibility of ARM is entirely determined by the magnitude of the particle’s coercivity parameter. The data of the previous approximate calculations of ARM value are compared with the exact results presented here and it is shown that the discrepancy between the exact data and the approximate estimates increases with the growth of g, but remains relatively small, within 23%. The proposed algorithm for the rapid calculation of kinetic equations allows us to analyze with physical rigor the method of pseudo-Thellier estimation of paleointensivity for an ensemble of single-domain particles, which is supposed to be done in subsequent works.



Determination of the Effective Electrical Conductivity of a Fluid–Saturated Core from Computed Tomography Data
摘要
Abstract—This paper proposes a technique for determining the effective specific electrical conductivity of rock samples when their digital models are used. A modified algorithm for reconstructing the internal structure of the sample from the core’s nondestructive imaging data can be used to construct a relevant discrete model that approximates the pore space with a high degree of accuracy. Unlike existing approaches, the reconstructed discrete geometric model of a heterogeneous medium sample is hierarchical and oriented to the application of parallel computational schemes of multiscale finite element methods for a forward mathematical simulation of electromagnetic processes. The paper presents the results of solving the problem of determining the effective specific electrical conductivity of fluid–saturated rock samples and compares them with the data from laboratory experiments.



An Ultralow Geomagnetic Field Intensity in The Mesoproterozoic Based on Studies of 1380 Ma Old Intrusive Bodies from the Udzha Aulacogen of the Siberian Platform
摘要
Abstract—Estimating the time of formation of the Earth’s solid inner core is a first–order problem in the thermal evolution of the Earth as a planet, which can be solved in particular by paleointensity (Banc) determinations. For this purpose, we have studied a collection of ~1380 Ma rocks sampled in the Udzha River valley within the Udzha aulacogen from the dolerite sill near the Khapchanyr River. The sill is an intrusion 5–7 meters thick that cuts through the siltstones and carbonates of the Udzha Formation. To obtain reliable paleointensity (Banc) determinations, magnetic and thermomagnetic properties of the studied rocks have been investigated, and their X–ray diffraction patterns and electron–microscopic images have been collected. It is shown that the carriers of the characteristic component of natural remanent magnetization are single–domain (SD) and small pseudo–single–domain (PSD) magnetite grains. Two methods have been used to determine the Banc: the Thellier–Coe procedure including reheating to lower temperatures (the pTRM checkpoints procedure) and the Wilson express method. Paleointensity determinations have been obtained for 9 samples (30 duplicates) that met the selection criteria. The mean value of Banc is extremely low and amounts to 4.54 ± 0.49 μT; the corresponding calculated value of the virtual dipole moment (VDM) is (11.1 ± 1.2) × 1021 Am2, which is almost an order of magnitude lower than the mean VDM in modern epoch (≈80 × 1021 Am2) and more than six times lower than the mean VDM in the Cenozoic (6.44 × 1022 Am2). The VDM data over the 350–3500 Ma interval presented in the world paleointensity database (WPD), Borok, have been analyzed. An alternation of periods of low and high paleointensity is observed in the Precambrian and Paleozoic, which indicates a large variability in the operation mode of the geomagnetic dynamo regardless of the existence or absence of an inner solid core of the Earth. It is important to note that the number of reliable VDM values over such a long interval, 350–3500 Ma, is too small for a complete statistical analysis to single out any time interval as the most probable for the formation of the inner core.



Reconstructions of Paleowind Directions in the Pleistocene: Evidence from the Anisotropy of Magnetic Susceptibility of the Loess–Paleosol Series of Tajikistan and the Azov Region
摘要
Abstract—The magnetic fabric of the reference Khonako-II (south Tajikistan), Beglitsa, and Chumbur-Kosa (Azov region) loess-paleosol sections was studied to reconstruct the paleowind directions during the Middle and Upper Pleistocene. The paleosols and loesses of the Khonako-II section represent a promising object for study of paleowinds. Our data indicate a prevailing ~ northwestward or ~northeastward wind directions during accumulation of the upper 40.8 m of the section. Thereby, a change from warm (pedocomplex 2) to cold (loess 2) epochs during the accumulation of the lower half of the studied sequence (end of the Middle Pleistocene) was accompanied by “switching” the predominant winds (from northwestern to northeastern ones), whereas accumulation of the upper part of the sequence was mainly controlled by winds of intermediate directions. The latter can be considered as evidence for the lesser contrast of wind rose during later epochs that were responsible for the formation of the upper part of the sequence. The study revealed that the magnetic fabric of loess and paleosol horizons of the Azov region sequences was deformed and cannot be used for paleowind reconstruction. It was demonstrated that sampling using plastic container technique widely applied in studying loess–paleosol deposits could led to the significant deformation of primary magnetic fabric.



Staging of Occurrence of Seismicity Anomalies before Earthquakes in Kamchatka, Japan and Iceland
摘要
Abstract—The paper presents the results of a study showing that anomalies in the seismic regime parameters before earthquakes of various magnitudes occur in stages. The occurrence in stages means the correlation between the times of formation and development of anomalies in various seismic regime parameters. Earthquakes in regions with two general types of tectonics are selected for analysis: in the subduction zone (Kamchatka and Japan) and in the rift zone (Iceland). The selection of regions is primarily based on the availability and quality of regional seismic catalogs. GR b-value and the composite parameter known as the RTL are used as the seismic regime parameters. The detection of spatiotemporal anomalies before the selected earthquakes is based on the known “precursory patterns” of the seismic regime parameters. Comparing the durations of the detected anomalies shows that the anomalies of b-value generally occur earlier than the RTL anomalies. Possible reasons why the anomalies occur in stages are suggested. In the vicinity of the studied earthquakes, a change in the seismogenic rupture concentration parameter within the corresponding seismic cycles is also estimated. Comparing the times at which the detected seismic regime anomalies occur with the values of the seismogenic rupture concentration parameter corresponding to these times shows that the formation of seismic regime anomalies occurs at a stage when the system of seismogenic ruptures accumulated during the seismic cycle has almost reached its critical value.



The Khubsugul Earthquake of January 12, 2021, ML = 6.9, in the Seismicity Structure of the Tuva-Mongolian Block
摘要
Abstract
—The paper presents the studies of the Khuvsgul earthquake on January 11, 2021 at 21:32 UTC (January 12, 2021 at 05:32 local time), MW = 6.7, ML = 6.9, and the seismicity structure in the aftershock period for the Altai-Sayan mountain region and the Baikal rift zone, where the epicenter of this earthquake was located. Two faults are seismically activated, diverging from the southern end of the aftershock area at an acute angle: one in the northeast and one in the northwest direction, as well as transverse faults between them. According to the epicenter position and studies of the source area by other authors, the main event corresponds to the northeastern fault, and large aftershocks occurred at the junction of the northwestern fault with transverse faults feathering from the east. The main event was immediately followed by a series of large aftershocks, the strongest of which occurred on March 31, 2021 with ML = 6.2 and on May 3, 2021 with ML = 6.4. Spatial changes in the seismic regime of the aftershock region led to the predominant activity of its southern end. The junction area of the collisional structures of the Altai-Sayan folded zone and the rift structures of the Baikal depressions system is distinguished in seismicity as a block structure with increased seismicity near the block boundaries. First of all, these are the Tuva-Mongolian block and the eastern part of the Sayano-Tuva block. After the Khuvsgul earthquake of 2021, a block structure with the activation of the epicentral zones of the 1991 Busingol earthquake, the 2011–2012 Tuva earthquakes, and other structures seismically active until 2021 has increased seismic activity. It is proved that the 2014 Khubsgul earthquake occurred under the basin of the same name and is associated with other faults than the 2021 earthquake and is not a direct precursor of the 2021–2022 activation.



On an Analytical Representation of an Integral Related to the Fock Integral That Appears in Calculations of the Electromagnetic Fields of Dipole Sources at the Interface between Two Half-Spaces
摘要
Abstract—
The Fock integral is called after Fock who introduced it for the theoretical analysis of the electromagnetic field of magnetic dipoles at the boundary of a uniform conducting (nonmagnetic) half-space and obtained its analytical expression in terms of cylindrical functions. Detailed analytical representations of integrals, where all components of the fields of the vertical and horizontal magnetic dipoles are expressed, are reported in [A.V. Veshev et al., 1971]. To obtain analytical expressions for similar integrals representing the components of the fields of electric dipoles in a similar model, it is necessary to consider not only the Fock integral but also another related integral conditionally called the Fock integral 1 whose analytical expression is still unknown. The aims of this work are to derive an inhomogeneous linear first-order differential equation for this integral with the corresponding boundary conditions and to obtain the analytical representation of the Fock integral 1 by solving this equation. The result of this work will allow one to simplify the simulation of fields in a uniform half-space and to improve the interpretation of electromagnetic data due to more accurate and reliable estimates of the normal field in such models of a host medium.



Geoelectric Interpretation of the Magnetic Monitoring Data in the Southwest Part of Baikal
摘要
Abstract—
The long-term gradient magnetic measurements performed in the southwest part of Baikal as supplementing the experiment on deep-water monitoring of the electric field vertical component, can be applied to study the deep geoelectric section. Baikal cannot be approximated by a two-dimensional model in the region under study. Therefore, a three-dimensional geoelectric model of the lake and adjacent territories was constructed based on the existing geological and geophysical data. However, these data are insufficient for the region where our observations are carried out. Refinement of the section upper part over the magnetic field gradients on long bases has been performed using the neural network inversion method. Then, proceeding from the obtained basic model, alternative models were developed, which reflect the known hypotheses about the structure of the deep part of the Baikal rift section: mantle diapir and asthenospheric upwelling. The phase differences between the base-averaged field gradient and the field itself at the base point, as well as the phase differences on the ends of this base, as the most noise-proof gradient functions, are calculated for them. Based on the comparison of the experimental and model frequency curves, it is concluded that the mantle diapir fits best the experimental data, at least, in the southwest part of the Baikal rift zone.



Deep Electrical Conductivity Anomalies in the Chaun Bay Region Based on Magnetic Variation Sounding Data
摘要
Abstract—This paper discusses the results of magnetic variation soundings at two sites in the eastern Arctic, in the Chaun Bay region, at the Pevek and Valkarkai weather stations, by using the ModEM program to perform a three–dimensional (3D) inversion of tippers. The inversion has produced a geoelectric model of the region in a subsurface area of 300 × 300 × 200 km. The moduli of tippers at both sites have values between 0.2 and 0.85, reaching the maximum ones in a period of 1000 s. At the Pevek weather station, the real induction arrow in the Parkinson convention is oriented to the west throughout the range of periods. At the second site, its azimuth changes from 30° to the NE to –30° to the NW as the period of variations increases. Throughout the range of depths, conductive inhomogeneities are located to the west and north of the Chaun folded zone. They extend as a narrow strip under the western and northern coastal parts of the zone at depths of 8–12 km. As the depth increases, they are split into blocks, which reach their maximum size in the horizontal plane at depths ranging between 20 and 30 km. The most prominent among them is the conductive block beneath the southern part of the Chaun Bay. The roots of these blocks are distinguishable at depths of up to 100 km. The Chaun folded zone is represented by a high–resistance block down to a depth of 150 km with an electrical resistivity of more than 1000 Ω⋅m. Comparison of the obtained geoelectric model with the geophysical studies previously conducted here reveals a correlation between the location of conductive formations and the location of weakly positive anomalies in the gravity field in the Bouguer and isostatic reductions in the coastal water area.



Joint Three-Dimensional Interpretation of AMTS and RMT-C Tensor Data in a Region Promising the Discovery of Primary Diamond Sources (Karelian Isthmus)
摘要
Abstract
—The paper presents the techniques and results of the joint interpretation of AMTS and RMT-C tensor data in the area of the Karelian Isthmus near the village of Yablonovka, promising the discovery of primary diamond sources. Earlier, according to drilling data, fluid-cataclastic breccias, fluidolites, were identified here, which, according to modern geological concepts, can be primary diamond sources. The site is located in the marginal part of the Pash-Ladoga rift zone and is characterized by a complex three-dimensional structure. Under these conditions, in the presence of a system of multidirectional faults, the joint interpretation of the AMTS and RMT-C tensor data was performed within a 3D model of the medium. To reduce the inversion execution time, an approach and a procedure for transforming the RMT-C data obtained in the intermediate zone of a controlled source to sounding curves corresponding to the plane wave model are proposed. The transformed RMT-C apparent resistivity curves were also used to correct the static shifts of the AMTS curves. A joint 3D inversion was performed taking into account the relief. The resulting model is in good agreement with the geological data obtained from the wells previously drilled near the work site. The geoelectric sections show zones of relatively increased resistivity, which may be associated with fluidolites. The developed approaches to the joint 3D interpretation of AMTS and RMT-C tensor data and the obtained information about the structure and properties of the basement rocks, rocks hosting fluidolites, and overlying sedimentary deposits, can be used in prospecting works for diamonds on the Karelian Isthmus and adjacent territories.



PMTOOLS: New Application for Paleomagnetic Data Analysis
摘要
Abstract
—This paper introduces PMTools (https://pmtools.ru), a novel cross-platform open-source web application designed for the analysis of paleomagnetic data. Our software offers a user-friendly interface and supports the most commonly used data formats in paleomagnetism, including .pmd, .dir, .pmm, .rs3, .squid, .vgp, and .gpml. It encompasses all the necessary functionalities for conducting principal component analysis of natural remanent magnetization and computing mean paleomagnetic directions along with corresponding virtual geomagnetic (paleomagnetic) poles. PMTools facilitates the application of various paleomagnetic tests to user data, such as the fold test, conglomerate test, and reversals test. Notably, the application allows for seamless import and export of paleomagnetic data (tables) in formats compatible with .csv and MS Excel. Additionally, all generated figures can be exported as high-quality vector graphics in .svg format, specifically designed for direct integration into publications and presentations. Furthermore, PMTools enables the export of sets of paleomagnetic poles in GMap and GPlates software formats, enabling researchers to promptly use their paleomagnetic data for paleotectonic reconstructions. PMTools has an intuitive interface, customizable hotkeys, and an extensive array of graph elements, all of which contribute to making PMTools an appealing, cutting-edge tool for processing and analyzing the results of paleomagnetic studies.



ДИСКУССИИ
Magnetic Effects of the Double Earthquake of March 16, 2022 (Japan). Observation Results
摘要
Abstract—The magnetic effects of strong earthquakes are considered using the example of a unique event: two underwater earthquakes of magnitude 6.0 and 7.3 that occurred on March 16, 2022 at close times (at 14:34 and 14:36 UTC, respectively) with a distance between sources of ~11 km. With the use of data from the Mikhnevo geophysical observatory of the IDG RAS (MHV) and a number of magnetic observatories of the INTERMAGNET international network, it is shown that in the absence of global geomagnetic disturbances, earthquakes were accompanied by characteristic variations of the Earth’s magnetic field recorded as a train of alternating geomagnetic field oscillations with an amplitude of ~2–8 nT and a period of ~30 min almost simultaneously at ~15:30 UTC by all magnetic observatories located at distances from 210 to ~13 000 km, as well as by variations in the periods of passage of the seismic and infrasonic signal in the MHV.



On the Magnetic Effects due to the Earthquake of March 16, 2022 in Japan
摘要
Abstract
—The magnetic effects of two similar underwater earthquakes with magnitudes of 6.0 and 7.3, which occurred on March 16, 2022, were considered in [Adushkin et al., 2023]. According to the data of INTERMAGNET magnetic observatories, these earthquakes were found to be accompanied (with a delay of ~55 min) by variations in the Earth’s magnetic field in the form of a train of quasi-periodic oscillations with an amplitude of ~2–8 nT and a period of ~30 min at distances of ~210 to ~3000 km from the epicenter. It was suggested in the aforementioned study that this magnetic effect is caused by a perturbation of the geodynamo as a result of the impact of seismic waves propagating deep into Earth. This interesting hypothesis requires a detailed discussion from different points of view. A more detailed analysis of the pattern of geomagnetic field perturbation at all latitudes, performed by us, leads to a conclusion that the found quasi-periodic perturbation is a mid-latitude response to auroral electrojet variations and is not related to the earthquake. According to our estimates, variations with a source at the core–mantle interface on a time scale less than 1 year cannot manifest themselves on the Earth’s surface at all.


