ALBERT

Description: Seismic waves sensitive to the outermost part of the Earth's liquid core seem to be affected by a stably stratified layer at the core–mantle boundary. Such a layer could have an observable signature in both long-term and short-term variations of the magnetic field of the Earth, which are used to probe the flow at the top of the core. Indeed, with the recent SWARM mission, it seems reasonable to be able to identify waves propagating in the core with period of several months, which may play an important role in the large-scale dynamics. In this paper, we characterize the influence of a stratified layer at the top of the core on deep quasi-geostrophic (Rossby) waves. We compute numerically the quasi-geostrophic eigenmodes of a rapidly rotating spherical shell, with a stably stratified layer near the outer boundary. Two simple models of stratification are taken into account, which are scaled with commonly adopted values of the Brunt–Väisälä frequency in the Earth's core. In the absence of magnetic field, we find that both azimuthal wavelength and frequency of the eigenmodes control their penetration into the stratified layer: the higher the phase speed, the higher the permeability of the stratified layer to the wave motion. We also show that the theory developed by Takehiro & Lister for thermal convection extends to the whole family of Rossby waves in the core. Adding a magnetic field, the penetrative behaviour of the quasi-geostrophic modes (the so-called fast branch) is insensitive to the imposed magnetic field and only weakly sensitive to the precise shape of the stratification. Based on these results, the large-scale and high-frequency modes (1–2 month periods) may be detectable in the geomagnetic data measured at the Earth's surface, especially in the equatorial area where the modes can be trapped.

Description: A 3-D magnetotelluric (MT) inversion code using unstructured tetrahedral elements has been developed in order to correct the topographic effect by directly incorporating it into computational grids. The electromagnetic field and response functions get distorted at the observation sites of MT surveys because of the undulating surface topography, and without correcting this distortion, the subsurface structure can be misinterpreted. Of the two methods proposed to correct the topographic effect, the method incorporating topography explicitly in the inversion is applicable to a wider range of surveys. For forward problems, it has been shown that the finite element method using unstructured tetrahedral elements is useful for the incorporation of topography. Therefore, this paper shows the applicability of unstructured tetrahedral elements in MT inversion using the newly developed code. The inversion code is capable of using the impedance tensor, the vertical magnetic transfer function (VMTF), and the phase tensor as observational data, and it estimates the subsurface resistivity values and the distortion tensor of each observation site. The forward part of the code was verified using two test models, one incorporating topographic effect and one without, and the verifications showed that the results were almost the same as those of previous works. The developed inversion code was then applied to synthetic data from a MT survey, and was verified as being able to recover the resistivity structure as well as other inversion codes. Finally, to confirm its applicability to the data affected by topography, inversion was performed using the synthetic data of the model that included two overlapping mountains. In each of the cases using the impedance tensor, the VMTF and the phase tensor, by including the topography in the mesh, the subsurface resistivity was determined more proficiently than in the case using the flat-surface mesh. Although the locations of the anomalies were not accurately estimated by the inversion using distorted impedance tensors due to the slightly undervalued gain, these locations were correctly estimated by using undistorted impedance tensors or adding VMTFs in the data. Therefore, it can be concluded that the inversion using the unstructured tetrahedral element effectively prevents the misinterpretation of subsurface resistivity and recovers subsurface resistivity proficiently by representing the topography in the computational mesh.

Description: Separating the contribution of different hematite coercivity grains to the magnetic fabric is a standing problem in rock magnetism because of the common occurrence of thermochemical alterations when measuring the anisotropy of thermal remanence. A technique that eliminates this bias is presented, which is useful when there is a need to separate the fabric of detrital from pigmentary hematite, for example. The method is based on stepwise thermal demagnetization of saturation isothermal remanent magnetizations (IRMs) applied orthogonally on three sister specimens, allowing calculation of the anisotropy tensor from the three components of each demagnetized IRM vector, avoiding the necessity of having to apply IRMs to thermochemically altered specimens. Vector subtraction allows determining the anisotropy tensor for specific unblocking-temperature ranges. The anisotropies of the pigmentary, specular and total hematite of the Mauch Chunk Formation red beds of Pennsylvania have been measured from an oriented block sample and results are compared to previous anisotropy measurements performed using the high-field anisotropy of isothermal remanence technique (hf-AIR), which measures total undifferentiated hematite. Experiments were conducted using non-saturating 1 T and fully saturating 5.5 T fields: both experimental sets seem capable of measuring the orientation of the specularite anisotropy principal axes, but 5.5 T are needed to capture the orientation of the higher coercivity pigmentary grains. The magnitudes of the principal axes, instead, are only faithfully measured using 5.5 T fields and yield somewhat higher anisotropies than those measured by hf-AIR. The fundamental requirement for this technique is homogeneous material among the three sister specimens, which is a significant limitation; homogeneity tests allow assessment of applicability of the method and reliability of the results.

Description: Archaeomagnetic field models cover longer timescales than historical models and may therefore resolve the motion of geomagnetic features on the core–mantle boundary (CMB) in a more meaningful statistical sense. Here we perform a detailed appraisal of archaeomagnetic field models to infer some aspects of the physics of the outer core. We characterize and compare the identification and tracking of reversed flux patches (RFPs) in order to assess the RFPs robustness. We find similar behaviour within a family of models but differences among different families, suggesting that modelling strategy is more influential than data set. Similarities involve recurrent positions of RFPs, but no preferred direction of motion is found. The tracking of normal flux patches shows similar qualitative behaviour confirming that RFPs identification and tracking is not strongly biased by their relative weakness. We also compare the tracking of RFPs with that of the historical field model gufm1 and with seismic anomalies of the lowermost mantle to explore the possibility that RFPs have preferred locations prescribed by lower mantle lateral heterogeneity. The archaeomagnetic field model that most resembles the historical field is interpreted in terms of core dynamics and core–mantle thermal interactions. This model exhibits correlation between RFPs and low seismic shear velocity in co-latitude and a shift in longitude. These results shed light on core processes, in particular we infer toroidal field lines with azimuthal orientation below the CMB and large fluid upwelling structures with a width of about 80° (Africa) and 110° (Pacific) at the top of the core. Finally, similar preferred locations of RFPs in the past 9 and 3 kyr of the same archaeomagnetic field model suggest that a 3 kyr period is sufficiently long to reliably detect mantle control on core dynamics. This allows estimating an upper bound of 220–310 km for the magnetic boundary layer thickness below the CMB.

Description: Post-depositional reductive diagenesis usually results in partial or entire cleansing of the pristine palaeomagnetic signal, therefore, its intensity is important to be assessed for sediments that are in the purpose of retrieving palaeomagnetic information. Grain size, rock magnetic and geochemical studies on the entire core, along with scanning electron microscope observations and X-ray diffraction analyses for representative samples were carried out on a Holocene sediment core retrieved from the deep water part of Huguangyan maar lake (HGY), southeast China. The pristine magnetic mineral assemblage of the studied core is domianted by superparamagnetic (SP) and stable single domain titanomagnetite, and high coercivity minerals are not detectable. Based on down-core variations of the average grain size ( M Z ), total organic carbon (TOC), detrital elements (Al, Ti, Fe and Mn) and the concentration and mineralogy of magnetic minerals, the studied core could be divided into three subsections. The uppermost subsection is the least affected by diagenesis, with detrital titanomagnetite as the dominant magnetic mineral. This is owing to low TOC contents, but high detrital input generated by weak Asian summer monsoon intensity during the late Holocene. The intermediate subsection shows down-core progressively enhanced dissolution of detrital titanomagnetite, and concomitant formation of authigenic pyrite and siderite, which indicates down-core progressively enhanced diagenesis generated by down-core progressive increasing TOC content, but decreasing detrital input as the result of down-core progressively strengthened Asian summer monsoon intensity. The pristine magnetic mineral assemblage has been profoundly modified in the lowermost subsection. At certain positions of the lowermost subsection, detrital titanomagnetite has been even completely dissolved via diagenesis, giving place to authigenic pyrite and siderite. High TOC content, but low detrital input generated from strong Asian summer monsoon intensity during the early Holocene are accountable for intensive diagenesis in the lowermost subsection. Complete erasing of detrital magnetic input signal at certain positions of the lowermost subsection, and considerable formation of authigenic siderite indicate that palaeomagnetic records of the studied core have been significantly compromised. The studied core has relatively higher TOC content, lower detrital matter content, calmer sedimentary environments, and less DO available at its water–sediment interface than the cores retrieved at relatively shallower water depths, which all contribute to its relatively stronger diagenesis. Progressive thickening of the upper two subsections with increasing water depth is owing to progressive increase in sedimentation rate with increasing water depth, which is the key factor in determining the thickness of each diagenetic subsection of cores from HGY. It would be better that lake sediments for palaeomagnetic investigations collected at a water depth shallower than the depth of its thermocline.

Description: The rocks in the crust and the upper mantle of the Earth are believed to exhibit electrical anisotropy to some extent. It is beneficial to further understand and recognize the propagation of the electromagnetic waves in the Earth by investigating the magnetotelluric (which is one of the main geophysical techniques to probe the deep structures in the Earth) responses of the media with anisotropic conductivity structures. In this study, we examine the magnetotelluric fields over an idealized 2-D model consisting of two segments with axially anisotropic conductivity structures overlying a perfect conductor basement by a quasi-static analytic approach. The resulting analytic solution could not only contribute to the electromagnetic induction theory in the anisotropic Earth but also serve as at least an initial standard solution which could be used to validate the reliability and accuracy of the numerical algorithms developed for modelling the magnetotelluric responses of the 2-D media with much more general anisotropic conductivity.

Description: Several works have reported that haematite has non-linear initial susceptibility at room temperature, like pyrrhotite or titanomagnetite, but there is no explanation for the observed behaviours yet. This study sets out to determine which physical property (grain size, foreign cations content and domain walls displacements) controls the initial susceptibility. The performed measurements include microprobe analysis to determine magnetic phases different to haematite; initial susceptibility (300 K); hysteresis loops, SIRM and backfield curves at 77 and 300 K to calculate magnetic parameters and minor loops at 77 K, to analyse initial susceptibility and magnetization behaviours below Morin transition. The magnetic moment study at low temperature is completed with measurements of zero field cooled–field cooled and AC susceptibility in a range from 5 to 300 K. The minor loops show that the non-linearity of initial susceptibility is closely related to Barkhausen jumps. Because of initial magnetic susceptibility is controlled by domain structure it is difficult to establish a mathematical model to separate magnetic subfabrics in haematite-bearing rocks.

Description: A series of important geological events occurred in the Tibetan Plateau area during the Jurassic, such as the collision of the Lhasa and Qiangtang Terranes, the closure of the Meso-Tethyan Ocean, the opening of the Neo-Tethyan Ocean and the cessation of the mega-monsoon. The ~3000 m thick Jurassic sedimentary sequence in the Qiangtang Basin on the central Tibetan Plateau, which is called the Yanshiping (YSP) Group, recorded these geological events. However, the chronology of the sequence is surprisingly poorly constrained. Here, we perform a detailed palaeomagnetic analysis on the ~1060 m thick middle and upper portions of the YSP Group (the Xiali and Suowa Formations) in the YSP section of the eastern Qiangtang Basin. Three bivalve zones at stratigraphic intervals of ~40–140, 640–800 and 940–1040 m are identified, which yield a Bathonian–Callovian age for the Lower Xiali Fm., a Callovian–Oxfordian age for the Lower Suowa Fm. and an Oxfordian–Kimmeridgian age for the Upper Suowa Fm. A total of 544 oriented palaeomagnetic samples were collected from the section. By combining thermal and alternating field demagnetizations, clear characteristic remanent magnetization (ChRM) directions are isolated for most of the samples. The robust ChRM directions pass fold and reversals tests, which support the primary nature of the ChRMs and yield a palaeopole at 76.8°N/297.2°E (dp = 2.2°, dm = 3.7°). A total of 27 normal and 26 reversed polarity zones were successfully recorded in the section. Combined with fossil age constraints, results suggest that the section is plausibly composed of a Callovian-Early Kimmeridgian age sedimentary sequence.

Description: This study investigates pressure effects on the magnetic properties of non-interacting single-domain (SD) magnetite. Using a high-pressure cell specially designed for a Magnetic Property Measurement System, magnetic hysteresis measurements were conducted under high pressures of up to 1 GPa on natural plagioclase crystals containing much acicular SD magnetite. Coercivity and saturation magnetization were nearly constant with pressure, while saturation remanent magnetization and coercivity of remanence decreased with pressure at moderate rates of –8 per cent GPa –1 and –18 per cent GPa –1 , respectively. These results suggest that temperature effects govern the magnetic behaviour of acicular SD magnetite grains in the middle and lower crusts.

Description: Characterization of geomagnetic field behaviour on timescales of centuries to millennia is necessary to understand the mechanisms that sustain the geodynamo and drive its evolution. As Holocene paleomagnetic and archeomagnetic data have become more abundant, strategies for regularized inversion of modern field data have been adapted to produce numerous time-varying global field models. We evaluate the effectiveness of several approaches to inversion and data handling, by assessing both global and regional properties of the resulting models. Global Holocene field models cannot resolve Southern hemisphere regional field variations without the use of sediments. A standard data set is used to construct multiple models using two different strategies for relative paleointensity calibration and declination orientation and a selection of starting models in the inversion procedure. When data uncertainties are considered, the results are similar overall regardless of whether we use iterative calibration and reorientation, or co-estimation of the calibration and orientation parameters as part of the inversion procedure. In each case the quality of the starting model used for initial relative paleointensity calibration and declination orientation is crucial and must be based on the best absolute information available. Without adequate initial calibration the morphology of dipole moment variations can be recovered but its absolute value will be correlated with the initial intensity calibrations, an effect that might be mitigated by ensuring an appropriate fit to enough high quality absolute intensity data with low uncertainties. The declination reorientation mainly impacts regional field structure and in the presence of non-zonal fields will result in a non-zero local average. The importance of declination orientation is highlighted by inconsistencies in the West Pacific and Australian sediment records in CALS10k.1b model. Great care must also be taken to assess uncertainties associated with both paleomagnetic and age data and to evaluate the effects of poor data distribution. New consistently allocated uncertainty estimates for sediment paleomagnetic records highlight the importance of adequate uncertainties in the inversion process, as they determine the relative weighting among the data and overall normalized misfit levels which in turn influence the complexity of the inferred field models. Residual distributions suggest that the most appropriate misfit measure is the L 1 norm (minimum absolute deviation) rather than L 2 (least squares), but this seems to have relatively minor impact on the overall results. For future Holocene field modelling we see a need for comprehensive methods to assess uncertainty in individual archeomagnetic data so that these data or models derived from them can be used for reliable initial relative paleointensity calibration and declination orientation in sediments. More work will be needed to assess whether co-estimation or an iterative approach to inversion is more efficient overall. This would be facilitated by realistic and globally consistent data and age uncertainties from the paleomagnetic community.

Description: Galvanic distortion of magnetotelluric (MT) data due to small-scale surficial bodies or due to topography is one of the major factors that prevents accurate imaging of the subsurface. We present a 3-D algorithm for joint inversion of MT impedance tensor data and a frequency-independent full distortion matrix that circumvents this problem. We perform several tests of our algorithm on synthetic data affected by different amounts of distortion. These tests show that joint inversion leads to a better conductivity model compared to the inversion of the MT impedance tensor without any correction for distortion effects. For highly distorted data, inversion without any distortion correction results in strong artefacts and we cannot fit the data to the specified noise level. When the distortion is reduced, we can fit the data to an RMS of one, but still observe artefacts in the shallow part of the model. In contrast, in both cases our joint inversion can fit the data within the assumed noise level and the resulting models are comparable to the inversion of undistorted data. In addition, we show that the elements of the full distortion matrix can be well resolved by our algorithm. Finally, when inverting undistorted data, including the distortion matrix in the inversion only results in a minor loss of resolution. We therefore consider our new approach a promising tool for the general analysis of field MT data.

Description: The magnetic field that originates in the earth's core is transformed across the electrically conducting mantle before being observed, at the earth's surface or above. Assuming that the conductivity depends only on radius, it has been customary to treat the mantle as a linear time-invariant filter for the core magnetic field, with properties (as a function of the frequency ) specified by the transfer function (). An high-frequency approximation to (), which is derived from a three terms WKBJ expansion with –1/2 as small parameter, is found here to reproduce adequately, for low harmonic degrees and/or thin conducting layers, the exact solution, which is evaluated numerically. It is contrasted with the low-frequency estimation of , which consists in a perturbation procedure and in writing () as a series in powers of ( -〉 0). The low-frequency theory is applied to the magnetic variations produced by the geostrophic core flows with about 6 yr period as the phase of these flows is independently determined from their effect on the length of the day. Apart from that, the low-frequency approximation overestimates the screening by the mantle of high-frequency signals, especially the low harmonic degree ones. In practice, the attenuating factor defined from the O ( 2 ) term in the expansion of as -〉 0 cannot be retrieved from analyses of geomagnetic time-series. Application of the mantle filter theory hinges on our knowledge about the time spectrum of the magnetic field at the core surface. The low-frequency theory had been previously applied to observatory series on the assumption that geomagnetic jerks occurring in the core are rare and isolated events. Rather than following up these earlier studies, I note that the spectral density function for the second time derivative of the main magnetic field coefficients is approximately independent of in a frequency range for which the mantle has undoubtedly negligible influence. In the absence of any other information, this scaling law is extrapolated to higher frequencies.

