Reliability of geomagnetic paleointensity data: the effects of the NRM fraction and concave-up behavior on paleointensity determinations by the Thellier method

To test the reliability of the Thellier method for paleointensity determinations, we studied six historic lavas from Hawaii and two Gauss-age lava flows from Raiatea Island (French Polynesia). Our aim is to investigate the effects of the NRM fraction and concave-up behavior of NRM–thermal remanent magnetization (TRM) diagrams on paleointensity determinations. For the Hawaiian samples, the paleointensity results were investigated at both sample and site levels. For consistency and confidence in the paleointensity results, it is important to measure multiple samples from each cooling unit. The results from the Raiatea Island samples confirm that reliable paleointensities can be obtained from NRM–TRM diagrams with concave-up curvature, provided the data are accompanied by successful partial TRM (pTRM) checks and no significant chemical remanent magnetization (CRM) production. We conclude that reliable determinations of the paleofield strength require analyses of linear segments representing at least 40–50% of the total NRM. This new criterion has to be considered for future studies and for evaluating published paleointensities for calculating average geomagnetic field models. Using this condition together with other commonly employed selection criteria, the observed mean site paleointensities are typically within 10% of the Definitive Geomagnetic Reference Field (DGRF). Our new results for the Hawaii 1960 lava flow are in excellent agreement with the expected value, in contrast to significant discrepancies observed in some earlier studies.

Overestimates of paleointensity determinations can arise from cooling-rate dependence of TRM acquisition, viscous remanent magnetization (VRM) at elevated temperatures, and TRM properties of multidomain (MD) particles. These outcomes are exaggerated at lower temperature ranges. Therefore, we suggest that, provided the pTRM checks are successful and there is no significant CRM production, it is better to increase the NRM fraction used in paleointensity analyses rather than to maximize correlation coefficients of line segments on the NRM–TRM diagrams.

We introduce the factor, Q = Nq, to assess the quality of the weighted mean paleointensity, H_w, for each cooling unit.     

Intensity of the geomagnetic field from recent Italian lavas using a new paleointensity method

A new method has been tested on Etna historic lavas for determining the geomagnetic field intensity (F) using the thermoremanent magnetization of volcanic rocks. The procedure involves a number of very short duration heatings above the Curie point, to produce successive laboratory TRM. Thus, it is possible to check the variations in the TRM-acquiring capacity of the samples with the time of heating (t). The curve J = f(t) is then extrapolated towards t = 0, leading to a virtual value of TRM without any laboratory heating, i.e., without the changes that currently occur when the lavas are heated to produce the TRM. Using such virtual values of TRM, satisfactory results of F had been derived from the majority of the samples studied. These results are consistent with Thellier's archeomagnetic data: they show a nearly constant intensity of the geomagnetic field during the last three centuries and, further back into the past, a significant increase of this.     

基于地磁场相对古强度变化的下蜀黄土年代序列及其古气候意义

分布于宁镇地区的下蜀黄土年代标尺研究薄弱.我们以镇江大港钻孔岩芯的下蜀黄土为例,探讨了下蜀黄土记录地磁场相对古强度的可靠性.岩石磁学的实验显示,下蜀黄土大体上符合估计相对古强度对沉积物的要求.以低频磁化率()作为天然剩磁(NRM300)的归一化因子,我们获得了下蜀黄土记录的地磁场相对古强度变化.通过与邻区和全球的单个及合成曲线的对比,我们发现大港钻孔岩芯的相对古强度记录展现出主要的全球偶极场的变化特征,因而可用于建立下蜀黄土的年代标尺.新的年代标尺表明,大港钻孔的磁化率年代序列与北方黄土无法直接对比,证实了该地点的下蜀黄土磁化率变化机制与北方的不同.大港钻孔下蜀黄土的沉积速率与磁化率相关,低风尘沉积速率对应低磁化率,是降水增加所导致的结果.沉积速率与磁化率年代序列显示,下蜀黄土记录的本区季风变化过程可分为4个阶段.阶段Ⅳ(819～700ka)对应中更新世转型阶段,东亚季风降水较少.阶段Ⅲ(700～412ka)对应中更新世的大间冰期,东亚季风降水最多.阶段Ⅱ(412～197ka)时东亚季风降水减少,较阶段Ⅳ略少.阶段Ⅰ(197～34ka)东亚季风降水最少.因此,中更新世气候转型后,长时间尺度东亚季风降水持续减少,可能受全球温度阶段性降低驱动.     