Description: Surface nuclear magnetic resonance (surface-NMR) is a promising technique for exploring shallow subsurface aquifer structures. Surface-NMR can be applied in environments that are characterized as a 1-D layered Earth. The technique utilizes a single loop and is referred to as magnetic resonance sounding. The technique referred to as magnetic resonance tomography (MRT) allows complex 2-D aquifer structures to be explored. Currently, MRT requires multiple loops and a roll along measurement scheme, which causes long survey time. We propose a loop layout using an elongated transmitter and an in-loop receiver arrays (ETRA) to conduct a 2-D survey with just one measurement. We present a comprehensive comparison between the new layout and the common approaches based on sensitivity and resolution analyses and show synthetic and field data. The results show that ETRA generates subsurface images at sufficient resolution with significantly lower survey times than other loop layouts.

Description: We studied ionospheric responses to the 2012 April 11 M w 8.6 North Sumatra earthquake using total electron content (TEC) measurements with the regional Global Navigation Satellite System network. This earthquake ruptured the oceanic lithosphere off the Indian Ocean coast of North Sumatra, and is known as the largest strike-slip earthquake ever recorded. Coseismic ionospheric disturbances (CIDs) with rapid TEC enhancement of a few TEC units propagated northward with a speed of acoustic waves (~1 km s –1 ). Resonant atmospheric oscillation with a frequency ~4 mHz have been found as monochromatic oscillation of TEC lasting for an hour after the main shock and the largest aftershock. We compared CID amplitudes of 21 earthquakes world-wide with moment magnitudes ( M w ) 6.6–9.2. They roughly obeyed a law such that CID amplitude increases by two orders of magnitude for the M w increase of three. The 2012 North Sumatra earthquakes slightly deviated negatively from the trend possibly reflecting their strike-slip mechanisms, that is small vertical crustal movements for their magnitudes.

Description: The magnetic field of the Earth's lithosphere arises from rock magnetization contrasts that were shaped over geological times. The field can be described mathematically in spherical harmonics or with distributions of magnetization. We exploit this dual representation and assume that the lithospheric field is induced by spatially varying susceptibility values within a shell of constant thickness. By introducing a statistical assumption about the power spectrum of the susceptibility, we then derive a statistical expression for the spatial power spectrum of the crustal magnetic field for the spatial scales ranging from 60 to 2500 km. This expression depends on the mean induced magnetization, the thickness of the shell, and a power law exponent for the power spectrum of the susceptibility. We test the relevance of this form with a misfit analysis to the observational NGDC-720 lithospheric magnetic field model power spectrum. This allows us to estimate a mean global apparent induced magnetization value between 0.3 and 0.6 A m –1 , a mean magnetic crustal thickness value between 23 and 30 km, and a root mean square for the field value between 190 and 205 nT at 95 per cent. These estimates are in good agreement with independent models of the crustal magnetization and of the seismic crustal thickness. We carry out the same analysis in the continental and oceanic domains separately. We complement the misfit analyses with a Kolmogorov–Smirnov goodness-of-fit test and we conclude that the observed power spectrum can be each time a sample of the statistical one.

Description: Data on the evolution of Earth's magnetic field intensity are important for understanding the geodynamo and planetary evolution. However, the paleomagnetic record in rocks may be adversely affected by many physical processes, which must be taken into account when analysing the palaeointensity database. This is especially important in the light of an ongoing debate regarding core thermal conductivity values, and how these relate to the Precambrian geodynamo. Here, we demonstrate that several data sets in the Precambrian palaeointensity database overestimate the true paleofield strength due to the presence of non-ideal carriers of palaeointensity signals and/or viscous re-magnetizations. When the palaeointensity overestimates are removed, the Precambrian database does not indicate a robust change in geomagnetic field intensity during the Mesoproterozoic. These findings call into question the recent claim that the solid inner core formed in the Mesoproterozoic, hence constraining the thermal conductivity in the core to ‘moderate’ values. Instead, our analyses indicate that the presently available palaeointensity data are insufficient in number and quality to constrain the timing of solid inner core formation, or the outstanding problem of core thermal conductivity. Very young or very old inner core ages (and attendant high or low core thermal conductivity values) are consistent with the presently known history of Earth's field strength. More promising available data sets that reflect long-term core structure are geomagnetic reversal rate and field morphology. The latter suggests changes that may reflect differences in Archean to Proterozoic core stratification, whereas the former suggest an interval of geodynamo hyperactivity at ca. 550 Ma.

Description: Red Clay underlying the loess-palaeosol sequences on the Chinese Loess Plateau is an eolian deposit. There is a controversy over whether magnetic susceptibility ( ) variations in Red Clay sequence can be used as an indicator of summer palaeomonsoon intensity. This study investigates the magnetic mineralogy, magnetic concentration and magnetic grain size distribution of Jiaxian Red Clay with multimagnetic methods. Our results indicate that the magnetic properties of Jiaxian Red Clay are similar to those of the Quaternary loess-palaeosol sequences, and ultrafine ferrimagnetic grains produced during pedogenesis are responsible for an increase in susceptibility, therefore the enhancement mechanism of Red Clay is similar to that of the overlying loess-palaeosol sequences. This paper explores variations in the Red Clay sequence through spatial and temporal analysis. The susceptibility variation of six sites along a NNE to SSW transect correlate to palaeoclimatic cycles, so can be used to trace the summer palaeomonsoon intensity from a spatial perspective. However, a simple loess-derived calibration function cannot be used to quantitative reconstruct the palaeomonsoon intensity variations thought time. An adjusted calibration function for palaeosols from Red Clay sequence needs to be developed, so that can be used to quantitative reconstruct palaeomonsoon intensity. Further study is necessary to develop such a transfer function.

Description: Palaeointensity experiments were carried out to a sample collection from two sections of basalt lava flow sequences of Pliocene age in north central Iceland (Chron C2An) to further refine the knowledge of the behaviour of the palaeomagnetic field. Selection of samples was mainly based on their stability of remanence to thermal demagnetization as well as good reversibility in variations of magnetic susceptibility and saturation magnetization with temperature, which would indicate the presence of magnetite as a product of deuteric oxidation of titanomagnetite. Among 167 lava flows from two sections, 44 flows were selected for the Königsberger–Thellier–Thellier experiment in vacuum. In spite of careful pre-selection of samples, an Arai plot with two linear segments, or a concave-up appearance, was often encountered during the experiments. This non-ideal behaviour was probably caused by an irreversible change in the domain state of the magnetic grains of the pseudo-single-domain (PSD) range. This is assumed because an ideal linear plot was obtained in the second run of the palaeointensity experiment in which a laboratory thermoremanence acquired after the final step of the first run was used as a natural remanence. This experiment was conducted on six selected samples, and no clear difference between the magnetic grains of the experimented and pristine sister samples was found by scanning electron microscope and hysteresis measurements, that is, no occurrence of notable chemical/mineralogical alteration, suggesting that no change in the grain size distribution had occurred. Hence, the two-segment Arai plot was not caused by the reversible multidomain/PSD effect in which the curvature of the Arai plot is dependent on the grain size. Considering that the irreversible change in domain state must have affected data points at not only high temperatures but also low temperatures, f v ≥ 0.5 was adopted as one of the acceptance criteria where f v is a vectorially defined fraction of the linear segment. A measure of curvature k ' was also used to check the linearity of the selected linear segment. It was avoided, however, to reject the result out of hand by the large curvature k of the entire data points because it might still include a linear segment with a large fraction. Combining with the results of Shaw's experiments, 52 palaeointensities were obtained out of 192 specimens, or 11 flow means were obtained out of the 44 lava flows. Most of the palaeointensities were from the upper part of the lava section (Chron C2An.1n) and ranged between 30 and 66 μT. Including two results from the bottom part of the lava section, the mean virtual dipole moment for 2.5–3.5 Ma is 6.3 ± 1.4  x  10 22 Am 2 ( N  = 11), which is ~19 per cent smaller than the present-day dipole moment.

Description: The sedimentary sequence deposited during the deglaciation phase following the last glacial maximum in the Storfjorden trough, on the northwestern Barents Sea south of Svalbard, was sampled with 10 piston and gravity cores during the SVAIS and EGLACOM cruises. Three cores (SV-02, SV-03 and SV-05) collected on the upper continental slope are characterized by a thin (20–40 cm) Holocene interval and a thick (up to 4.5 m in core SV-03) late Pleistocene sequence of finely laminated fine-grained sediments that have been interpreted as plumites deposited during the Melt Water Pulse 1a (MWP-1a). Radiocarbon ages obtained at the top and bottom of this stratigraphic interval revealed that deposition occurred during less than two centuries at around 15 ka ago, with a very high sedimentary rate exceeding 3 cm a –1 . We studied the palaeomagnetic and rock magnetic properties of this interval, by taking magnetic measurements at 1 cm spacing on u-channel samples collected from the three cores. The data show that this sequence is characterized by good palaeomagnetic properties and the palaeomagnetic and rock magnetic trends may be correlated at high resolution from core to core. The obtained palaeomagnetic data therefore offer the unique opportunity to investigate in detail the rate of geomagnetic palaeosecular variation (PSV) in the high northern latitudes at a decadal scale. Notwithstanding the palaeomagnetic trends of the three cores may be closely matched, the amplitude of directional PSV and the consequent virtual geomagnetic pole (VGP) scatter (S) is distinctly higher in one core (SV-05) than in the other two cores (SV-02 and SV-03). This might result from a variable proportion of two distinct populations of magnetic minerals in core SV-05, as suggested by the variable tendency to acquire a gyromagnetic remanent magnetization at high fields during the AF demagnetization treatment. For the plumite interval of cores SV-02 and SV-03, where the magnetic mineralogy is uniform and magnetite is the main magnetic carrier, a S value of about 9° is obtained. We consider this value as a reliable approximation of palaeomagnetic secular variation at a latitude of 75°N over a time interval spanning a couple of centuries around 15 ka ago.

Description: A recent study of the Matuyama–Brunhes (M-B) geomagnetic field reversal recorded in exposed lacustrine sediments from the Sulmona Basin (Italy) provided a continuous, high-resolution record indicating that the reversal of the field direction at the terminus of the M-B boundary (MBB) occurred in less than a century, about 786 ka ago. In the sediment, thin (4–6 cm) remagnetized horizons were recognized above two distinct tephra layers—SUL2-19 and SUL2-20—that occur ~25 and ~35 cm below the MBB, respectively. Also, a faint, millimetre-thick tephra (SUL2-18) occurs 2–3 cm above the MBB. With the aim of improving the temporal resolution of the previous Sulmona MBB record and understanding the possible influence of cryptotephra on the M-B record in the Sulmona Basin, we performed more detailed sampling and analyses of overlapping standard and smaller samples from a 50 cm-long block that spans the MBB. The new data indicate that (i) the MBB is even sharper than previously reported and occurs ~2.5 cm below tephra SUL2-18, in agreement with the previous study; (ii) the MBB coincides with the rise of an intensity peak of the natural remanent magnetization (NRM) intensity, which extends across SUL2-18; (iii) except for a 2-cm-thick interval just above tephra SUL2-18, the rock magnetic parameters ( k , ARM, M r , M s , B c , B cr ) indicate exactly the same magnetic mineralogy throughout the sampled sequence. We conclude that either SUL2-18 resulted in the remagnetization of an interval of about 6 cm (i.e. during the NRM intensity peak spanning ~260 ± 110 yr, according to the estimated local sedimentation rate), and thus the detailed MBB record is lost because it is overprinted, or the MBB is well recorded, occurred abruptly about 2.5 cm below SUL2-18 and lasted less than 13 ± 6 yr. Both hypotheses challenge our understanding of the geomagnetic field behaviour during a polarity transition and/or of the NRM acquisition process in the Sulmona lacustrine sediment.

Description: We test the ability of large-scale velocity fields inferred from geomagnetic secular variation data to produce the global magnetic field of the Earth. Our kinematic dynamo calculations use quasi-geostrophic (QG) flows inverted from geomagnetic field models which, as such, incorporate flow structures that are Earth-like and may be important for the geodynamo. Furthermore, the QG hypothesis allows straightforward prolongation of the flow from the core surface to the bulk. As expected from previous studies, we check that a simple QG flow is not able to sustain the magnetic field against ohmic decay. Additional complexity is then introduced in the flow, inspired by the action of the Lorentz force. Indeed, on centennial timescales, the Lorentz force can balance the Coriolis force and strict quasi-geostrophy may not be the best ansatz. When our columnar flow is modified to account for the action of the Lorentz force, magnetic field is generated for Elsasser numbers larger than 0.25 and magnetic Reynolds numbers larger than 100. This suggests that our large-scale flow captures the relevant features for the generation of the Earth's magnetic field and that the invisible small-scale flow may not be directly involved in this process. Near the threshold, the resulting magnetic field is dominated by an axial dipole, with some reversed flux patches. Time dependence is also considered, derived from principal component analysis applied to the inverted flows. We find that time periods from 120 to 50 yr do not affect the mean growth rate of the kinematic dynamos. Finally, we note that the footprint of the inner core in the magnetic field generated deep in the bulk of the shell, although we did not include one in our computations.

Description: The Auca Mahuida volcano (2.03–0.88 Ma) located east of the Andean thrust front in the Neuquén basin (Argentina) hosts an oil system of thermogenic origin and is affected by the NW–SE striking-faults. Intrusive bodies and the underlying Jurassic sediments constitute the reservoir rocks. Aeromagnetic data collected in the Auca Mahuida area detected multiple dipolar magnetic anomalies, many of which have reverse polarity. Palaeomagnetic measurements on rock samples collected in the field together with available age determinations indicate that the reversely magnetized sources were mainly emplaced during the Matuyama reverse polarity chron while the normal polarity sources were emplaced during the Olduvai and/or Jaramillo subchrons. The location and geometry of the intrusive bodies is poorly known and the customary magnetic inversion is rendered difficult because of multiple natural remanent magnetization directions. To address these difficulties, a voxel inversion was applied to model the vector residual magnetic intensity (VRMI) transformation of the observed total magnetic intensity data. The modelling showed a 1.5 km deep, subcircular ring-shaped intrusion below the summit of the volcano and a series of NW–SE elongated, fault-controlled intrusive bodies to depths up to 3–4 km. Our results show that magnetic data and VRMI modelling help resolve the geometry of the shallow plumbing system of volcanoes with remanently magnetized sources, and estimate the depth and geometry of potential oil reservoirs in volcanic areas.

Description: Was the Yukon-Tanana Terrane (YTT), a California-sized part of south-central Yukon, an autochthonous or para-autochthonous part of northern British Columbia in the Early Cretaceous or was it part of a proposed allochthonous ‘Baja B.C.’ continent offshore of southern California? To answer this fundamental question, a paleomagnetic study has been completed on 347 specimens from 24 sites in the 114.7 ± 1.1 Ma Quiet Lake batholith. This 1300 km 2 pluton is composed mostly of massive medium-to-coarse grained biotite quartz monzonite that exhibits no evidence of either deformation or metamorphism, and that intrudes metamorphosed pre-Cretaceous basement rocks of the YTT in southern Yukon. The paleomagnetic analysis utilized thermal and alternating field step demagnetization, and saturation isothermal remanence methods. A well-defined characteristic remanent magnetization (ChRM) direction was isolated throughout the 500–585 °C temperature range at Decl. = 340.6°, Incl. = 77.4° ( N = 14 sites, k = 51.2, A 95 = 5.6°). The ChRM resides in magnetite with a low titanium content and is interpreted to be a primary thermoremanent magnetization. After correction for 490 km of geologically demonstrable dextral displacement on the inboard Tintina fault zone, the Quiet Lake batholith's paleopole is not significantly different at 95 per cent confidence from the co-eval 115 Ma reference paleopole for North America, giving non-significant translation and rotation estimates of 1.4° ± 5.1° (1) northwestwards and 10° ± 13° (1) clockwise, respectively. Thus, this is the first Early Cretaceous paleopole to show clearly that the YTT in Yukon is a para-autochthon that was part of North America's continental margin at that time. Further, after correction for Tintina fault displacement, the eight available Mesozoic YTT paleopoles agree closely with the North American apparent polar wander path (APWP). In contrast, the 22 paleopoles from the Intermontane Belt show the expected behaviour of an allochthonous thin-skin of terranes (i.e. ‘Baja B.C.’ behaviour) and record ~8° of northward translation and ~50° of clockwise rotation relative to both the YTT and North American APWPs during the Mesozoic and Paleogene.

Description: There exists a fundamental as well as practical interest in being able to accurately forecast the future evolution of Earth's magnetic field at decadal to secular ranges. This work enables such forecasts by combining geomagnetic data with an Earth-like numerical model of a convection-driven fluid dynamo. The underlying data assimilation framework builds on recent progress in inverse geodynamo modelling, a method which estimates an internal dynamic structure for Earth's core from a snapshot of the magnetic field and its instantaneous rate of change at the surface, and takes advantage of linear relationships and long-range correlations between observed and hidden state variables. Here the method is further evolved into a single-epoch ensemble Kalman filter, in order to initialise at a given epoch an ensemble of states compatible with the observations and representative of the uncertainties in the estimation of hidden quantities. The ensemble dynamics, obtained by subsequent numerical integration of the prognostic model equations, are found to be governed by a thermal wind balance or equilibrium between buoyancy forces, the Coriolis force and the pressure gradient. The resulting core fluid flow pattern is a quasi-steady eccentric gyre organised in a column parallel to Earth's rotation axis, in equilibrium with a longitudinal hemispheric convective density anomaly pattern. The flow provides induction for the magnetic field, which also undergoes a realistic amount of diffusion. Predictions of the present magnetic field from data taken within the past century show that the ensemble has an average retaining good consistency with the true geomagnetic evolution and an acceptable spread well representative of prediction errors, up to at least a secular range. The predictability of the geodynamo thus appears to significantly exceed previous theoretical expectations based on the chaotic divergence of ensemble members. The assimilation generally outperforms the linear mathematical extrapolations from a 30-yr prediction range onwards, with a 40 per cent improvement in Earth-surface error at a secular range. The geomagnetic axial dipole decay observed over the past two centuries is predicted to continue at a similar pace in the next century, with a further loss of $1.1\pm 0.3 \hspace{2.84526pt}\mathrm{\mu T}$ by year 2115. The focal (or minimum intensity) point of the South Atlantic geomagnetic anomaly is predicted to enter the South Pacific region in the next century, with the anomaly itself further deepening and widening. By year 2065, the minimum intensity is predicted to decrease by $1.46\pm 0.4 \hspace{2.84526pt}\mathrm{\mu T}$ at the Earth surface and the focal point to move 12.8 ± 1.4 deg westwards with a slight northward component. This corresponds to a drift rate of 0.26 deg yr –1 , similar to the typical geomagnetic westward drift observed over the past four centuries. The same drift rate is also predicted until 2115 with a further (but more uncertain) intensity decrease.

Description: Understanding the contribution of biogenic magnetic particles into sedimentary assemblages is a current challenge in palaeomagnetism. It has been demonstrated recently that magnetic particles produced through biologically controlled mineralization processes, such as magnetosomes from magnetotactic bacteria, contribute to the recording of natural remanent magnetization in marine and lacustrian sediments. Contributions from other, biologically induced, mineralization types, which are known from multiple laboratory experiments to include magnetic minerals, remain largely unknown. Here, we report magnetic properties of iron minerals formed by a community of iron- and manganese-reducing bacteria isolated from a natural groundwater deposit during a 2 yr long incubation experiment. The main iron phases of the biomineralized mass are lepidocrocite, goethite and magnetite, each of which has environmental significance. Unlike the majority of the previous studies that reported superparamagnetic grain size, and thus no remanence carrying capacity of biologically induced magnetite, hysteresis and first-order reversal curves measurements in our study have not detected significant superparamagnetic contribution. The biomineralized mass, instead, contains a mixture of single-domain to pseudo-single-domain and multidomain magnetite particles that are capable of carrying a stable chemical remanent magnetization. Isothermal remanent magnetization acquisition parameters and first-order reversal curves signatures of the biomineralized samples deviate from previously proposed criteria for the distinction of extracellular (biologically induced) magnetic particles in mixtures. Given its potential significance as a carrier of natural remanent magnetization, environmental requirements, distribution in nature and the efficiency in the geomagnetic field recording by biologically induced mineralization need comprehensive investigation.

Description: Monthly means of the magnetic field measurements at ground observatories are a key data source for studying temporal changes of the core magnetic field. However, when they are calculated in the usual way, contributions of external (magnetospheric and ionospheric) origin may remain, which make them less favourable for studying the field generated by dynamo action in the core. We remove external field predictions, including a new way of characterizing the magnetospheric ring current, from the data and then calculate revised monthly means using robust methods. The geomagnetic secular variation (SV) is calculated as the first annual differences of these monthly means, which also removes the static crustal field. SV time-series based on revised monthly means are much less scattered than those calculated from ordinary monthly means, and their variances and correlations between components are smaller. On the annual to decadal timescale, the SV is generated primarily by advection in the fluid outer core. We demonstrate the utility of the revised monthly means by calculating models of the core surface advective flow between 1997 and 2013 directly from the SV data. One set of models assumes flow that is constant over three months; such models exhibit large and rapid temporal variations. For models of this type, less complex flows achieve the same fit to the SV derived from revised monthly means than those from ordinary monthly means. However, those obtained from ordinary monthly means are able to follow excursions in SV that are likely to be external field contamination rather than core signals. Having established that we can find models that fit the data adequately, we then assess how much temporal variability is required. Previous studies have suggested that the flow is consistent with torsional oscillations (TO), solid body-like oscillations of fluid on concentric cylinders with axes aligned along the Earth's rotation axis. TO have been proposed to explain decadal timescale changes in the length-of-day. We invert for flow models where the only temporal changes are consistent with TO, but such models have an unacceptably large data misfit. However, if we relax the TO constraint to allow a little more temporal variability, we can fit the data as well as with flows assumed constant over three months, demonstrating that rapid SV changes can be reproduced by rather small flow changes. Although the flow itself changes slowly, its time derivative can be locally (temporally and spatially) large, in particular when and where core surface secular acceleration peaks. Spherical harmonic expansion coefficients of the flows are not well resolved, and many of them are strongly correlated. Averaging functions, a measure of our ability to determine the flow at a given location from the data distribution available, are poor approximations to the ideal, even when centred on points of the core surface below areas of high observatory density. Both resolution and averaging functions are noticeably worse for the toroidal flow component, which dominates the flow, than the poloidal flow component, except around the magnetic equator where averaging functions for both components are poor.

Description: The electrical rock conductivity is a sensitive indicator for carbon dioxide (CO 2 ) injection and migration processes. For a reliable balancing of the free CO 2 in pore space with petrophysical models such as Archie's law or for the detection of migrating CO 2 , detailed knowledge of the pore water conductivity during interaction with CO 2 is essential but not available yet. Contrary to common assumptions, pore water conductivity cannot be assumed constant since CO 2 is a reactive gas that dissolves into the pore water in large amounts and provides additional charge carriers due to the dissociation of carbonic acid. We consequently carried out systematic laboratory experiments to quantify and analyse the changes in saline pore water conductivity caused by CO 2 at thermodynamic equilibrium. Electrical conductivity is measured on pore water samples for pressures up to 30 MPa and temperatures up to 80 °C. The parameter range covers the gaseous, liquid and supercritical state of the CO 2 involved. Pore water salinities from 0.006 up to 57.27 g L –1 sodium chloride were investigated as well as selective other ion species. At the same time, the CO 2 concentration in the salt solution was determined by a wet-chemical procedure. A two-regime behaviour appears: for small salinities, we observe an increase of up to more than factor 3 in the electrical pore water conductivity, which strongly depends on the solution salinity (low-salinity regime). This is an expected behaviour, since the additional ions originating from the dissociation of carbonic acid positively contribute to the solution conductivity. However, when increasing salinities are considered this effect is completely diminished. For highly saline solutions, the increased mutual impeding causes the mobility of all ions to decrease, which may result in a significant reduction of conductivity by up to 15 per cent despite the added CO 2 (high-salinity regime). We present the data set covering the pressure, temperature, salinity and ion species dependence of the CO 2 effect. Furthermore, the observations are analysed and predicted with a semi-analytical formulation for the electrical pore water conductivity taking into account the species’ interactions. For the applicability of our results in practice of exploration and monitoring, we additionally provide a purely empirical formulation to compute the impact of CO 2 on pore water conductivity at equilibrium which only requires the input of pressure, temperature and salinity information.

Description: Although the Matuyama–Brunhes boundary (MBB) in the Chinese Loess Plateau (CLP) is very important in reliably correlating Quaternary loess with other sediments in the world, particularly with marine and polar ice cores, its exact stratigraphic position remains controversial. Previous investigations usually placed the MBB between paleosol unit S8 and loess unit L8 in various locations. To better understand the spatial differences in the MBB position, a high-resolution paleomagnetic study was conducted in a loess section of the Lantian Basin at the southern margin of the CLP. The results show that the MBB is situated in the middle of the relatively weak paleosol unit S7, consistent with a recent report on the MBB based on a 10 Be study from the Xifeng and Luochuan loess sections of the central CLP. However, the regional anomalously low magnetic susceptibility in paleosols S7 and S8 indicates that it is more reliable to determine the paleoclimate boundaries between loess and paleosol horizons of this segment with median grain size. Then, the MBB in the Yushan section can be correlated with the bottom of paleosol S7, corresponding to the older part of interglacial marine isotope stage 19. This result temporally reconciles the striking discrepancy of the position of the MBB recorded in between loess and other typical sedimentary sequences, and further confirms that the stratigraphic position of the MBB could spatially vary to a certain extent due to regional sedimentary or paleoclimatic conditions in the marginal areas of the CLP. In the Yushan section, the high-frequency variations of paleomagnetic directions during a long period of ~31 ka before the MBB, however, could not be attributed to a genuine response to the true geomagnetic behaviour. Moreover, the climate offset defined by the magnetic susceptibility and median grain size of the section can be preliminarily attributed to the regional geology and paleoenvironment background. A multiproxy-based stratigraphic division is considered very necessary when paleomagnetic and climatic boundaries are defined exactly in a specific area of the southern CLP.

Description: We study rotating thermal convection in spherical shells as prototype for flow in the cores of terrestrial planets, gas planets or in stars. We base our analysis on a set of about 450 direct numerical simulations of the (magneto)hydrodynamic equations under the Boussinesq approximation. The Ekman number ranges from 10 –3 to 10 –5 . The supercriticality of the convection reaches about 1000 in some models. Four sets of simulations are considered: non-magnetic simulations and dynamo simulations with either free-slip or no-slip flow boundary conditions. The non-magnetic setup with free-slip boundaries generates the strongest zonal flows. Both non-magnetic simulations with no-slip flow boundary conditions and self-consistent dynamos with free-slip boundaries have drastically reduced zonal-flows. Suppression of shear leads to a substantial gain in heat-transfer efficiency, increasing by a factor of 3 in some cases. Such efficiency enhancement occurs as long as the convection is significantly influenced by rotation. At higher convective driving the heat-transfer efficiency tends towards that of the classical non-rotating Rayleigh–Bénard system. Analysis of the latitudinal distribution of heat flow at the outer boundary reveals that the shear is most effective at suppressing heat-transfer in the equatorial regions. Simulations with convection zones of different thickness show that the zonal flows become less energetic in thicker shells, and, therefore, their effect on heat-transfer efficiency decreases. Furthermore, we explore the influence of the magnetic field on the non-zonal flow components of the convection. For this we compare the heat-transfer efficiency of no-slip non-magnetic cases with that of the no-slip dynamo simulations. We find that at E = 10 –5 magnetic field significantly affects the convection and a maximum gain of about 30 per cent (as compared to the non-magnetic case) in heat-transfer efficiency is obtained for an Elsasser number of about 3. Our analysis motivates us to speculate that convection in the polar regions in dynamos at E = 10 –5 is probably in a ‘magnetostrophic’ regime.

Description: ‘Anomalous’ magnetic fabrics in dikes that appear to indicate flow into the wall confound many workers. Here, we present extensive magnetic data on five dikes from Tenerife, Canary Islands, and use these to interpret the causes of the anomalous fabrics. Comparison of the anisotropy of magnetic susceptibility (AMS) and anhysteretic magnetization (AARM) results show that, in some cases, the anomalous fabrics are caused by single-domain grains, which produce AMS fabrics perpendicular to the grain elongation, whereas AARM fabrics are parallel. To check this, hysteresis experiments were used to characterize the domain state. These show most are mixtures of pseudo-single-domain or single-domain plus multi-domain particles, but many have wasp-waisted hysteresis loops, likely indicating mixed populations of stable single-domain and superparamagnetic grains. First-order reversal curves were used to better characterize this and show mixtures of stable single-domain and superparamagnetic grains dominate the magnetic signal. Magnetic particles at the stable single-domain/superparamagnetic threshold are unstable at timespans relevant to the analytical techniques, so they produce complicated results. This suggests that anomalous AMS fabrics in dikes cannot simply be attributed to elongated stable single-domain particles and that mixtures of the different grain types can produce hybrid fabrics, in which the fabrics are neither perpendicular or parallel to the dike plane, that are difficult to interpret without extensive magnetic analysis.

Description: We derive scaling relationships for planetary dynamos based on a balance between energy production and Joule dissipation, and between the curl of the buoyancy and Coriolis forces. These scaling relationships are deduced for the particular case of dynamos driven by helical waves, but are shown to have a much broader applicability. They are consistent with the evidence of the numerical dynamos, yielding predictions consistent with published empirical scaling laws and also with the observed transition from dipolar to multipolar dynamos. A direct comparison with the observational evidence for the planets is hampered by the fact that we do not know what sets the smallest scale of the motion in the planets. Nevertheless, we use our scaling relationships to show that the traditional assumption that the Elsasser number is of order unity is inconsistent with the observation that the gas-giant dynamos are dipolar dynamos, as is the more recent suggestion that the strength of the dipole is independent of rotation rate and controlled by the buoyancy flux alone. On the other hand, we show that the observational data is consistent with the hypothesis that a dipolar dynamo saturates at the lowest permissible magnetic energy compatible with a given buoyancy flux.

Description: We present a detailed palaeomagnetic study from 35 sites on Holocene lava flows of the Tongariro Volcanic Centre, central North Island, New Zealand. Prior to the study the eruption ages of these flows were constrained to within a few thousand years by recently published high-precision 40 Ar/ 39 Ar geochronological data and tephrostratigraphic controls. Correlation of flow mean palaeomagnetic directions with a recently published continuous sediment record from Lake Mavora, Fiordland, allows us to reduce the age uncertainty to 300–500 yr in some cases. Our refined ages significantly improve the chronology of Holocene effusive eruptions of the volcanoes of the Tongariro Volcanic Centre. For instance, differences in the palaeomagnetic directions recorded by lavas from the voluminous Iwikau and Rangataua members suggest that individual effusive periods lasted up to thousands of years and that these bursts have been irregularly spaced over time. While over the last few millennia the effusive eruptive activity from Mt Ruapehu has been relatively quiet, the very young age (200–500 BP) of a Red Crater sourced flow suggests that effusive activity around Mt Tongariro lasted into the past few centuries. This adds an important hazard context to the historical record, which has otherwise comprised frequent relatively small, tephra producing, explosive eruptions without the production of lava flows.

Description: The inversion of magnetotelluric data into subsurface electrical conductivity poses an ill-posed problem. Smoothing constraints are widely employed to estimate a regularized solution. Here, we present an alternative inversion scheme that estimates a sparse representation of the model in a wavelet basis. The objective of the inversion is to determine the few non-zero wavelet coefficients which are required to fit the data. This approach falls into the class of sparsity constrained inversion schemes and minimizes the combination of the data misfit in a least-squares 2 sense and of a model coefficient norm in an 1 sense ( 2 - 1 minimization). The 1 coefficient norm renders the solution sparse in a suitable representation such as the multiresolution wavelet basis, but does not impose explicit structural penalties on the model as it is the case for 2 regularization. The presented numerical algorithm solves the mixed 2 - 1 norm minimization problem for the nonlinear magnetotelluric inverse problem. We demonstrate the feasibility of our algorithm on synthetic 2-D MT data as well as on a real data example. We found that sparse models can be estimated by inversion and that the spatial distribution of non-vanishing coefficients indicates regions in the model which are resolved.

Description: We examined the mathematical and computational aspects of the magnetic potential, vector and gradient tensor fields of a tesseroid in a geocentric spherical coordinate system (SCS). This work is relevant for 3-D modelling that is performed with lithospheric vertical scales and global, continent or large regional horizontal scales. The curvature of the Earth is significant at these scales and hence, a SCS is more appropriate than the usual Cartesian coordinate system (CCS). The 3-D arrays of spherical prisms (SP; ‘tesseroids’) can be used to model the response of volumes with variable magnetic properties. Analytical solutions do not exist for these model elements and numerical or mixed numerical and analytical solutions must be employed. We compared various methods for calculating the response in terms of accuracy and computational efficiency. The methods were (1) the spherical coordinate magnetic dipole method (MD), (2) variants of the 3-D Gauss–Legendre quadrature integration method (3-D GLQI) with (i) different numbers of nodes in each of the three directions, and (ii) models where we subdivided each SP into a number of smaller tesseroid volume elements, (3) a procedure that we term revised Gauss–Legendre quadrature integration (3-D RGLQI) where the magnetization direction which is constant in a SCS is assumed to be constant in a CCS and equal to the direction at the geometric centre of each tesseroid, (4) the Taylor's series expansion method (TSE) and (5) the rectangular prism method (RP). In any realistic application, both the accuracy and the computational efficiency factors must be considered to determine the optimum approach to employ. In all instances, accuracy improves with increasing distance from the source. It is higher in the percentage terms for potential than the vector or tensor response. The tensor errors are the largest, but they decrease more quickly with distance from the source. In our comparisons of relative computational efficiency, we found that the magnetic potential takes less time to compute than the vector response, which in turn takes less time to compute than the tensor gradient response. The MD method takes less time to compute than either the TSE or RP methods. The efficiency of the (GLQI and) RGLQI methods depends on the polynomial order, but the response typically takes longer to compute than it does for the other methods. The optimum method is a complex function of the desired accuracy, the size of the volume elements, the element latitude and the distance between the source and the observation. For a model of global extent with typical model element size (e.g. 1 degree horizontally and 10 km radially) and observations at altitudes of 10s to 100s of km, a mixture of methods based on the horizontal separation of the source and observation separation would be the optimum approach. To demonstrate the RGLQI method described within this paper, we applied it to the computation of the response for a global magnetization model for observations at 300 and 30 km altitude.