Test determinations of paleointensity in historical lavas of Kamchatka

The reliability of the Thellier method for determining the paleointensity of a geomagnetic field is explored on recent igneous rocks of Kamchatka. The main magnetic mineral in the studied rocks is titanomagnetite with different degree of oxidation. It is obtained that the reliability of the results can be assessed based on the deviations of the check points of the partial thermoremanent magnetization (pTRM) during the Thellier experiment. Besides, for different rocks, it is found that the stability of titanomagnetites to heating during the experiments can be insufficient for validating the reliability of the results of paleointensity determination; however, at the same time, the reliability may depend on the initial (oxidation) state of the magnetic minerals of the studied rocks.     

A discussion of the Alternating Field method and new paleointensity measurements of the Miocene geomagnetic field in Oregon

Summary Thirty eight cores from five Miocene lavas and their underlying baked zones from South Eastern Oregon have provided a test for the Alternating Field method of determining geomagnetic paleointensities, by allowing external consistency tests as well as internal consistency tests. All specimens were run regardless of the reliability tests results. The susceptibility change reliability test is useful as some of the specimens which failed it nevertheless gave misleadingly plausible results, but is shown to be unsufficient, as some specimens which were retained on that test failed to yield valid results. Different specimens from the same core behaved differently during the various experiments, proving that two specimens from a same core cannot reliably be assumed to be identical. A relation between the oxidation state and a change of susceptibility upon heating is seen.The paleofield intensities varied from 0.17±0.04 Oe and 0.31±0.08 Oe, corresponding to Virtual Dipole Moments of 1.23×10²⁵ to 2.15×10²⁵ Oe cm³. These low values are not accompanied by Virtual Geomagnetic Pole instability, but suggest an instability of the strength of the main dipole.     

Solar flare effect on the geomagnetic field and ionosphere

This paper studies the ionospheric and geomagnetic response to an X6.2 solar flare recorded at 14:30 UT on December 13, 2001, in quiet geomagnetic conditions which allow the variations in the geomagnetic field and ionosphere measurements to be easily related to the solar flare radiation.By using measurements from the global positioning system (GPS) and geomagnetic observatories, the temporal evolution of ionospheric total electron content variation, vTEC_V, and geomagnetic field variations, δB, as well as their rates of variation, were obtained around the subsolar point at different solar zenith angles. The enhancement of both parameters was recorded one to three minutes later than the Geostationary Operational Environmental Satellite (GOES) programme recording; such delay tends to depend on the latitude, longitude, and solar zenith angle of the observatory's observations.The vTEC_V is related to the local time and the δB to the intensity and position of the ionospheric currents.The vTEC_V′s maximum value is always recorded later than the maximum values reached by δB and the X-ray intensity. The maximum δB is larger in the local morning than in the afternoon.The rates of vTEC_V and δB have two maximum values at the same time as the maximum values recorded by Hα (for each ribbon).This work shows the quantitative and qualitative relations between a solar flare and the ionospheric and geomagnetic variations that it produces.     

Relative paleointensity of the geomagnetic field during the past 200 ka from the West Philippine Sea and its chronological significance

Our rock magnetic analysis of core Ph05 from the West Philippine Sea demonstrates that the core preserves a strong, stable remanent magnetization and meets the magnetic mineral criteria for relative paleointensity (RPI) analyses. The magnetic minerals in the sequence are dominated by pseudosingle-domain magnetite, and the concentration of magnetic minerals is at the same scale. Both the conventional normalizing method and the pseudo-Thellier method were used in conjunction with the examination of the rock magnetic properties and natural remanent magnetization. Susceptibility (χ), anhysteretic remnant magnetization (ARM) and saturation isothermal remnant magnetization (SIRM) were used as the natural remanent magnetization normalizer. However, coherence analysis indicated that only ARM is more suitable for paleointensity reconstruction. The age model of core is established based on oxygen isotope data and AMS¹⁴C data, which is consistent with the age model estimated from RPI records. The relative paleointensity data provide a continuous record of the intensity variation during the last 200 ka, which correlates well with the global references RPI stacks. Several prominent low paleointensity values are identified and are correlated to the main RPI minima in the SINT-200 record, suggesting that the sediments have recorded the real changes of geomagnetic field. Supported by National Natural Science Foundation of China (Grant No. 90411014) and Pilot Project of the Knowledge Innovation Program of Chinese Academy of Sciences (No. KZCX2-YW-211)     