Description: Kamacite particles (Fe–Ni, Ni 〈 5 per cent), are very common in extra-terrestrial materials, such as meteorites. It is normally assumed that for kamacite particles to be reliable recorders of magnetic fields, they need to be magnetically uniform (single domain, SD) and thermally stable. Larger particles subdivide into non-uniform multidomain (MD) magnetic structures that produce weaker magnetic signals, while small SD particles become magnetically unstable due to thermal fluctuations and exhibit superparamagnetic behaviour. In this paper we determine the first micromagnetic calculation of the stable SD range domain-state phase diagram for metallic iron; previous calculations were analytical. There is a significant increase in the critical size for the SD/MD threshold size, for example, for cube-shaped iron particles, the critical SD/MD threshold has now been estimated to be 25 nm, compared to 17 nm for previous estimates. The larger critical SD/MD threshold size for iron, agrees better with previously published nanometric observations of domain state for FeNi particles, then early analytical models.

Description: This paper presents a distributed magnetotelluric inversion scheme based on adaptive finite-element method (FEM). The key novel aspect of the introduced algorithm is the use of automatic mesh refinement techniques for both forward and inverse modelling. These techniques alleviate tedious and subjective procedure of choosing a suitable model parametrization. To avoid overparametrization, meshes for forward and inverse problems were decoupled. For calculation of accurate electromagnetic (EM) responses, automatic mesh refinement algorithm based on a goal-oriented error estimator has been adopted. For further efficiency gain, EM fields for each frequency were calculated using independent meshes in order to account for substantially different spatial behaviour of the fields over a wide range of frequencies. An automatic approach for efficient initial mesh design in inverse problems based on linearized model resolution matrix was developed. To make this algorithm suitable for large-scale problems, it was proposed to use a low-rank approximation of the linearized model resolution matrix. In order to fill a gap between initial and true model complexities and resolve emerging 3-D structures better, an algorithm for adaptive inverse mesh refinement was derived. Within this algorithm, spatial variations of the imaged parameter are calculated and mesh is refined in the neighborhoods of points with the largest variations. A series of numerical tests were performed to demonstrate the utility of the presented algorithms. Adaptive mesh refinement based on the model resolution estimates provides an efficient tool to derive initial meshes which account for arbitrary survey layouts, data types, frequency content and measurement uncertainties. Furthermore, the algorithm is capable to deliver meshes suitable to resolve features on multiple scales while keeping number of unknowns low. However, such meshes exhibit dependency on an initial model guess. Additionally, it is demonstrated that the adaptive mesh refinement can be particularly efficient in resolving complex shapes. The implemented inversion scheme was able to resolve a hemisphere object with sufficient resolution starting from a coarse discretization and refining mesh adaptively in a fully automatic process. The code is able to harness the computational power of modern distributed platforms and is shown to work with models consisting of millions of degrees of freedom. Significant computational savings were achieved by using locally refined decoupled meshes.

Description: A spatially averaged mean-field model for fully or partially ordered members of the ilmenite–hematite solid solution series is rigorously derived from the Heisenberg Hamiltonian by first assuming no temporal correlation of atomic spins, and then by spatially averaging over spins at equivalent atomic positions. The model is based on the geometry of exchange interactions between nearest and next-nearest neighbours and predicts magnetization curves in homogenous solid solutions with variable degree of order. While the general framework presented can also be applied to atomic scale models, and to other solid solution series, here the symmetries of the ilmenite–hematite lattice are exploited to show that four different sublattice magnetizations and six independent combinations of exchange constants determine the temperature variation of the magnetization curves. Comparing measured Curie temperatures T C and M s ( T ) curves to model predictions results in accurate constraints for these combinations. It is also possible to calculate predictions for high-field magnetization slopes HF , which not only improve accurate experimental determination of the Curie temperature but also provide a new magnetic method to estimate the order parameter for ilmenite–hematite solid solution samples.

Description: Tibetan Plateau uplift and global cooling have generally been thought to have caused the drying of the Asian inland, but how and when these factors drove the aridification is unknown. The Xining Basin at the NE Tibetan Plateau received continuous Eocene–Miocene fine-grained sediments, providing an excellent opportunity to address this question. Here we present detailed rock magnetic and diffuse reflectance spectroscopy (DRS) analyses for a well-dated Cenozoic sedimentary record from the Xiejia section in the basin. Magnetic susceptibility ( ), saturation magnetization (Ms) and saturation isothermal remanent magnetization (SIRM) in this section show a long-term decreasing trend from ~52 to ~25 Ma, well coinciding with global cooling and drying in the region, and an increasing trend since ~25 Ma, which is in contrast to the further progressing aridification of the basin. Thermomagnetic results and DRS-determined hematite contents suggest that the relative content of magnetite and hematite is the main control on the , Ms and SIRM values. We argue that the long-term Eocene–Oligocene global cooling increased the drying of the Asian inland, lowering the lake level and exposing a larger area to low temperature oxidation for longer times, thus producing more hematite and leading to the decreasing trend of magnetic concentration parameters from ~52 to ~25 Ma. An intensive uplift of the NE Tibetan Plateau since ~25 Ma, associated with a change in the sedimentary source, might be responsible for the increase of , Ms and SIRM after 25 Ma.

Description: The Mt Amram igneous complex (AIC) represents northern tip of the Neoproterozoic Arabian Nubian Shield (ANS). For the first time the AIC deep structure was studied using the gravity, aero and ground magnetic, magnetic susceptibility and density measurements and geological data. Analysing all available data at the Amram area we concluded what only monzonite body can be reason for gravity high and coinciding reduced to pole (RTP) maximum. Geological knowledge allowed suggesting its intrusive character and compact body form. Cluster of inverse solutions (Werner deconvolution) localized this body as initial model for forward modelling. Further iterations (2 3/4 -D forward modelling) clarified the monzonite geometry and properties; the modelling allowed also to investigate the non-uniqueness and estimate also the confident intervals for final solution. The research consists three interconnected stages. At the detailed scale, ground magnetic data suggested three magmatic blocks of few hundred meters shifted dextral about 100 m along the Zefunut fault. Estimated accuracy for geometry of the magnetic bodies is a few tens metres. At the middle scale, quantitative gravity and magnetic interpretations provide model of the monzonite body, which is an order of magnitude more than the volume of the felsic rhyolites and granite rocks. Boundary of the whole monzonite body was estimated with accuracy as a hundred meters. As a result we suggest that the parent magma for the AIC is the monzonite, similar to the model suggested for the Timna Igneous Complex 12 km north of the AIC. The model developed can be applied to evaluate the subsurface volumes of the mafic magmatic rocks in adjacent locations. At the regional scale for exposed the Sinai and Arab Saudi Precambrian crystalline shield our approach allows to understand the apparent contradiction between geological predominantly granite composition (low magnetic rocks) and magnetic data. The aeromagnetic data show number strong magnetic anomalies suggesting the presence large volume of high magnetic (mainly basic) rocks at the depth. This problem is proposed for future research.

Description: 3-D electromagnetic (EM) studies of the Earth have advanced significantly over the past decade. Despite a certain success of the 3-D EM inversions of real data sets, the quantitative assessment of the recovered models is still a challenging problem. It is known that one can gain valuable information about model uncertainties from the analysis of Hessian matrix. However, even with modern computational capabilities the calculation of the Hessian matrix based on numerical differentiation is extremely time consuming. Much more efficient way to compute the Hessian matrix is provided by an ‘adjoint sources’ methodology. The computation of Hessian matrix (and Hessian-vector products) using adjoint formulation is now well-established approach, especially in seismic inverse modelling. As for EM inverse modelling we did not find in the literature a description of the approach, which would allow EM researchers to apply this methodology in a straightforward manner to their scenario of interest. In the paper, we present formalism for the efficient calculation of the Hessian matrix using adjoint sources approach. We also show how this technique can be implemented to calculate multiple Hessian-vector products very efficiently. The formalism is general in the sense that it allows to work with responses that arise in EM problem set-ups either with natural- or controlled-source excitations. The formalism allows for various types of parametrization of the 3-D conductivity distribution. Using this methodology one can readily obtain appropriate formulae for the specific sounding methods. To illustrate the concept we provide such formulae for two EM techniques: magnetotellurics and controlled-source sounding with vertical magnetic dipole as a source.

Description: Geomagnetic jerks are sudden changes of trend in the geomagnetic secular variation. The Earth's mantle behaves as a filter for the jerks, causing a delayed and a smoothed signal at the Earth's surface. Backus' mantle filter theory relies on approximating the impulse response function (IRF) of the mantle by a Gaussian. The advantage of this theory is the linear relation between jerks’ delay times and the mantle electrical conductivity, as expressed by kernels. However, the limitations of this theory arise when negative delay and/or smoothing times occur. The applicability of the mantle filter theory is examined by analysing the validity of the Gaussian as an approximation for the composite IRF (CIRF) at a given location. We show that the electrical conductivity of the lower mantle is mostly responsible for the jerk delay time. Alternating sign CIRFs might cause negative delay and/or smoothing times which prevents the use of the mantle filter theory. Adequate/inadequate Gaussian approximations to the CIRFs give small/large differences in the convolved jerk occurrence times. Most observatories yield positive time constants, but in most cases the difference in the jerk occurrence times exceeds 0.5 yr.

Description: A new method of estimating fluid flows near the core–mantle boundary is presented. Both the geomagnetic diffusion and dynamics in a boundary layer are taken into account under the no-slip condition at the core–mantle boundary. Below the boundary layer, the tangentially geostrophic constraint is imposed on the flow to be estimated. The radial component of the geomagnetic field just within the core is calculated from the radial component and its partial derivatives with respect to the radius, first, derived from continuity of the solenoidal field across the finite jump in electrical conductivity and, second, the geomagnetic diffusion. Below the boundary layer, fluid motions thus derived show some features common to those estimated on the basis of the frozen-flux approximation, and their accompanying vortices at mid and high latitudes are well correlated with upwellings or downwellings inside the boundary layer. The positional relation of upwellings and downwellings in equatorial regions inside and below the boundary layer suggests existence of columnar convective cells there. Intense magnetic flux spots seen in equatorial regions are likely to be produced not by downwellings at tops of cyclonic columns but by flux expulsion due to upwellings between adjacent columnar flows. The present new method leads to a new insight on core surface flows that have not been noted so far.

Description: In most numerical simulations of the Earth's core the dynamo resides outside the tangent cylinder and may be crudely classified as being of the α 2 type. In this region the flow comprises a sea of thin columnar vortices aligned with the rotation axis, taking the form of alternating cyclones and anticyclones. The dynamo is thought to be driven by these columnar vortices within which the flow is observed to be highly helical, helicity being a crucial ingredient of planetary dynamos. As noted in Davidson, one of the mysteries of this dynamo cartoon is the origin of the helicity, which is observed to be positive in the south and negative in the north. While Ekman pumping at the mantle can induce helicity in some of the overly viscous numerical simulations, it is extremely unlikely to be a significant source within planets. In this paper we return to the suggestion of Davidson that the helicity observed in the less viscous simulations owes its existence to helical wave packets, launched in and around the equatorial plane where the buoyancy flux is observed to be strong. Here we show that such wave packets act as a potent source of planetary helicity, constituting a simple, robust mechanism that yields the correct sign for h north and south of the equator. Since such a mechanism does not rely on the presence of a mantle, it can operate within both the Earth and the gas giants. Moreover, our numerical simulations show that helical wave packets dispersing from the equator produce a random sea of thin, columnar cyclone/anticyclone pairs, very like those observed in the more strongly forced dynamo simulations. We examine the local dynamics of helical wave packets dispersing from the equatorial regions, as well as the overall nature of an α 2 -dynamo driven by such wave packets. Our local analysis predicts the mean emf induced by helical waves, an analysis that rests on a number of simple approximations which are consistent with our numerical experiments, while our global analysis yields exact integral relationships between the mean emf induced by the wave packets and the large-scale dipole and azimuthal field. Combining these local and integral equations yields a kinematic model for an α 2 -dynamo driven by helical waves. Order-of-magnitude estimates based on these equations suggest that such a dynamo is indeed feasible.

Description: The Lower Danube Basin is one of the most important loess regions from Europe, which have provided excellent archives for long-term high-resolution palaeoclimate studies. The aim of this paper is to derive new information on the Middle–Late Pleistocene palaeoenvironment from a high resolution multiproxy assessment of the iron mineralogical composition at the Costinesti loess-palaeosol sequence located on the western Black Sea shore. It is the easternmost loess section in the Romanian loess region studied and its distinct pattern of the proxy records can be used to correlate the lower Danube loess to other key sites of the Moldavia and Ukraine loess regions. To investigate the climatic control on soft and hard ferromagnetic minerals we used several types of rock magnetic properties: magnetic susceptibility and its frequency dependence, anhysteretic remanent magnetization, isothermal remanent magnetization, hysteresis properties and FORC distributions, an unmixing model for isothermal remanent magnetization curves and high field (up to 8 T) isothermal remanence measurements. Our results show that the palaeosol horizons, formed during interglacials and climatically more favored periods of the Pleistocene, experienced pedogenic alteration, resulting in high amounts of superparamagnetic, single domain and pseudosingle domain magnetite/maghemite grains and hematite. The loess layers, formed during glacial periods, are mainly dominated by multidomain and/or pseudosingle domain oxidized magnetite and some hematite, all probably of aeolian origin. Goethite contribution is probably minor and constant both in loess and palaeosol horizons. We review the correlation of the loess sections from the lower Danube basin concluding that the new results can be interpreted as a support for the transition of a Mediterranean type climate to a steppe type climate in the last two interglacial periods in the western Black Sea. Because the pattern of magnetic susceptibility data from the lower Danube basin is changing relative fast with distance from the Black Sea shore, it probably reflects the local influence of the Black Sea on continental scale climatic oscillations during the last 600 ka. The values of background magnetic susceptibility of the Romanian loess-palaeosol sections indicate that the main source area of the dust changed during this climatic transition. Our analysis also shows that the age of the loess-palaeosol sections from the Eastern European low lands (Moldavia and Ukraine) must be revised to be in agreement with the chronostratigraphy of the sections from the Lower Danube Basin loess area.

Description: Statistical models of geomagnetic secular variations in terms of a Giant Gaussian Process (GGP) are compared with palaeontensity databases (Brunhes epoch, effusive rocks). We tested two GGP models TK and QC (suggested by Tauxe-Kent and Quidelleur-Courtillot) against 392 volcanic palaeointensity data, and found both models to be incompatible with the data. Even though at present there is no accepted GGP model of the Brunhes chron, we found interesting to clarify the details of this negative test. While such an incompatibility is clear enough for the TK model, because of its too small dipole coefficient, the incompatibility of the QC model was due the 12 per cent of the data with very small intensities; remaining 88 per cent of data is fully compatible with QC model QC. The troublesome, low-palaeointensity records belong to different sites and ages so their excess presence is not an effect of particular location or time. We found no evidence that the marked discrepancy between the empirical data and theoretical calculations on the QC model can be explained by the possibility that the low-intensity values belong to known excursions of the geomagnetic field. As an alternative, we conjecture that selected low-intensity values can possibly stem from acquisition of chemical remanent magnetization instead of the thermoremanence.

Description: A new method of estimating fluid flows near the core–mantle boundary is presented. Both the geomagnetic diffusion and dynamics in a boundary layer are taken into account under the no-slip condition at the core–mantle boundary. Below the boundary layer, the tangentially geostrophic constraint is imposed on the flow to be estimated. The radial component of the geomagnetic field just within the core is calculated from the radial component and its partial derivatives with respect to the radius, first, derived from continuity of the solenoidal field across the finite jump in electrical conductivity and, second, the geomagnetic diffusion. Below the boundary layer, fluid motions thus derived show some features common to those estimated on the basis of the frozen-flux approximation, and their accompanying vortices at mid and high latitudes are well correlated with upwellings or downwellings inside the boundary layer. The positional relation of upwellings and downwellings in equatorial regions inside and below the boundary layer suggests existence of columnar convective cells there. Intense magnetic flux spots seen in equatorial regions are likely to be produced not by downwellings at tops of cyclonic columns but by flux expulsion due to upwellings between adjacent columnar flows. The present new method leads to a new insight on core surface flows that have not been noted so far.

Description: Practical applications of controlled-source electromagnetic (EM) modelling require solutions for multiple sources at several frequencies, thus leading to a dramatic increase of the computational cost. In this paper, we present an approach using block Krylov subspace solvers that are iterative methods especially designed for problems with multiple right-hand sides (RHS). Their main advantage is the shared subspace for approximate solutions, hence, these methods are expected to converge in less iterations than the corresponding standard solver applied to each linear system. Block solvers also share the same preconditioner, which is constructed only once. Simultaneously computed block operations have better utilization of cache due to the less frequent access to the system matrix. In this paper, we implement two different block solvers for sparse matrices resulting from the finite-difference and the finite-element discretizations, discuss the computational cost of the algorithms and study their dependence on the number of RHS given at once. The effectiveness of the proposed methods is demonstrated on two EM survey scenarios, including a large marine model. As the results of the simulations show, when a powerful preconditioning is employed, block methods are faster than standard iterative techniques in terms of both iterations and time.