The intensity of the geomagnetic field at 3.5 Ga: paleointensity results from the Komati Formation, Barberton Mountain Land, South Africa

Paleointensity measurements have been carried out on 3.5 Ga samples from the Komati Formation type locality using both the Thellier and Van Zijl methods. These samples contain a single steeply-directed negative TRM component acquired during metamorphism of the Komati lavas. Thellier experiments yielded values ranging from 12 to 37 μT but an average paleofield intensity for the four best determinations is 20 ± 3 μT. A slightly lower average paleointensity of 15 ± 3 μT was obtained using Van Zijl experiments. Preheating was used to chemically stabilize seven samples used in Van Zijl determinations and these produced nearly ideal plots with an average paleointensity of 13 ± 2 μT. A single basaltic komatiite sample gave a nearly ideal Van Zijl plot indicating about 21 μT, nearly the same paleointensity as the peridotitic komatiite samples even though its NRM intensity was several orders of magnitude lower. Since the Komati characteristic remanence was acquired during a slow cooling, the data must be reduced by a factor of 1.55 to account for the difference between laboratory and natural cooling rates. Calculation of an equivalent equatorial paleointensity using the paleolatitude implied by the steep Komati characteristic remanence then gives value of 5 μT for the intensity of the geomagnetic field at 3.5 Ga, lower than the present value of about 30 μT.     

The effect of solar activity on the secular variation of the geomagnetic field in Romania

The 11-year solar cycle effect in the geomagnetic components H and Z is made clear for Surlari Observatory and 19 repeat stations for the interval 1952–1974. The correlation with Wolf number and its time derivative is discussed in terms of the effects of the external and induced current systems.The $\text{H}\text{?}$ data available for solar cycle 20 (1964–1976) were processed to give the geographical distribution of the secular variation impulse for epoch 1969.5 in Romania. It is suggested that this distribution might reflect the deep internal structure of the area considered.A qualitative correlation is noted between long-period solar activity and variation of the horizontal component of the geomagnetic field at some repeat stations.     

The fine structure of the geomagnetic field

It was shown in our previous works that the dipole part of the geomagnetic field direction is mainly represented by the 1200-yr variation. According to some indications, the residual of the dipole part of the field, the so-called δ variation, may be classified as standing waves. Average values of δ are calculated for each hemisphere. In this work, the difference values equal to the δ variations for each territory minus δ averages over hemispheres (western or eastern) are calculated. The resulting values, fine structure (FS) variations, characterize the FS of the geomagnetic field. The study of the activity of the FS variations and specific features of their N-S and E-W behavior and the comparative analysis of dynamic characteristics (activity and rate) give grounds to classify the FS variations (as well as the δ variations) as standing waves of the nondipole field.     

The effect of observational errors on geomagnetic field models based solely on total-intensity measurements

An investigation is made of the changes produced in the spherical harmonic coefficients of a geomagnetic main field model when values of the total intensity synthesized from the model are systematically perturbed. Corresponding differences between the values of the field components from the original model and those from the perturbed model are also derived. The effects of varying the altitude of the perturbed point and of including a small proportion of component data are discussed. It is shown that significant errors can arise in the region of the dip equator when only scalar data are used, but that these errors can be eliminated if the data set contains 10% of well-distributed component data, some of which refers to the equatorial zone.     

Laboratory verification of submicron magnetite production in pseudotachylytes: relevance for paleointensity studies

Pseudotachylytes generally possess stable remanent magnetizations but the processes by which pseudotachylytes are magnetized remain poorly understood. Magnetic hysteresis and scanning electron microscope studies reveal that experimental frictional melting of granites produces dispersed submicron inclusions of weakly interacting pseudo-single-domain (PSD) magnetite, in artificial pseudotachylyte. The magnetite inclusions are absent in the undeformed granite protolith and result from oxidation of Fe in melt-susceptible mafic minerals during the melt-quenched event. The pseudotachylytes acquired a stable thermal remanence in fine-grained PSD magnetites during the rapid cooling of the melt, implying that fine-grained magnetite has the potential for paleointensity determinations of contemporaneous magnetic fields with co-seismic faulting in granitoids.     

地磁正常场的选取与地磁异常场的计算

根据2003年中国地磁观测数据（包括135年地磁测点和35个地磁台）以及我国邻近地区38个IGRF计算点的地磁数据,计算中国地磁异常场的分布。选取两种地磁场模型作为地磁正常场,一是国际参考地磁场的球谐模型,二是中国地磁场泰勒多项式模型。根据各个测点的地磁异常值（观测值减去模型计算值）,用球冠谐分析方法计算地磁异常场的球冠谐模型,并绘制2003年中国地磁异常（△D,△I,△F,△X,△Y,△Z）。分析和讨论了中国地磁异常场。     