Description: Loess deposits in the arid Central Asia contain valuable information on the evolution of local aridification and dust sources in the Northern Hemisphere. Xinjiang is located in the eastern part of Central Asia and previous researches have revealed the complex enhancement of magnetic susceptibility in loess–paleosol sequences. However, systematic magnetic archives of loess deposit in this arid Asian interior are still far from adequate. In this study, magnetic parameters combined with nonmagnetic properties (granulometry and chromaticity) were analysed on a loess section in Shawan (SW), northwestern China. The section shares a similar magnetic composition with those in the Chinese Loess Plateau (CLP) as well as other sites in Xinjiang. Ferrimagnetic components (magnetite and maghemite) dominate the magnetic signal while the contribution of antiferromagnetic phases (like hematite and goethite) and paramagnetic portions are relatively low. There is no specific correlation between magnetic concentration and pedogenic intensity in the SW section. In general, magnetic enhancement was largely influenced by the paleowind intensity. However, a positive correlation between magnetic susceptibility and pedogenesis is observed in the upper part (0–3.5 m depths), which is characterized by a moderate wind intensity. Moreover, pedogenesis might be responsible for the enhancement of fine magnetic particles in paleosols. Magnetic properties are controlled by coarse magnetic particles in the pseudo-single domain state, but a coarse stable single domain phase was found in certain paleosol samples. The input of detrital fractions from a nearby dust source probably controlled the magnetic properties while a superparamagnetic fraction, which has been deemed as a product of pedogenesis in the CLP, is limited in the SW section. Caution is needed to employ magnetic susceptibility directly for paleoclimatic assessment because of its uncertainty in the Xinjiang loess. However, the ARM /SIRM ratio enhanced in paleosols and has more potential as an alternative index for the weakly pedogenic degree in this area.

Description: Mantle control on the geodynamo is often simulated using numerical dynamos with imposed outer boundary heat flux inferred from lower mantle tomography, assuming that seismic and thermal anomalies in the lowermost mantle are highly correlated. However, non-thermal effects might perturb this idealized linear seismic–thermal mapping. Here we use a probabilistic tomography model to isolate the thermal part of the seismic anomaly in order to impose a more realistic core–mantle boundary heat flux pattern on the outer boundary of numerical dynamo simulations. We demonstrate that on time average these dynamo models have more low-latitude convective and magnetic activity than corresponding models with conventional tomographic heat flux. In addition, the low-latitude magnetic flux and kinetic energy contributions are more time-dependent in the dynamo models with a probabilistic tomography heat flux, and thus may recover the observed latitudinal distribution of geomagnetic flux on the core–mantle boundary, which we propose as a morphological criterion for Earth-like dynamo models.

Description: The statistics of directional data on a sphere can be modelled either using the Fisher distribution that is conditioned on the magnitude being unity, in which case the sample space is confined to the unit sphere, or using the latitude–longitude marginal distribution derived from a trivariate Gaussian model that places no constraint on the magnitude. These two distributions are derived from first principles and compared. The Fisher distribution more closely approximates the uniform distribution on a sphere for a given small value of the concentration parameter, while the latitude–longitude marginal distribution is always slightly larger than the Fisher distribution at small off-axis angles for large values of the concentration parameter. Asymptotic analysis shows that the two distributions only become equivalent in the limit of large concentration parameter and very small off-axis angle.

Description: With advancing computational resources, 3-D inversion techniques have become feasible in recent years and are now a more widely used tool for magnetotelluric (MT) data interpretation. Galvanic distortion caused by small-scale near-surface inhomogeneities remains an obstacle for 3-D MT inversion which so far has experienced little attention. If not considered properly, the effect on 3-D inversion can be immense and result in erroneous subsurface models and interpretations. To tackle the problem we implemented inversion of the distortion-free phase tensor into the ModEM inversion package. The dimensionless phase tensor components describe only variations of the conductivity structure. When inverting these data, particular care has to be taken of the conductivity structure in the a priori model, which provides the reference frame when transferring the information from phase tensors into absolute conductivity values. Our results obtained with synthetic data show that phase tensor inversion can recover the regional conductivity structure in presence of galvanic distortion if the a priori model provides a reasonable assumption for the regional resistivity average. Joint inversion of phase tensor data and vertical magnetic transfer functions improves recovery of the absolute resistivity structure and is less dependent on the prior model. We also used phase tensor inversion for a data set of more than 250 MT sites from the central San Andreas fault, California, where a number of sites showed significant galvanic distortion. We find the regional structure of the phase tensor inversion results compatible with previously obtained models from impedance inversion. In the vicinity of distorted sites, phase tensor inversion models exhibit more homogeneous/smoother conductivity structures.

Description: Earth's liquid core hosts a diverse set of waves with periods ranging from days to thousands of years. One class of waves with periods of several decades is known to arise from an interplay between magnetic, Archimedes and Coriolis forces. These so-called MAC waves are thought to be relevant for interpreting historical fluctuations in the geomagnetic field. In this study, we show that MAC waves provide a good description of time-dependent zonal flow at the top of the core. The same collection of waves also offers a simple explanation for observed fluctuations in the dipole field. Both of these predictions require a stratified layer at the top of the core with a thickness of 130–140 km and a buoyancy frequency comparable to Earth's rotation rate. We extend these predictions to include changes in the length of day (LOD) and find that MAC waves can account for about half of the observed fluctuation at decadal periods. Larger fluctuations are possible when electromagnetic stresses couple MAC waves to flow in the interior of the core. In fact, an idealized model for the coupled motion overestimates the LOD fluctuations, probably reflecting limitations in this idealized model. Our results offer support for stable stratification at the top of the core and suggest a common origin for decadal fluctuations in the geomagnetic field and the LOD.

Description: A detailed magnetic analysis has been conducted on the borehole CJ-1 (172.3 m in length) from the Changjiang (Yangtze) River delta. Results show that the dominant magnetic carriers are magnetite and hematite. Palaeomagnetic results reveal that high-frequency changes in palaeomagnetic inclinations are tied significantly to abnormal anisotropy of magnetic susceptibility (AMS) patterns due to effects of the high-energy depositional environment in this region. On the basis of AMS patterns, doubtful palaeomagnetic directional anomalies can be distinguished from authentic palaeomagnetic excursions. Magnetostratigraphic results indicate that the Matuyama–Brunhes boundary (MBB) was recorded at a depth of ~152.5 m. The presence of several short-lived inclination anomalies implied that the sedimentation could be continuous even at the millennial timescale at certain depth intervals bracketing these fast geomagnetic events. In summary, our study provides new insights into constructing reliable magnetostratigraphy in the delta region. Moreover, our new magnetostratigraphy of the Changjiang River delta deposits will facilitate studies on the relevant long-term palaeoenvironmental evolution of the delta.

Description: Non-stoichiometric monoclinic 4C pyrrhotite (Fe 7 S 8 ) is a major magnetic remanence carrier in the Earth's crust and in extraterrestrial materials. Because of its low-temperature magnetic transition around 30 K also known as Besnus transition, which is considered to be an intrinsic property, this mineral phase is easily detectable in natural samples. Although the physical properties of pyrrhotite have intensively been studied, the mechanism behind the pronounced change in magnetization at the low-temperature transition is still debated. Here we report magnetization experiments on a pyrrhotite crystal (Fe 6.6 S 8 ) that consists of a 4C and an incommensurate 5C* superstructure that are different in their defect structure. The occurrence of two superstructures is magnetically confirmed by symmetric inflection points in hysteresis measurements above the transition at about 30 K. The disappearance of the inflection points and the associated change of the hysteresis parameters indicate that the two superstructures become strongly coupled to form a unitary magnetic anisotropy system at the transition. From this it follows that the Besnus transition in monoclinic pyrrhotite is an extrinsic magnetic phenomenon with respect to the 4C superstructure and therefore the physics behind it is in fact different from that of the well-known Verwey transition.

Description: Surface nuclear magnetic resonance technique, also called magnetic resonance sounding (MRS), is an emerging geophysical method that can detect the presence and spatial variations of the subsurface water content directly. In this paper, we introduce the MRS central loop geometry, in which the receiver loop is smaller than the transmitter loop and placed in its centre. In addition, using a shielded receiver coil we show how this configuration greatly increases signal-to-noise ratio and improves the resolution of the subsurface layers compared to the typically used coincident loop configuration. We compare sensitivity kernels for different loop configurations and describe advantages of the MRS central loop geometry in terms of superior behaviour of the sensitivity function, increased sensitivity values, reduced noise level of the shielded receiver coil, improved resolution matrix and reduced instrument dead time. With no extra time and effort in the field, central-loop MRS makes it possible to reduce measurement time and to measure data in areas with high anthropogenic noise. The results of our field example agree well with the complementary data, namely airborne electromagnetics, borehole data, and the hydrologic model of the area.

Description: Frictional heating during earthquake rupture reveals important information on earthquake mechanisms and energy dissipation. The amount of annealing varies widely and is, as yet, poorly constrained. Here we use magnetic susceptibility versus temperature measurements during cycling to increasingly elevated temperatures to constrain the maximum temperature a slip zone has experienced. The case study comprises sheared clay cored from the Japan Trench subduction plate-boundary fault zone (décollement), which accommodated the large slip of the 2011 M w 9.0 Tohoku-oki earthquake. The décollement was cored during the Integrated Ocean Drilling Program (IODP) Expedition 343, the Japan Trench Fast Drilling Project (JFAST). Heating signatures with estimated maximum temperatures ranging from ~300 to over 500 °C are determined close to the multiple slip surfaces within the décollement. Since it is impossible to tie a specific slip surface to a certain earthquake, thermal evidence for the cumulative effect of several earthquakes is unveiled. This as yet preliminary rock magnetic ‘geothermometer’ would be a useful tool to detect seismic heating along faults that experienced medium temperature rise, a range which is difficult to assess with other approaches.

Description: Feldspars are the most abundant rock-forming minerals in the Earth's crust, but their magnetic properties have not been rigorously studied. This work focuses on the intrinsic magnetic anisotropy of 31 feldspar samples with various chemical compositions. Because feldspar is often twinned or shows exsolution textures, measurements were performed on twinned and exsolved samples as well as single crystals. The anisotropy is controlled by the diamagnetic susceptibility and displays a consistent orientation of principal susceptibility axes; the most negative or minimum susceptibility is parallel to [010], and the maximum (least negative) is close to the crystallographic [001] axis. However, the magnetic anisotropy is weak when compared to other rock-forming minerals, 1.53 x 10 –9 m 3 kg –1 at maximum. Therefore, lower abundance minerals, such as augite, hornblende or biotite, often dominate the bulk paramagnetic anisotropy of a rock. Ferromagnetic anisotropy is not significant in most samples. In the few samples that do show ferromagnetic anisotropy, the principal susceptibility directions of the ferromagnetic subfabric do not display a systematic orientation with respect to the feldspar lattice. These results suggest that palaeointensity estimates of the geomagnetic field made on single crystals of feldspar will not be affected by a systematic orientation of the ferromagnetic inclusions within the feldspar lattice.

Description: Feldspars are the most abundant rock-forming minerals in the Earth's crust, but their magnetic properties have not been rigorously studied. This work focuses on the intrinsic magnetic anisotropy of 31 feldspar samples with various chemical compositions. Because feldspar is often twinned or shows exsolution textures, measurements were performed on twinned and exsolved samples as well as single crystals. The anisotropy is controlled by the diamagnetic susceptibility and displays a consistent orientation of principal susceptibility axes; the most negative or minimum susceptibility is parallel to [010], and the maximum (least negative) is close to the crystallographic [001] axis. However, the magnetic anisotropy is weak when compared to other rock-forming minerals, 1.53 x 10 –9 m 3 kg –1 at maximum. Therefore, lower abundance minerals, such as augite, hornblende or biotite, often dominate the bulk paramagnetic anisotropy of a rock. Ferromagnetic anisotropy is not significant in most samples. In the few samples that do show ferromagnetic anisotropy, the principal susceptibility directions of the ferromagnetic subfabric do not display a systematic orientation with respect to the feldspar lattice. These results suggest that palaeointensity estimates of the geomagnetic field made on single crystals of feldspar will not be affected by a systematic orientation of the ferromagnetic inclusions within the feldspar lattice.

Description: Palaeomagnetic data obtained from archaeological materials are used for reconstructions of the Earth's magnetic field of the past millennia. While many studies tested the reliability of this recorder for palaeointensity only a few studies did this for direction. The study presents an archaeomagnetic and rock magnetic investigation applied to an experimental pottery kiln, which was operated in 2003 to produce stone ware. This kind of high-quality pottery needs a temperature of at least 1160 °C. Shortly before heating of the kiln direct absolute measurements of the absolute geomagnetic field vector have been carried out close to it. After cooling of the kiln 24 oriented palaeomagnetic samples have been taken. Although Curie temperatures are about 580 °C, that is the typical temperature for magnetite, thermal as well as alternating field demagnetisations reveal also a considerable amount of hematite as magnetic carrier. This mixture of magnetite and hematite is dominated by pseudo-single domain grains. Demagnetisation removed in some cases weak secondary components, but in most cases the specimens carried a single component thermoremanent magnetisation. The mean characteristic remanent magnetisation direction agrees on 95 per cent confidence level with the directly measured field direction. Archaeointensity was obtained from five specimens with the Thellier–Coe method and with the multiple-specimen palaeointensity domain-state corrected method. Six of these specimens also provided a result of the Dekkers–Böhnel method, which overestimated the archaeointensity by about 9 per cent compared to the direct value, while after correction for fraction the value agrees very well. For the multiple-specimen palaeointensity domain-state corrected method only fractions between 25 and 75 per cent have been used and specimens showing alteration have been excluded. Above 450 °C many specimens showed alteration of the magnetic grains. Because median destructive temperatures were often above this value in most cases the fraction was less than 50 per cent. Nevertheless the obtained intensity (48.48 ± 0.24 μ) is on 95 per cent confidence level in agreement with the direct observation. Behaviour of the specimens during the Thellier-experiments was not ideal because of narrow unblocking temperature spectra and alteration. Nevertheless, the obtained mean archaeointensity is also in agreement with the direct field observation. Here the relative palaeointensity error is about 6 per cent and very high compared the multiple-specimen palaeointensity domain-state corrected method. The investigation demonstrates that a pottery kiln can provide a very precise estimate of the ancient geomagnetic field vector.

Description: We present a new model of the radial (1-D) conductivity structure of Earth's mantle. This model is derived from more than 10 yr of magnetic measurements from the satellites Ørsted, CHAMP, SAC-C and the Swarm trio as well as the global network of geomagnetic observatories. After removal of core and crustal field as predicted by a recent field model, we fit the magnetic data with spherical harmonic coefficients describing ring current activity and associated induction effects and estimate global C -responses at periods between 1.5 and 150 d. The C -responses are corrected for 3-D effects due to induction in the oceans and inverted for a 1-D model of mantle conductivity using both probabilistic and deterministic methods. Very similar results are obtained, consisting of a highly resistive upper mantle, an increase in conductivity in and beneath the transition zone and a conductive lower mantle. Analysis of the Hessian of the cost function reveals that the data are most sensitive to structures at depths between 800 and 1200 km, in agreement with the results obtained from the probabilistic approach. Preliminary interpretation of the inverted conductivity structure based on laboratory-based conductivity profiles shows that the recovered structure in the lower mantle either requires higher temperatures or the presence of material of high conductivity related to ponding of carbonate melts below the transition zone.

Description: A prominent gravity and geoid low lies just south of the Lofoten peninsula in northern Norway, partly coinciding with the location of Proterozoic granites of the Transscandinavian Igneous Belt and being offset by ca . 100 km to the highest topography of northern Norway. The study area extends both onshore and offshore and lies at the transition between Archaean and Proterozoic lithosphere. The Palaeoproterozoic basement has been overthrusted by the Palaeozoic nappes of the Caledonian orogen and now forms the passive margin of the NE Atlantic. We investigate the gravity anomaly performing combined 3-D geophysical-petrological forward modelling of the lithosphere and sublithospheric upper mantle using the interactive modelling program LitMod3D. We include variations in thickness and composition of the lithospheric mantle in order to include the effects on the rifted margin adjoining the Baltic craton. We compare three possible origins of the anomaly: (i) a low-density upper crust, representing the northward extension of the Transscandinavian Igneous Belt, (ii) a lower crustal source formed by a Moho depression and (iii) a thick, depleted lithospheric mantle of possibly Archaean origin. A similar, yet wider and stronger gravity anomaly is found on the conjugate margin in northeastern Greenland. A shallow crustal source is most consistent with the geophysical data sets. A respective source of the granitic belt, however, is difficult to reconcile with the regional geology both in Fennoscandia and Greenland. An additional contribution form a deeper source is suggested.

Description: Analysis of more than 10 yr of vertical magnetic transfer function (VTF) estimates obtained at 12 mid-latitude sites, located in different continents and tectonic settings, reveals significant temporal variations for a period range between approximately 250 and 2000 s. The most ubiquitous pattern is a seasonal modulation of the VTF element that relates the vertical to the horizontal north–south magnetic components ( Tx ), which shows a high peak around the June solstice (and a low peak around the December solstice) regardless of the location of the site. To quantify the influence of this source effect on the amplitude of VTFs, we modelled the temporal variations of VTFs using a function with dependence on season and magnetic activity indexes. The model shows that differences between VTF estimates obtained at seasonal peaks can reach 0.08 of Tx absolute values and that the effect increases with latitude and period. Seasonal variations are observed also in the VTF component relating vertical to horizontal east–west magnetic components ( Ty ), but here the pattern with respect to the geographic distribution is less clear. In addition to seasonal trends, we observe long-term modulations correlating with the 11-yr solar cycle at some sites. The influence of these external source effects should be taken into account, before attempting a geological interpretation of the VTFs. It can be misleading, for example, to combine or compare VTFs obtained from long-period geomagnetic data acquired at different seasons or years. An effective method to estimate and remove these source effects from VTFs is by comparison with temporal variations of VTFs from synchronously recorded data at sites located at similar latitude (〈5° of difference) and longitude (〈10° of difference). Source effects in temporal variations of VTFs can be identified as those patterns that exhibit similar amplitudes and significant correlation with the geomagnetic activity at all compared sites. We also provide a second-order polynomial which can be used to estimate the amplitude of the seasonal variations in the Tx component globally as a function of latitude.

Description: Directions recovered from palaeomagnetic samples are usually archived with some quantitative information about their precision, most often in the form of a so-called α 95 angle. Such angles are classically co-estimated with the recovered palaeomagnetic direction from a collection of samples providing individual estimates of this direction. In some instances, however, palaeomagnetic directions have to be inferred from a single sample in which case no α 95 angle can be recovered in this way. Fortunately, the progressive demagnetization techniques and principal component analysis universally used to recover directional information from single samples provide alternative measures of the error affecting the recovered direction, known as Maximum Angular Deviation (MAD) angles. These have so far only been considered as rough quality indicators. Here, however, we show that directions recovered in this way can be assumed to satisfy a Fisher distribution, and that the corresponding MAD angles can be rescaled into α 95 estimates by multiplying it by an appropriate factor, which only depends on the number of demagnetization steps used in the principal component analysis and on whether one relies on a standard or a so-called ‘anchored’ principal component analysis. These coefficients have been tabulated and practical recommendations for taking advantage of them outlined in the final section of the text. They provide simple means for users to produce much needed error bars on declination and inclination time series recovered from sedimentary long sequences.

Description: Broad-band magnetic susceptibility (MS) measurement, a novel magnetic method capable of quantifying a narrow grain size distribution (GSD) of superparamagnetic (SP) particles by measuring low-field MS at a number of frequency steps spanning four orders of magnitude, has been tested in a loess/palaeosol section at Luochuan in the Chinese Loess Plateau. The studied succession consists of sequences from the latest palaeosol unit (S0) to the upper part of the loess unit (L2), spanning the last glacial–interglacial cycle. Reconstructed GSDs consist of volume fractions on the order of 10 –24 m 3 , and the mean GSDs are modal but with distinctive skewness among the loess, the weakly developed palaeosols (weak palaeosols), and the mature palaeosols. This indicates that the mean volume of SP particles in this loess/palaeosol sequence tends to increase during the transition from loess -〉 weak palaeosol -〉 palaeosol, an indication of grain growth as pedogenesis progresses. Total frequency dependence, or TFD(per cent), the difference between 130 at the lowest (130 Hz) and 500k at the highest (500 kHz) frequencies normalized to 130 , is judged to be a more suitable index than previous frequency dependence parameters for the concentration of SP particles. TFD(per cent) has a strong correlation with 130 , showing a continuous ‘growth curve’ with the rate of increase being highest for the loess, moderate for the weak palaeosols, and saturated for the palaeosols. The characteristic curve suggests that smaller SP particles are preferentially formed in the earlier stage of pedogenesis rather than the later phase when even larger particles are formed in mature palaeosols. These results demonstrate that the broad-band MS measurement method will be useful for the quantitative assessment of magnetic nanoparticles in soils and sediments.

Description: Environmental magnetism can provide simple and fast methods to semi-quantify the contamination levels by heavy metals in urban areas from the relationships between the magnetic properties and the heavy metal concentrations. The aim of this study is to explore the interpretation of selected magnetic parameters as proxies for the evaluation of heavy metal contamination in urban dusts. Dust samples were collected from different districts of Nanjing, Southeast China. The magnetic properties of dusts were analysed and heavy metal (Fe, Al, As, Ba, Cd, Cr, Co, Cu, Hg, Mn, Mo, Ni, Pb, Sr, V and Zn) concentrations were measured. Magnetic minerals, mainly from anthropogenic activities, are dominated by coarse stable single domain-grained ferrimagnetic minerals. The – T curves indicated that the main magnetic minerals are magnetite and hematite. Magnetic properties exhibited similar spatial distribution and significant positive correlations with the concentrations of most metals and with the pollution load index. Magnetic properties had a strong link with the abundance of heavy metals derived from industrial emission and traffic activities, whereas there was a poor correlation with the concentrations of metals from multi-sources such as commercial/domestic sources and natural processes. Besides the low-frequency magnetic susceptibility ( LF ) and the saturation isothermal remanent magnetization, other magnetic parameters can be considered as efficient indicators in the assessment of heavy metal contamination. The linkage between the magnetic properties and heavy metal concentrations in street dusts depends on the sources and nature of their magnetic fraction. Environmental magnetism analyses may provide simple and rapid methods for the assessment of heavy metal contamination in urban street dusts.

Description: We present in this paper an alternative approach to modelling and inversion of surface nuclear magnetic resonance (sNMR) data that has numerous attractive characteristics. By considering the forward and inverse problem in the frequency-domain (FD) it is possible to heavily truncate the band-limited data set with no appreciable loss of information. Furthermore, it introduces a natural method of signal demodulating and sidesteps the rotating frame transformation, which can be difficult to apply to large volumes. By instead remaining in the laboratory frame, phase shifts due to off-resonance transmission and static variation in the Larmor frequencies become simpler to account for. For these reasons, we present the first practical scheme for comprehensively inverting the complex sNMR data set. The use of complex data allows for improved depth and decay constant resolution. However, it requires the knowledge of subsurface electrical conductivity, and by corollary the ability to accurately model those effects in the sNMR record. Finally, the complex FD signal is sensitive to dephasing of the sNMR record due to static magnetic field inhomogeneity. We illustrate the ability to detect dephasing in synthetic and field data and demonstrate an effective inversion which is able to solve for static dephasing effects. This results in better estimation of meaningful decay constants as well as reduced error in extrapolated initial signal amplitude. Both of these parameters are critical for using sNMR in applications requiring estimation of hydraulic parameters.

Description: We report a palaeomagnetic investigation of the last full geomagnetic field reversal, the Matuyama-Brunhes (M-B) transition, as preserved in a continuous sequence of exposed lacustrine sediments in the Apennines of Central Italy. The palaeomagnetic record provides the most direct evidence for the tempo of transitional field behaviour yet obtained for the M-B transition. 40 Ar/ 39 Ar dating of tephra layers bracketing the M-B transition provides high-accuracy age constraints and indicates a mean sediment accumulation rate of about 0.2 mm yr –1 during the transition. Two relative palaeointensity (RPI) minima are present in the M-B transition. During the terminus of the upper RPI minimum, a directional change of about 180 ° occurred at an extremely fast rate, estimated to be less than 2 ° per year, with no intermediate virtual geomagnetic poles (VGPs) documented during the transit from the southern to northern hemisphere. Thus, the entry into the Brunhes Normal Chron as represented by the palaeomagnetic directions and VGPs developed in a time interval comparable to the duration of an average human life, which is an order of magnitude more rapid than suggested by current models. The reported investigation therefore provides high-resolution integrated palaeomagnetic and radioisotopic data that document the fine details of the anatomy and tempo of the M-B transition in Central Italy that in turn are crucial for a better understanding of Earth's magnetic field, and for the development of more sophisticated models that are able to describe its global structure and behaviour.

Description: The anisotropy of magnetic susceptibility (AMS) of soft rocks was measured in order to distinguish between the effect of remote and local strain fields, determine the size of the related inelastic damage zone and resolve the fault-plane solutions of past earthquakes. The AMS fabrics were explored next to late Pleistocene syndepositional normal faults (total displacement up to ~3.5 m) that cross soft lacustrine rocks within the seismically active Dead Sea basin. ‘Deposition fabrics’ prevail meters away from the fault planes and are characterized by scattered maximum and intermediate principal AMS axes. ‘Deformation fabrics’ are detected up to tens of centimetres from the fault planes and are characterized by well-grouped AMS axes, in which one of the principal axes is parallel to the strike of the nearby fault. Variations in the AMS fabrics and magnetic lineations define the size of the inelastic damage zone around the faults. The results demonstrate that the deformation-driven magnetic fabrics and the associated inelastic damage zones are compatible with coseismic dynamic faulting and the effects of the local strain field during earthquakes. Most of the AMS fabrics show a conspicuous similarity to that of the fault-plane solutions, i.e., the principal AMS axes and instantaneous strain ellipsoids are coaxial. These results suggest a novel application of the AMS method for defining the shape and size of the damage zones surrounding dynamic faults and determining the full tensor of the local strain field.

Description: Very low frequency (VLF) electromagnetic waves that penetrate conductive magma-filled dykes generate secondary fields on the surface that can be used to invert for dyke properties. The model used for the interpretation calculates currents induced in a conductive strip by an inducing field that decays exponentially with depth due to the conductivity of the surrounding medium. The differential equations are integrated to give an inhomogeneous Fredholm equation of the second kind with a kernel consisting of a modified Bessel function of the second kind. Numerical methods are typically used to solve for the induced currents in the strip. In this paper, we apply a modified Galerkin–Chebyshev method, which involves separating the kernel into source and field spectra and integrating the source terms to obtain a matrix equation for the unknown coefficients. The incident wave is expressed as a Chebyshev series. The modified Bessel function is separated into a logarithmic singularity and a non-singular remainder, both of which are expanded in complex Chebyshev polynomials. The Chebyshev coefficients for the remainder are evaluated using a fast Fourier transform, while the logarithmic term and incident field have analytic series. The deconvolution then involves a matrix inversion. The results depend on the ratio of strip-size to skin-depth. For infinite skin-depth and a singular conductivity distribution given by $\tau _0 a/\sqrt{a^2 - z^2 }$ (where 0 is the conductance, a is the half-length and z the distance from the centre), Parker gives an analytic solution. We present a similar analytic series solution for the finite skin-depth case, where the size to skin depth ratio is small. Results are presented for different ratios of size to skin depth that can be compared with numerical solutions. We compare full-space and half-space solutions. A fit of the model to VLF data taken above a magma filled dykes in Hawaii and Mt Etna demonstrates that while properties such as depth to top, conductivity ratio, tilt and dip can be determined, the depth to bottom is indeterminate due to the exponential decay of the inducing field.

Description: Correlated data errors (noise) are common in magnetotelluric (MT) data, but MT inversions typically neglect error correlations without investigating the impact of this simplification on inversion results. This paper examines effects of neglecting frequency- and spatially correlated noise on MT inversion, based on a nonlinear Bayesian formulation which quantifies the uncertainties of inversion results in terms of marginal posterior probability densities and credibility intervals. To do so, data with frequency- and spatially correlated noise of differing degrees are generated for several layered (1-D) synthetic cases. Bayesian MT inversions are carried out for these data sets with and without accounting for error correlation (i.e. applying full and diagonal covariance matrices, respectively, in the inversion), and the results are compared. For cases with noise that is strongly correlated over frequency or space, parameter uncertainties estimated using the diagonal-covariance simplification (neglecting error correlations) are found to often be significantly underestimated compared with results computed using the full covariance matrix.

Description: We report a new method to infer continuous time-series of the declination, inclination and intensity of the magnetic field from archaeomagnetic data. Adopting a Bayesian perspective, we need to specify a priori knowledge about the time evolution of the magnetic field. It consists in a time correlation function that we choose to be compatible with present knowledge about the geomagnetic time spectra. The results are presented as distributions of possible values for the declination, inclination or intensity. We find that the methodology can be adapted to account for the age uncertainties of archaeological artefacts and we use Markov chain Monte Carlo to explore the possible dates of observations. We apply the method to intensity data sets from Mari, Syria and to intensity and directional data sets from Paris, France. Our reconstructions display more rapid variations than previous studies and we find that the possible values of geomagnetic field elements are not necessarily normally distributed. Another output of the model is better age estimates of archaeological artefacts.

Description: We present a rigorous method for interpolation of electric and magnetic fields close to an interface with a conductivity contrast. The method takes into account not only a well-known discontinuity in the normal electric field, but also discontinuity in all the normal derivatives of electric and magnetic tangential fields. The proposed method is applied to marine 3-D controlled-source electromagnetic modelling (CSEM) where sources and receivers are located close to the seafloor separating conductive seawater and resistive formation. For the finite-difference scheme based on the Yee grid, the new interpolation is demonstrated to be much more accurate than alternative methods (interpolation using nodes on one side of the interface or interpolation using nodes on both sides, but ignoring the derivative jumps). The rigorous interpolation can handle arbitrary orientation of interface with respect to the grid, which is demonstrated on a marine CSEM example with a dipping seafloor. The interpolation coefficients are computed by minimizing a misfit between values at the nearest nodes and linear expansions of the continuous field components in the coordinate system aligned with the interface. The proposed interpolation operators can handle either uniform or non-uniform grids and can be applied to interpolation for both sources and receivers.

Description: According to field observations, electromagnetic (EM) signals accompanying seismic waves can be recorded. The orders of magnitude of observed coseismic electric and magnetic signals associated with earthquakes are usually around 1 to 10 1 mV km –1 and 10 –2 to 1 nT, respectively. In this paper, we carry out numerical simulation of coseismic EM signals associated with seismic waves due to electrokinetic effect and compare with field observations. The seismic source is represented by a finite fault measuring 15 x 9 km 2 with a max slip displacement 1.5 m, corresponding to a M w 5.9 earthquake. While using the EM surface boundary condition of continuous horizontal EM components, the magnetic signals only accompany the late-arriving S waves at receiver near the ground surface. This is obviously different from field observations. Thus, we adopt another EM surface boundary condition which assumes the ground surface carries surface charge. For the used half-space model, a surface-charge density magnitude | Q sc | in excess of 10 –4 C m –2 is sufficient to make horizontal magnetic components clearly show up at the whole time duration of seismic waves. When | Q sc | increases, the contribution of surface-charge density to coseismic EM signals becomes more and more dominant comparing with that of the seismically induced streaming-current. We estimate the Q sc expected at the Earth's surface might be a value between –5 x 10 –4 and –0.1 C m –2 by the comparison between numerical results and field observations. The vertical magnetic signals only accompany the late-arriving seismic waves, because they are theoretically only induced by SH wave. The field observation results of vertical magnetic signals may be resulted from the scattering effect or the seismic dynamo effect. We conclude that electrokinetic effect combined with surface-charge assumption is one possible generation mechanism of the observed coseismic EM signals.

Description: The downward continuation of the observed geomagnetic field from the Earth's surface to the core–mantle boundary (CMB) is complicated due to induction and diffusion processes in the electrically conducting Earth mantle, which modify the amplitudes and morphology of the geomagnetic field. Various methods have been developed to solve this problem, for example, the perturbation approach by Benton & Whaler, or the non-harmonic downward continuation by Ballani et al. In this paper, we present a new approach for determining the geomagnetic field at the CMB by reformulating the ill-posed, one-sided boundary-value problem with time-variable boundary-value function on the Earth's surface into an optimization problem for the boundary condition at the CMB. The reformulated well-posed problem is solved by a conjugate gradient technique using the adjoint gradient of a misfit. For this purpose, we formulate the geomagnetic adjoint-state equations for efficient computations of the misfit gradient. Beside the theoretical description of the new adjoint-state method (ASM), the first applications to a global geomagnetic field model are presented. The comparison with other methods demonstrates the capability of the new method to determine the geomagnetic field at the CMB and allows us to investigate the variability of the determined field with respect to the applied methods. This shows that it is necessary to apply the ASM when investigating the effect of the Earth's mantle conductivity because the difference between the results of approximate methods (harmonic downward continuation, perturbation approach) and the rigorous ASM are of the same order as the difference between the results of the ASM applied for different mantle conductivities.

Description: Recent first principles calculations of the Earth's outer core thermal and electrical conductivities have raised their values by a factor of three. This has significant implications for geodynamo operation, in particular, forcing the development of a stably stratified layer at the core–mantle boundary (CMB). This study seeks to test the hypothesis of a stably stratified layer in the uppermost core by analysing geomagnetic observations made by the CHAMP satellite. An inversion method is utilized that jointly solves for the time-dependent main field and the core surface flow, where we assume the temporal variability of the main field, its secular variation (SV), to be entirely due to advective motion within the liquid outer core. The results show that a large-scale pure toroidal flow, consistent with a stably stratified layer atop the outer core, is not compatible with the observed magnetic field during the CHAMP era. However, allowing just a small amount of poloidal flow leads to a model fitting the observations satisfactorily. As this poloidal flow component is large scale, within a predominantly toroidal, essentially tangentially geostrophic flow, it is compatible with a stably stratified upper outer core. Further, our assumption of little or no diffusive SV may not hold, and a small amount of SV generated locally by diffusion might lead to a large-scale pure toroidal flow providing an acceptable fit to the data.