Spatial power spectra of the crustal geomagnetic field and core geomagnetic field

Lowes (1966, 1974) has introduced the function R_n defined by $\text{R}{}^{\mathrm{n}}\text{=(n + 1)}\text{∑}\text{m=0}\text{∞}\text{[(g}{}^{\mathrm{m}}{}_{\mathrm{n}}\text{)}{}^2\text{+ (h}{}^{\mathrm{m}}{}_{\mathrm{n}}\text{)}{}^2\text{]}\text{where}\text{g}{}_{\mathrm{n}}{}^{\mathrm{m}}\text{and}\text{h}{}_{\mathrm{n}}{}^{\mathrm{m}}$ are the coefficients of a spherical harmonic expansion of the scalar potential of the geomagnetic field at the Earth's surface. The mean squared value of the magnetic field B = ??V on a sphere of radius r > α is given by $\text{〈}\text{B}\text{·〉 =}\text{∑}\text{n=1}\text{∞}\text{R}{}^{\mathrm{n}}\text{(}\text{a/}\text{r}\text{)}{}^{\mathrm{2n=4}}\text{where}\text{a}$ is the Earth's radius. We refer to R_n as the spherical harmonic spatial power spectrum of the geomagnetic field.In this paper it is shown that Rⁿ = R^M_n = R^C_n where the components R_n^M due to the main (or core) field and R_n^C due to the crustal field are given approximately by $\text{R}{}^{\mathrm{M}}{}_{\mathrm{n}}\text{= [}\text{(n =1)/}\text{(n + 2)}\text{](1.142 × 10}{}^9\text{)(0.288}{}^{\mathrm{n}}\text{Λ}{}^2\text{R}{}^{\mathrm{C}}{}_{\mathrm{n}}\text{= [}\text{(n =1){[1 — exp(-n/290)]/}\text{(n/290)} 0.52 Λ}{}^2\text{where}\text{Iγ = 1 nT}$. The two components are approximately equal for n = 15.Lowes has given equations for the core and crustal field spectra. His equation for the crustal field spectrum is significantly different from the one given here. The equation given in this paper is in better agreement with data obtained on the POGO spacecraft and with data for the crustal field given by Alldredge et al. (1963).The equations for the main and crustal geomagnetic field spectra are consistent with data for the core field given by Peddie and Fabiano (1976) and data for the crustal field given by Alldredge et al. The equations are based on a statistical model that makes use of the principle of equipartition of energy and predicts the shape of both the crustal and core spectra. The model also predicts the core radius accurately. The numerical values given by the equations are not strongly dependent on the model.Equations relating average great circle power spectra of the geomagnetic field components to R_n are derived. The three field components are in the radial direction, along the great circle track, and perpendicular to the first two. These equations can, in principle, be inverted to compute the R_n for celestial bodies from average great circle power spectra of the magnetic field components.     

Direction and relative paleointensity of the geomagnetic field in the middle and later pleistocene from paleomagnetic data of the Roxolany section (Ukraine)

The lower part of the Roxolany section (Ukraine) is studied, and paleomagnetic and petromagnetic characteristics of rock samples are obtained for the time interval 300-75 ka. Detailed curves of variations in angular parameters of the geomagnetic field are constructed in the entire time interval, and curves of variations in the relative paleointensity are obtained for the interval 300-180 ka. Using the values of the angular parameters and the deviations of the virtual geomagnetic pole from the position of the stationary field, anomalous directions compatible with the Jamaica excursion are identified in the intervals 250-249 and 221-220 ka. The geomagnetic field evolution is studied by methods of wavelet analysis, and the field generation process is shown to vary in the interval 300-180 ka.     

第10代国际地磁参考场

国际地磁与超高层大气物理学协会（IAGA）于2004年12月12日发表了第10代国际地磁参考场（IGRF）。这是对地球主磁场标准数学描述的最新版本，在对地球深部、地壳、电离层及磁层的研究中有广泛应用。国际地磁参考场系数由国际地磁与超高层大气物理学协会的一个特别工作组最终确定下来。国际地磁参考场是地磁场模型工作者和全世界从事卫星及台站地磁观测数据收集和发布的单位共同努力的结果。