Description: Fifteen sites of lava flows from Martinique Island (FWI) have been selected to document the geomagnetic field in the Caribbean area over the past 2.5 Ma and further constrain the time-averaged field during this period. Identical characteristic directions were isolated using both AF and thermal stepwise demagnetization techniques in all flows. Nine mean-site directions have a normal polarity, while three others are reversed. The mean geomagnetic pole position obtained after reducing all directions to the same polarity is indistinguishable from the present north geographic pole. The dispersion is at least 8° larger than the values derived from the time-averaged field models and remains unexplained otherwise than resulting from the relatively small number of directions. The other three flows are characterized by large deviations from the expected north–south direction. One lava flow dated at 1.69 Ma (±0.02 Ma) is likely associated with a transitional field during the Gilsà subchron. The lava flow dated at 770 ka (±11 ka) coincides with the age of the Brunhes-Matuyama geomagnetic reversal and is also coeval with another intermediate flow of the same age found at Guadeloupe Island. The 617 ka (±52 ka) old unit is characterized by reversed directions that are evidently not related to the last reversal, but with other reversed polarity and transitional lava flows of the same age recorded, respectively at Mexico and La Palma island. We infer that the presence of reversed directions with the same age at distinct localities confirms that a short episode of reversed polarity has occurred during this period.

Description: We present the first inversion of geomagnetic Sq data in a framework of 3-D conductivity models. This problem has been considered as immensely difficult due to the complex spatial structure of the Sq source which, in addition, varies with season and solar activity. Recently, we developed a 3-D electromagnetic (EM) inversion solution that allows one to work in a consistent manner with data that originates from sources, irrespective of their spatial complexity. In this paper, we apply our 3-D EM inversion scheme to Sq data collected during the Australian Wide Array of Geomagnetic Stations project. Within this project, three components of the geomagnetic field were recorded between 1989 November and 1990 December with the use of 53 portable vector magnetometers. The instruments were distributed over the Australian mainland with an average spacing of 275 km between sites. Inverting this unique—in a sense of its spatial regularity, density and long operational time—data set, we recovered the 3-D conductivity distribution beneath Australia at upper mantle depths (100–520 km). This depth range was justified in the paper from resolution studies using checkerboard tests. In addition, we performed extensive modelling to estimate quantitatively the influence of various factors on Sq signals, namely from hypothetical anomalies, inaccuracy in the source, ocean, and model discretization. As expected, the ocean (coastal) effect appeared to be the largest so that it has to be accounted for during 3-D inversion as accurately as possible. Our 3-D inversions—of data from either single or multiple days—revealed a strong offshore conductor near the south-east coast of Australia, which persists at all considered depths. Varying in details, this anomaly is remarkably robust irrespective of the considered day(s). We compared our results to those obtained from a different inversion scheme and an independent induction data set, and observed encouraging similarity. Combination of the two results suggests, that this conductor continues to the base of the mantle transition zone at 660 km. The nature of this anomaly is not fully understood but one possible explanation is that it is attributed to a reservoir responsible for three hotspots in the region.

Description: A detailed palaeomagnetic, rock-magnetic and palaeointensity study has been carried out on a Miocene volcanic sequence which consists of 39 consecutive lava flows recording a polarity transition in La Gomera (Canary Islands, Spain). In addition, new 40 Ar/ 39 Ar ages were obtained in two flows, yielding 9.63 ± 0.06 Ma in the lower and 9.72 ± 0.08 Ma in the upper part of the sequence. Palaeomagnetic results allowed determining a ChRM direction in all studied lavas: The 25 lowermost flows of the sequence display normal polarity directions and above, a sequence of 14 flows correspond to a transitional geomagnetic regime. If considered together with palaeomagnetic results from a previous study, which were obtained on the flows immediately overlying the upper part of the sequence analysed in this work, these results indicate that the reversal recorded in the Hermigua sequence corresponds to the normal to reverse C4Ar2n to C4Ar3r polarity transition. The lower-lying 25 normal-polarity flows yield a mean direction D  = 359.6°, I  = 42.4° ( α 95  = 5.1°; k  = 33) which agrees well with the expected values. Above, a sequence of 14 flows displays a more irregular directional behaviour, including several transitional directions, suggesting the occurrence of a precursor to the transition. The presence of a virtual geomagnetic pole (VGP) cluster in the western Atlantic Ocean observed in this study coincides with previous records of Miocene transitions. Interestingly, this region corresponds to a near-radial flux centre of the present-day non axial dipole field. Angular dispersion of VGPs calculated for the 25 lowermost normal polarity flows of the sequence shows a lower than expected result. Palaeointensity determinations were carried out using a Thellier type double heating method. 27 of the 48 analysed samples measured yielded successful results. Mean VDMs mean values range from 1.1 ± 0.5 to 8.8 ± 0.9 x 10 22 Am 2 . The intensity values decrease significantly on approaching the directional transitional zone, suggesting an earlier start of the polarity transition in the intensity record, typical of a decreasing axial dipole.

Description: Presently, most loess/palaeosol magnetic susceptibility records are interpreted as following either the wind-vigour model or the pedogenic enhancement model. However redoxomorphic processes induced by waterlogging, often referred to gleying in the loess literature, are also known to alter loess deposits but their impact on loess/palaeosol magnetic susceptibility records has received little attention. The reported rock magnetic study aims to characterize the mineral magnetic response of loess to waterlogging-induced redoxomorphic processes, thus improving our understanding of mineral magnetic changes within loess deposits with respect to environmental and climate conditions. The Nussloch loess-palaeosol deposit (Rhine Valley, Germany) was targeted because it is one of the best-studied Pleniglacial deposits for Western Europe in which numerous tundra gley intervals have been identified. Moreover, a comprehensive high-resolution environmental magnetism study has never been undertaken for this site. Various rock magnetism experiments were conducted at both room and low temperatures to characterise the composition, concentration and relative magnetic grain size of the mineral magnetic assemblage. The relative changes in magnetic parameters within the investigated loess interval are primarily controlled by (1) varying concentrations of coarse-grained ferrimagnetic particles of detrital (aeolian) origin and (2) dissolution of fine-grained ferrimagnetic particles related to in situ post-depositional alteration promoted by waterlogging-induced redoxomorphic processes. Goethite is found to be ubiquitous throughout the studied interval and is argued to have both a primary (aeolian) and secondary ( in situ ) origin. We conclude, that redoxomorphic processes induced by waterlogging, if present, will hinder the interpretation of magnetic susceptibility variations within loess and palaeosol deposits following the expected relationships dictated by the wind-vigour and the pedogenic enhancement magnetism models.

Description: Palaeosecular variation (PSV), as estimated from the scatter of remanent magnetization directions or poles, can be used to shed light to processes in the geodynamo, and potentially, to model the growth of the solid inner core. To understand the temporal aspects of the geomagnetic field behavior in the far past, we have calculated the scatter of palaeomagnetic poles (S) from a set of 55 high-quality observations derived from a new Precambrian paleomagnetic data compilation. Our analysis supports former Phanerozoic and Precambrian analyses of PSV, which favour a lower reversal rate, a higher stability of the geodynamo and a substantially smaller size of the inner core prior to 1.5 billion years ago.

Description: A palaeomagnetic study has been carried out in 29 cores drilled at six different sites from the volcanic products of Lomo Negro eruption (El Hierro, Canary Islands, Spain). Systematic thermal and alternating field demagnetization of the samples’ natural remanent magnetization revealed a northward, stable palaeomagnetic direction similar in all the samples. Rock magnetic experiments indicate that this palaeomagnetic component is carried by a mixture of high-Ti and low-Ti titanomagnetite crystals typical of basaltic lithologies that have experienced a significant degree of oxyexsolution during subaerial cooling. The well constrained palaeomagnetic direction of Lomo Negro lavas was used to perform a palaeomagnetic dating of the volcanic event, using the SHA.DIF.14k global geomagnetic model restricted for the last 3000 yr. It can be unambiguously concluded that Lomo Negro eruption occurred well before the previously proposed date of 1793 AD, with three different age ranges being statistically possible during the last 3 ka: 115 BC-7 AD, 410–626 AD and 1499–1602 AD. The calibration of a previously published non-calibrated 14 C dating suggests a XVI c. date for Lomo Negro eruption. This conclusion leaves open the possibility that the seismic crisis occurred at El Hierro in 1793 AD was related to an intrusive magmatic event that either did not reach the surface or either culminated in an unregistered submarine eruption similar to the one occurred in 2011–2012 at the southern off-shore ridge of the island.

Description: We have performed an extended palaeomagnetic study of the Oslo Graben volcanics, compared to the study of half a century ago by van Everdingen, using modern techniques and a four times larger amount of sites, plus additional rock magnetic experiments. We conclude that the average direction ( D = 204.0, I = –37.9, k = 46.9, α 95 = 2.0) and associated palaeomagnetic pole ( = 48.3, = 155.5, K = 52.2, A 95 = 1.9) of the Krokskogen and Vestfold volcanics together are statistically identical to those of the earlier study. This gives confidence in the fact that older palaeomagnetic studies can be reliable and robust, even though methods have improved. Our larger number of samples, and better age constraints, enable us to separate the data into two major intervals: the younger, on average, Krokskogen area and the older Vestfold area. The results show firstly that palaeolatitudes are slightly higher than predicted by the latest apparent polar wander path (APWP) for Eurasia by Torsvik et al. These data support an early Permian Pangaea A configuration and do not necessitate a Pangaea B configuration. The larger data set also allows us to assess the distribution of the characteristic remanent magnetization directions of the Oslo Graben in terms of geomagnetic field behaviour, which were acquired during a long period of dominantly single polarity the Permo-Carboniferous Reversed Superchron (PCRS). The distributions show a significantly lower virtual geomagnetic pole (VGP) scatter at the observed (low) latitudes than expected from a compilation from lavas of the last 5 Myr. The data do however show excellent agreement with the scatter observed both during the Cretaceous Normal Superchron and the PCRS. A comparison of the directional distributions in terms of elongation is less discriminating, since the large errors in all cases allow a fit to the predicted elongation/inclination behaviour of the TK03.GAD model.

Description: We illustrate a new imaging method to estimate the depth to the sources of potential fields and the structural index. The method consists of applying the Depth from EXtreme Point (DEXP) transformation to the ratio (R) between two different-order partial derivatives of the field. While the scaling function of the potential field depends on the structural index, we show that the scaling function of R merely depends on the difference between the two used orders of differentiation. This allows three main features to be established for the DEXP transformation of R: (1) it is independent from the structural index; (2) the estimation of the source depths is fully automatic, simply consisting in the search of position of the extreme points of the DEXP image and (3) the structural index of each source is finally determined from the scaling function or the extreme points using the estimated depth. Besides the well-known characteristics of the DEXP transformation, such as high-resolution and stability, the DEXP transformation of R enjoys one more relevant feature: it can be applied to multisource cases, yielding simultaneously correct estimations of structural index and depth for each source in the same image. However, while the DEXP transformation is a linear transformation of the field, the DEXP transformation of R is non-linear, and a procedure is described to circumvent the non-linear effects. The method is tested with synthetic examples and the estimated source parameters show a good agreement with the true values. The method was applied also to real magnetic data from the Pima copper mine, Arizona, USA, Hamrawien area, Egypt and Cataldere, Bala district of Turkey. The results are consistent with the known information about the causative sources.

Description: Variational data assimilation (4DVar) is a powerful technique for tuning dynamic models to observations, in order not only to forecast future time evolution of the system, but to make inferences about quantities that are otherwise unconstrained by observation. We apply this technique, well-grounded in meteorology and oceanography, to the Earth's core where incompressible fluid motions in an electrically conducting medium are responsible for magnetic field generation. Our dynamic model's momentum equation neglects inertia such that the entire evolution depends only on the structure of the initial magnetic field; time evolution of the system is solely governed by the equation of magnetic induction. Nevertheless the dynamic system encompasses the effects of rotation, Lorentz forces and viscosity and aims to mimic a reasonable force-balance in the Earth's core. Building on the work of Li et al. , in order to optimize the data-fit subject to the dynamics, we further develop the mathematical structure of the adjoint equations of the system. We address the feasibility of recovering 3-D spatial properties of the system using only time-varying 2-D observations of different character. Using closed-loop testing, we demonstrate the retrievability of the initial state (and thus the entire trajectory) of the system over convective timescales, when sampling in regions in which magnetic induction dominates over diffusion. The results suggest the possibility of retrieving the entire trajectory of the dynamo system of the Earth using the 400-yr model of secular variation gufm1 .

Description: We investigate polarity reversals in the geodynamo using a rotating, convection-driven dynamo model. As the flow in rapidly rotating convection is dominated by columns aligned with the axis of rotation, the focus is on the dynamics of columnar vortices. By studying the growth of a seed magnetic field to a stable axial dipole field, we show that the magnetic field acts in ways that significantly enhance the relative helicity between cyclonic and anticyclonic vortices. This flow asymmetry is the hallmark of a dipolar dynamo. Strong buoyancy, on the other hand, offsets the effect of the magnetic field, establishing parity between positive and negative vortices. As the dipole field is deprived of the helicity required to support itself, the dynamo is pushed into a reversing state. This is a likely regime for polarity reversals in the Earth's core. The integral lengthscale at which buoyancy injects energy is not significantly different from the convective flow lengthscale, which implies that buoyancy does not feed vortices at the small scales where non-linear inertia is present. The lengthscale at which the Lorentz force acts in the reversing dynamo is small, which may allow the passive presence of non-linear inertia in the small scales.

Description: A new methodology to estimate magnetotelluric (MT) tensor relationships, called Estimation of Local transfer-functIons by Combining Interstation Transfer-functions (ELICIT), is proposed whereby the MT tensor relationships of the local site are derived using only interstation transfer functions. The MT impedance tensor and the geomagnetic transfer function at the local site are characterised by combining interstation tensor relationships between electric and magnetic fields at the local site with the horizontal magnetic fields acquired at a neighbouring site. The main property of the proposed method is that the employed interstation transfer functions are independently constrained, without the need to acquire the electric and the magnetic fields at the local site simultaneously to recover the local MT tensor relationships. Due to this property, the ELICIT method offers new possibilities for MT data acquisition and processing, providing significant improvements when the magnetic time-series at the local site are affected by local noise or are truncated. Error analysis shows that, even when magnetic fields are truncated, the quality of the results obtained following the ELICIT method are similar to those we would obtain if the magnetic fields had not been truncated. Another important property is that different neighbouring sites can be used to recover the tensor relationships at the local site. Averaging of results obtained using different neighbouring sites can be performed to improve the statistics. For our example data, when the ELICIT method is used to improve the statistics, errors of the estimates for periods between 1000 and 20 000 s periodicities are clearly reduced. All interstation transfer functions are calculated doing remote reference and the bootstrap method is used to compute the errors, when necessary. Long period magnetotelluric data acquired in the Pyrenees and in the Atlas Mountains in Morocco, and magnetic data provided by Fürstenfeldbruck magnetic observatory have been used to test the proposed ELICIT method, with positive results. Due to the lack of requirement that the electric and the magnetic fields of the local site be acquired simultaneously, the proposed method also offers new possibilities for MT data acquisition, optimizing the available instrumentation.

Description: The robust statistical model of a Gaussian core contaminated by outlying data that underlies robust estimation of the magnetotelluric (MT) response function has been re-examined. The residuals from robust estimators are systematically long tailed compared to a distribution based on the Gaussian, and hence are inconsistent with the robust model. Instead, MT data are pervasively described by the alpha stable distribution family whose variance and sometimes mean are undefined. A maximum likelihood estimator (MLE) that exploits the stable nature of MT data is formulated, and its two-stage implementation in which stable parameters are first fit to the data and then the MT responses are solved for is described. The MLE is shown to be inherently robust, but differs from the conventional robust estimator because it is based on a model derived from the data, while robust estimators are ad hoc , being based on the robust model that is inconsistent with actual data. Propriety versus impropriety of the complex MT response was investigated, and a likelihood ratio test for propriety and its null distribution was established. The Cramér-Rao lower bounds for the covariance matrix of proper and improper MT responses were specified. The MLE was applied to exemplar long period and broad-band data sets from South Africa. Both are shown to be significantly stably distributed using the Kolmogorov–Smirnov goodness of fit and Ansari-Bradley non-parametric dispersion tests. Impropriety of the MT responses at both sites is pervasive, hence the improper Cramér-Rao bound was used to estimate the MLE covariance. The MLE is shown to be nearly unbiased and well described by a Gaussian distribution based on bootstrap simulation. The MLE was compared to a conventional robust estimator, establishing that the standard errors of the former are systematically smaller than for the latter and that the standardized differences between them exhibit excursions that are both too frequent and too large to be described by a Gaussian model. This is ascribed to pervasive bias of the robust estimator that is to some degree obscured by their systematically large confidence bounds. Finally, a series of topics for further investigation is proposed.

Description: Broad-band and long period magnetotelluric (MT) data were acquired at 39 stations along five NNW-SSE profiles crossing the Iapetus Suture Zone (ISZ) in Ireland. Regional strike analyses indicate that the vast majority of the MT data is consistent with an assumption of a 2-D geo-electric strike direction. Strike is N52°E for the three easternmost profiles and N75°E for the two westernmost profiles; these directions correlate well with the observed predominant geological strike of the study region. 2-D inversions of the galvanic distortion-corrected TE and TM mode data from each profile are shown and discussed. As mapped geological variations between the neighbouring profiles suggest a heterogeneous subsurface, it is important to verify the robustness of the presence and geometries of prominent conductivity anomalies by employing 3-D forward and inverse modelling. A high conductivity layer (resistivity of 1–10 m), found at middle to lower crustal depths and presumed to be indicative of metamorphosed graphitic sediments rich in sulphides deposited during the convergence of the Laurentian and Avalonian continents, essentially constitutes the electrical signature of the ISZ. Shallow conductors observed are probably due to black shales that were widely deposited within the sedimentary accretionary wedge during Ordovician time. We interpret the moderately low resistivity at shallow depths from west to east across Ireland as indicative of an increase in maturity of the black shales in the easterly direction. From our conductivity models the southern extent of the ISZ is inferred to lie between the Navan Silvermines Fault and the Navan Tipperary Line, and shows clear resistivity contrast along all the profiles at the southern MT stations. The change in resistivity deduced from the 2-D models is spatially related to the composition of Lower Palaeozoic Ordovician, Silurian, Devonian and Carboniferous rocks. At upper mantle depths of about 60 km, a high conductivity block below the central MT stations is found to lie within the accretionary wedge of the Iapetus suture, and the location of the conductive anomaly corroborates reasonably well with the inferred spreading head of the putative Iceland plume-related magmatic intrusion. The low resistivity upper crust beneath the ISZ is indeed rich in Ordovician rocks with black shale content in the eastern as well as the central part; the western part is largely underlain by a highly resistive block of volcanic and metamorphosed rocks forming crystalline basement.

Description: This paper presents new directional archeomagnetic data from nine Meso-/Neolithic fireplaces, sampled in a cave shelter, at Arconciel, in western Switzerland. Rock magnetic measurements indicate a homogenous magnetic mineralogy in all fireplaces, with magnetite as the main magnetic carrier. The remanent magnetization is stable and generally shows one characteristic directional component. Nine new directions, which were obtained from Arconciel, are combined with 356 other archeomagnetic data from a circular area with a radius of 700 km around this site, to obtain a penalized least square spline fit for the past 9000 yr. We found in general good agreement with other local compilations, such as the Balkan curve, the regional SCHA.DIF.8k model and with lake sediments from UK, Fennoscandia and Switzerland. Nevertheless, a time lag of several centuries is observed for a declination maximum between the archeomagnetic spline fit and the other European data records around 5900 BC. This time lag is also observed in the Swiss lake sediment record; therefore we interpret this shift as a local feature of the Earth's magnetic field.

Description: In Christian religion the sunrise is of great symbolic importance. Therefore, many churches constructed in the Middle Ages point towards geographic East. Although ‘easting’ of churches actually refers to the alignment towards the azimuth of sunrise on the individual churches’ patron's day, deviation of nave alignment from the geographic East direction is often assumed to be caused by the use of magnetic compasses. Therefore, the church alignment could provide information about historical magnetic declination. We investigate 124 churches in Lower Austria and 68 in northern Germany to clarify this question as well as the ‘easting’ hypothesis. Church orientations are determined from georeferenced satellite images. Metadata such as the construction year, possible reconstructions and the church patron are gathered to determine the date when current church direction was appointed, and to perform sunrise calculations. However, due to uncertainties of construction years and the declining importance of orientation tradition after the 15th century several churches are excluded from the study. Thus, 32 churches with reliable metadata remain for evaluation in each region. The analysis reveals a preferred alignment of naves towards geographic East in Lower Austria and northern Germany. The construction and alignment of churches was often affected by the pre-existing buildings and streets or topography and natural surroundings. Therefore, deviations from geographic East are more likely caused by town or landscape. The mean deviations from magnetic East in both regions are large compared to the mean deviations from geographic East and sunrise azimuths. Hence, the use of compasses cannot be confirmed. Despite a few churches indicating orientation according to their patron's day, a general trend cannot be observed in the data reviewed.

Description: A new finite difference solution to the global geomagnetic induction problem is developed and tested, based on a modified Lorenz gauge of the magnetic vector and electric scalar potentials and implementing a novel, overset ‘Yin-Yang’ grid that avoids unnecessary mesh refinement at the geographic poles. Previously used in whole-earth mantle convection models, the overset grid is built from a pair of partially overlapping mid-latitude latitude–longitude (lat/lon) grids, one of which is rotated with respect to the other for complete coverage of the sphere. Because of this symmetry, only one set of finite difference templates is required for global discretization of the governing Maxwell equations, a redundancy that is exploited for computational efficiency and multithreaded parallelization. Comparisons between solutions obtained by the proposed method show excellent agreement with those obtained by independent integral equation methods for 1-D, 2-D and 3-D problem geometries. The computational footprint of the method is minimized through a (non-symmetric) matrix-free BiCG-STAB iterative solver which computes finite difference matrix coefficients ‘on the fly’ as needed, rather than pulling stored values from memory. Scaling of the matrix-free BiCG-STAB algorithm with problem size shows behaviour similar to that seen with the (symmetric) QMR algorithm used in the Cartesian case from which the present algorithm is based. The proposed method may therefore provide a competitive addition to the existing body of global-scale geomagnetic induction modelling algorithms, allowing for resource-efficient forward modelling as the kernel for large-scale computing such as inversion of geomagnetic response functions, computational hypothesis testing and parametric studies of mantle geodynamics and physiochemical state.

Description: We present a robust and scalable solver for time-harmonic Maxwell's equations for problems with large conductivity contrasts, wide range of frequencies, stretched grids and locally refined meshes. The solver is part of the fully distributed adaptive 3-D electromagnetic modelling scheme which employs the finite element method and unstructured non-conforming hexahedral meshes for spatial discretization using the open-source software deal.II. We use the complex-valued electric field formulation and split it into two real-valued equations for which we utilize an optimal block-diagonal pre-conditioner. Application of this pre-conditioner requires the solution of two smaller real-valued symmetric problems. We solve them by using either a direct solver or the conjugate gradient method pre-conditioned with the recently introduced auxiliary space technique. The auxiliary space pre-conditioner reformulates the original problem in form of several simpler ones, which are then solved using highly efficient algebraic multigrid methods. In this paper, we consider the magnetotelluric case and verify our numerical scheme by using COMMEMI 3-D models. Afterwards, we run a series of numerical experiments and demonstrate that the solver converges in a small number of iterations for a wide frequency range and variable problem sizes. The number of iterations is independent of the problem size, but exhibits a mild dependency on frequency. To test the stability of the method on locally refined meshes, we have implemented a residual-based a posteriori error estimator and compared it with uniform mesh refinement for problems up to 200 million unknowns. We test the scalability of the most time consuming parts of our code and show that they fulfill the strong scaling assumption as long as each MPI process possesses enough degrees of freedom to alleviate communication overburden. Finally, we refer back to a direct solver-based pre-conditioner and analyse its complexity in time. The results show that for multiple right-hand sides the direct solver-based pre-conditioner can still be faster for problems of medium size. On the other hand, it also shows non-linear growth in memory, whereas the auxiliary space method increases only linearly.

Description: The configuration of the Pangea supercontinent has been a topic of intense debate for almost half a century, a controversy that stems from discrepancies between the geology-based Pangea-A and the palaeomagnetically based Pangea-B. Recent palaeomagnetic compilations aimed at resolving this controversy have identified the poor quality of palaeomagnetic data from Gondwana for Permian times as a major obstacle. Specifically, the vast majority of Gondwanan poles come from sedimentary rocks that are prone to biases from compaction or are poorly dated. Here, we present a new palaeomagnetic pole for cratonic South America based on impact melts from the 254.7 ± 2.5 Ma Araguainha impact structure (AIS). The granite basement, the impact-generated melt sheet and veins were sampled at 28 sites (169 specimens) and provided a reliable palaeomagnetic record similar to that of volcanic rocks. Alternating field and thermal demagnetization indicate a stable characteristic remanent magnetization carried by both magnetite and haematite. All sites but one show a single palaeomagnetic direction of normal polarity with a mean direction of Dec = 357.4°; Inc = –38.9°; N  = 28; k  = 62.35; α 95  = 3.5°, yielding a palaeomagnetic pole (AIS) at Lat = –84.2; Lon = 326.6; K  = 83.5; A 95  = 3.6°; S B  = 9.6°. The new pole provides a firm constraint on the position of Gondwana which is consistent with the Pangea A configuration.

Description: Reconstructions of the Holocene geomagnetic field and how it varies on millennial timescales are important for understanding processes in the core but may also be used to study long-term solar-terrestrial relationships and as relative dating tools for geological and archaeological archives. Here, we present a new family of spherical harmonic geomagnetic field models spanning the past 9000 yr based on magnetic field directions and intensity stored in archaeological artefacts, igneous rocks and sediment records. A new modelling strategy introduces alternative data treatments with a focus on extracting more information from sedimentary data. To reduce the influence of a few individual records all sedimentary data are resampled in 50-yr bins, which also means that more weight is given to archaeomagnetic data during the inversion. The sedimentary declination data are treated as relative values and adjusted iteratively based on prior information. Finally, an alternative way of treating the sediment data chronologies has enabled us to both assess the likely range of age uncertainties, often up to and possibly exceeding 500 yr and adjust the timescale of each record based on comparisons with predictions from a preliminary model. As a result of the data adjustments, power has been shifted from quadrupole and octupole to higher degrees compared with previous Holocene geomagnetic field models. We find evidence for dominantly westward drift of northern high latitude high intensity flux patches at the core mantle boundary for the last 4000 yr. The new models also show intermittent occurrence of reversed flux at the edge of or inside the inner core tangent cylinder, possibly originating from the equator.

Description: Several recent studies have used palaeomagnetic estimates of the virtual axial dipole moment to construct a quantitative stochastic model for fluctuations and reversals in the Earth's dipole field. We investigate the physical significance of the terms in a standard stochastic (Langevin) model using output from a numerical geodynamo model. The first term, known as the drift term, characterizes the slow adjustment of the dipole field toward a time-averaged state. We find that the timescale for this slow adjustment is set by the magnetic decay time of dipole fluctuations. These fluctuations are typically be represented by the first few decay modes. The second term is often called the noise term because it characterizes the influence of short-period convective fluctuations in the core. We establish a connection between the noise term and the rms variation in magnetic induction. Applying these results to the palaeomagnetic field suggests that the rms variation in dipole generation exceeds the mean rate of generation. Such large fluctuations may be necessary to permit magnetic reversals. Palaeomagnetic estimates of the drift term favour a high electrical conductivity in the core. A lower bound on electrical conductivity is 0.6 x 10 6 S m –1 . Similarly, we establish an upper bound on turbulent magnetic diffusivity (0.8 m 2 s –1 ), although realistic estimates may be much less.

Description: We demonstrate that first differences of polar orbiting satellite magnetic data in the along-track direction can be used to obtain high resolution models of the lithospheric field. Along-track differences approximate the north–south magnetic field gradients for non-polar latitudes. In a test case, using 2 yr of low altitude data from the CHAMP satellite, we show that use of along-track differences of vector field data results in an enhanced recovery of the small scale lithospheric field, compared to the use of the vector field data themselves. We show that the along-track technique performs especially well in the estimation of near zonal spherical harmonic coefficients. Moreover, lithospheric field models determined using along-track differences are found to be less sensitive to the presence of unmodelled external field contributions and problems associated with the polar gap are ameliorated. Experiments in modelling the Earth's lithospheric magnetic field with along-track differences are presented here as a proof of concept. We anticipate that use of such along-track differences in combination with east–west field differences, as are now provided by the Swarm satellite constellation, will be important in building the next generation of lithospheric field models.

Description: We develop an singular value decomposition-based compression of the Green's function matrix of an electromagnetic integral equation forward solver for global geomagnetic induction, on top of an fast Fourier transform reduction of the system to a block-diagonal form. With this approach, the memory usage and CPU time of Krylov subspace iterative solutions are significantly reduced at a very small cost of accuracy, making the accelerated forward solver well suited for 3-D inversions as well as forward simulations with multiple sources.

Description: The Fisher distribution is central to palaeomagnetism but presents several problems when used to characterize geomagnetic field directions as observed in sequences of volcanic rocks. First, it introduces a shallowing effect when used to define the mean of any group of directional unit vectors. This is problematic because it can suggest the presence of persistent non-axial dipole components when none are present. More importantly, it fails to capture the observed ‘long tail’ in distributions of both directions and associated virtual geomagnetic poles in terms of angular distance from a central direction. To achieve a good fit to data, it therefore requires the introduction of a second distribution (and therefore the estimation of additional parameters) or the arbitrary removal of data. Here we present a new distribution to describe palaeomagnetic directions and demonstrate that it overcomes both of these problems, generating robust indicators of both the central direction (or pole position) and the spread of palaeomagnetic data as defined by unit vectors. Starting from the assumption that poles (or directions) have an expected colatitude, rather than a mean location, we derive the spherical exponential distribution. We demonstrate that this new distribution provides a good fit to palaeomagnetic data sets from seven large igneous provinces between 15 and 65 Ma and also those produced by numerical dynamo models. We also use it to derive a new shape parameter which may be used as a diagnostic tool for testing goodness of fit of models to data and use this to argue for a shift in geomagnetic behaviour between 5 and 15 Ma. Furthermore, we point out that this new statistic can be used to determine the most appropriate distribution to be used when constructing confidence limits for poles.

Description: Haematite pigment is a common constituent of sedimentary rocks, but its contribution to the natural remanent magnetization of rocks is poorly understood. Here, we describe magnetic properties of two distinct pigment types that produce a characteristic decorative ‘print stone’ found in the ~2.5 Ga Mount McRae Shale Formation, Hamersley Province, Western Australia. Distinct magnetic remanence directions observed in the Print Stone can be correlated to each pigment type. By comparison with the Australian apparent polar wander path, the remanence carried by uniformly distributed pigment can be dated to ~15–25 Ma, while two age options, the Mesoproterozoic (~1.5 Ga) or the middle Carboniferous (~320–310 Ma), are permissible for the remanence carried by the pigment responsible for the distinctive ‘newsprint’ pattern. Magnetic properties and demagnetization characteristics of the different pigment types overlap significantly, and thus are not predictive of the dominant remanence carrier. Magnetic characteristics of the uniformly distributed pigment vary significantly on short spatial scales. Strong local control on pigment formation raises the possibility that a primary remanent magnetization may survive locally in pockets within sedimentary red bed formations.

Description: Remanent magnetizations of magnetites between single-domain (SD) threshold size (0.1 μm) and 20 μm have SD-like intensities and coercivities. This paper shows for the first time that magnetite's induced magnetization also has pseudo-single-domain behaviour. The first part of the paper reports temperature-dependent initial susceptibility data, k 0 ( T ), of sized magnetites and assesses their granulometric potential. The second part transforms coercive force data, H c ( T ), for the same magnetites into simulated k 0 ( T ) curves. The third part considers k 0 ( T ) results of coarse-grained mafic rocks as candidate sources of deep-seated magnetic anomalies. High-temperature susceptibility k 0 measured with a Kappabridge for eight fractions of crushed natural magnetites (median sizes of 0.6, 1, 3, 6, 9, 14, 110 and 135 μm) shows a progressive increase in the height of the Hopkinson peak below the Curie point as grain size decreases. The trend is systematic and has granulometric potential in the 1–14 μm range. Self-demagnetization should produce almost flat k 0 ( T ) in grains larger than SD size but experimentally, well-defined Hopkinson peaks are not limited to the finest grains. 1-μm magnetites have a peak 1.5 times k 0 at 20 °C and 14-μm grains have a peak of 1.25. Only 110 and 135 μm grains have T -independent k 0 . Using an empirical relationship between coercive force H c and k 0 , H c ( T ) data for the sized magnetites were used to simulate k 0 ( T ) results. A hump in the k 0 heating curve around 250 °C was traced to annealing out internal strains, evident in H c data measured in first heatings. For sizes ≤6 μm, observed Hopkinson peaks were smaller than predicted, possibly because of a previously unrecognized grain-size dependence of the empirical constant relating H c and k 0 . Two crystalline rocks, a gabbro and a diabase, combine SD-like Hopkinson peaks and multidomain (MD) flat ramps in their k 0 ( T ) data. In the diabase, a Hopkinson peak is prominent in separated plagioclase grains containing submicron magnetite, but is masked in whole-rock data. The gabbro has a clear superposition of SD and MD k 0 ( T ) functions in its whole-rock data, with a Hopkinson peak of 1.35. If oceanic layer-3 gabbros have similar susceptibility enhancement above 500 °C, they could be more important magnetic anomaly sources than room-temperature k 0 measurements on dredged or fault-uplifted samples would suggest.