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: 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: 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: 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: 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 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: 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: 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: 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: 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: We present a new grand mean palaeomagnetic pole (Plong: 222.1°, Plat: –64.0°, A 95 : 2.6°, N = 49) for the ca. 1110 Ma Umkondo large igneous province (LIP) of the Kalahari Craton. New palaeomagnetic data from 24 sills in Botswana and compiled reprocessed existing data are used to develop a palaeomagnetic pole as the Fisher mean of cooling unit virtual geomagnetic poles (VGPs). The mean and its associated uncertainty provide the best-constrained pole yet developed for the province. Comparing data from individual cooling units allows for evaluation of palaeosecular variation at this time in the Mesoproterozoic. The elongation of the population of VGPs is consistent with that predicted by the TK03.GAD model lending support to the dipolar nature of the field in the late Mesoproterozoic. In our new compilation, 4 of 59 (~7 per cent) of the igneous units have northerly declinations while the rest are south-directed indicating that a geomagnetic reversal occurred during magmatic activity. Interpreting which of these polarities corresponds with a normal or reversed geomagnetic field relative to other continents can constrain the relative orientations between cratons with time-equivalent data. This interpretation is particularly important in comparison to Laurentia as it bears on Kalahari's involvement and position in the supercontinent Rodinia. The dominance of south-directed declinations within the Umkondo Province was previously used to suggest that these directions are the same polarity as reversed directions from the early magmatic stage of the Keweenawan Midcontinent Rift of Laurentia. Two Umkondo sills with northerly declinations have U-Pb baddeleyite ages of ca. 1109 Ma that are temporally close to dated Midcontinent Rift units having reversed directions. Based on this comparison, and palaeomagnetic data from younger units in the Kalahari Craton, we favour the option in which the sites with northerly declinations from the Umkondo Province correspond to the reversed polarity directions from the early magmatic stage in the Midcontinent Rift. This interpretation allows for the Namaqua-Natal metamorphic belt of Kalahari to be a conjugate to the Grenville margin of North America and for Kalahari to have become conjoined with Laurentia within the supercontinent Rodinia subsequent to Umkondo LIP magmatic activity.

Description: The northeastern Caribbean island arc, which materializes the boundary between the North American and Caribbean plates, is particularly exposed to large earthquakes and tsunamis. The low level of preparedness of a large part of its population and the lack of risk reduction provisions in public policies in many countries of the region put their population and economy at high risk in case of large telluric events. Here, we investigate the impact of three possible earthquake scenarios, consistent with the regional seismotectonic setting, on northern Haiti through inundation by tsunami waves. These scenarios simulate the effect of a M 8.0 earthquake on the Septentrional strike-slip fault (possibly similar to the 1842 earthquake), a M 8.1 earthquake on the offshore thrust fault system north of Haiti, and an earthquake rupturing a large portion of the offshore thrust fault system north of Haiti and the Dominican Republic. We calculate run-up heights along the northern coast of Haiti, in particular in the densely populated Cap Haitien. We find that the rupture of the offshore North Hispaniola thrust fault could result in wave heights up to 10 m with inundation up to 4 km inland, with only 10–15 min between ground shaking and the first wave arrivals. The city of Cap Haitien is particularly exposed, with potential flooding of most of the city and its suburbs, including the international airport. We also find that the historical reports available for the 1842 earthquake, when compared to our simulations, favor a rupture of the North Hispaniola thrust fault, although much uncertainty remains. If the 1842 earthquake did not rupture the Septentional fault offshore Haiti, then it is currently capable of at least a M w 7.7 earthquake, significantly larger than previously thought. The simulations presented here provide a basis for developing conservative maps of run-up heights that can be transferred, with added factors of safety, into practical implementation for tsunami preparedness and protection.

Description: We present the results of palaeomagnetic analysis on Late Bronge Age pottery from Santorini carried out in order to estimate the thermal effect of the Minoan eruption on the pre-Minoan habitation level. A total of 170 specimens from 108 ceramic fragments have been studied. The ceramics were collected from the surface of the pre-Minoan palaeosol at six different sites, including also samples from the Akrotiri archaeological site. The deposition temperatures of the first pyroclastic products have been estimated by the maximum overlap of the re-heating temperature intervals given by the individual fragments at site level. A new statistical elaboration of the temperature data has also been proposed, calculating at 95 per cent of probability the re-heating temperatures at each site. The obtained results show that the precursor tephra layer and the first pumice fall of the eruption were hot enough to re-heat the underlying ceramics at temperatures 160–230 °C in the non-inhabited sites while the temperatures recorded inside the Akrotiri village are slightly lower, varying from 130 to 200 °C. The decrease of the temperatures registered in the human settlements suggests that there was some interaction between the buildings and the pumice fallout deposits while probably the buildings debris layer caused by the preceding and syn-eruption earthquakes has also contributed to the decrease of the recorded re-heating temperatures.

Description: We have developed an algorithm, which we call HexMT, for 3-D simulation and inversion of magnetotelluric (MT) responses using deformable hexahedral finite elements that permit incorporation of topography. Direct solvers parallelized on symmetric multiprocessor (SMP), single-chassis workstations with large RAM are used throughout, including the forward solution, parameter Jacobians and model parameter update. In Part I, the forward simulator and Jacobian calculations are presented. We use first-order edge elements to represent the secondary electric field ( E ), yielding accuracy O ( h ) for E and its curl (magnetic field). For very low frequencies or small material admittivities, the E -field requires divergence correction. With the help of Hodge decomposition, the correction may be applied in one step after the forward solution is calculated. This allows accurate E -field solutions in dielectric air. The system matrix factorization and source vector solutions are computed using the MKL PARDISO library, which shows good scalability through 24 processor cores. The factorized matrix is used to calculate the forward response as well as the Jacobians of electromagnetic (EM) field and MT responses using the reciprocity theorem. Comparison with other codes demonstrates accuracy of our forward calculations. We consider a popular conductive/resistive double brick structure, several synthetic topographic models and the natural topography of Mount Erebus in Antarctica. In particular, the ability of finite elements to represent smooth topographic slopes permits accurate simulation of refraction of EM waves normal to the slopes at high frequencies. Run-time tests of the parallelized algorithm indicate that for meshes as large as 176 x 176 x 70 elements, MT forward responses and Jacobians can be calculated in ~1.5 hr per frequency. Together with an efficient inversion parameter step described in Part II, MT inversion problems of 200–300 stations are computable with total run times of several days on such workstations.

Description: Following the creation described in Part I of a deformable edge finite-element simulator for 3-D magnetotelluric (MT) responses using direct solvers, in Part II we develop an algorithm named HexMT for 3-D regularized inversion of MT data including topography. Direct solvers parallelized on large-RAM, symmetric multiprocessor (SMP) workstations are used also for the Gauss–Newton model update. By exploiting the data-space approach, the computational cost of the model update becomes much less in both time and computer memory than the cost of the forward simulation. In order to regularize using the second norm of the gradient, we factor the matrix related to the regularization term and apply its inverse to the Jacobian, which is done using the MKL PARDISO library. For dense matrix multiplication and factorization related to the model update, we use the PLASMA library which shows very good scalability across processor cores. A synthetic test inversion using a simple hill model shows that including topography can be important; in this case depression of the electric field by the hill can cause false conductors at depth or mask the presence of resistive structure. With a simple model of two buried bricks, a uniform spatial weighting for the norm of model smoothing recovered more accurate locations for the tomographic images compared to weightings which were a function of parameter Jacobians. We implement joint inversion for static distortion matrices tested using the Dublin secret model 2, for which we are able to reduce nRMS to ~1.1 while avoiding oscillatory convergence. Finally we test the code on field data by inverting full impedance and tipper MT responses collected around Mount St Helens in the Cascade volcanic chain. Among several prominent structures, the north–south trending, eruption-controlling shear zone is clearly imaged in the inversion.

Description: Palaeointensity variation is investigated for an inferred time period spanning from 2.34 to 1.96 Ma. Twenty-nine consecutive lava flows are sampled along cliffs 350 m high generated by normal faulting on the Dobi section of Afar depression, Ethiopia. Magnetostratigraphy and K–Ar measurements indicate a lava sequence of R–N–R–N geomagnetic field polarities in ascending order; the lower normal polarity is identified as the Réunion Subchron. Reliability of palaeomagnetic data is ascertained through careful thermal demagnetization and by the reversal test. The Tsunakawa–Shaw method yielded 70 successful palaeointensity results from 24 lava flows and gave 11 acceptable mean palaeointensities. Reliability in palaeointensity data is ascertained by the similar values obtained by the IZZI–Thellier method and thus 11 reliable mean values are obtained from our combined results. After the older reverse polarity with the field intensity of 19.6 ± 7.8 μT, an extremely low palaeointensity period with an average of 6.4 μT is shown to occur prior to the Réunion Subchron. During the Réunion Subchron, the dipole field strength is shown to have returned to an average of 19.5 μT, followed by second extreme low of 3.6 μT and rejuvenation with 17.1 ± 5.3 μT in the younger reverse polarity. This ‘W-shape’ palaeointensity variation is characterized by occurrences of two extremely weak fields lower than 8 μT prior to and during the Réunion Subchron and a relatively weak time-averaged field of approximately 15 μT. This feature is also found in sedimentary cores from the Ontong Java Plateau and the north Atlantic, indicative of a possibly global geomagnetic field phenomenon rather than a local effect on Ethiopia. Furthermore, we estimate a weak virtual axial dipole moment of 3.66 (±1.85) x 10 22 Am 2 during early stage of the Matuyama Chron (inferred time period of 2.34–1.96 Ma).

Description: Adopting a model with two half-spaces that consist of solid and porous materials, we numerically investigate the seismoelectric conversion at the solid–porous interface. First, the wave fields in a low-porosity two-layer model are compared with those in a homogeneous full-space model. The consistency of seismic waves is a validation of our program. We are interested in the quasi-coseismic electromagnetic (EM) signals recorded in the solid area near the interface because they seemingly accompany seismic waves. Then, further numerical simulations on an ordinary two-layer model are conducted. On the basis of time slice snapshots and theoretical analysis, we determine that quasi-coseismic EM signals are essentially non-coseismic EM fields, which include radiation and evanescent EM waves. Evanescent EM waves are induced by the seismic waves that arrive at the interface with the incident angle greater than the critical angle. These waves decay faster than radiation EM waves when moving away from the interface. In the porous layer, evanescent EM waves can hardly be recognized unless they are separated from coseismic EM signals. This finding can be the reason why evanescent EM waves have not been identified in previous seismoelectric studies. Awareness of the fact that seismoelectric conversion at an interface can generate evanescent and EM waves is likely to result in a comprehensive understanding and improved interpretation of the seismoelectric coupling phenomenon.

Description: We present a new technique for modelling the global lithospheric magnetic field at Earth's surface based on the estimation of equivalent potential field sources. As a demonstration we show an application to magnetic field measurements made by the CHAMP satellite during the period 2009–2010 when it was at its lowest altitude and solar activity was quiet. All three components of the vector field data are utilized at all available latitudes. Estimates of core and large-scale magnetospheric sources are removed from the measurements using the CHAOS-4 model. Quiet-time and night-side data selection criteria are also employed to minimize the influence of the ionospheric field. The model for the remaining lithospheric magnetic field consists of magnetic equivalent potential field sources (monopoles) arranged in an icosahedron grid at a depth of 100 km below the surface. The corresponding model parameters are estimated using an iteratively reweighted least-squares algorithm that includes model regularization (either quadratic or maximum entropy) and Huber weighting. Data error covariance matrices are implemented, accounting for the dependence of data variances on quasi-dipole latitude. The resulting equivalent source lithospheric field models show a degree correlation to MF7 greater than 0.7 out to spherical harmonic degree 100. Compared to the quadratic regularization approach, the entropy regularized model possesses notably lower power above degree 70 and a lower number of degrees of freedom despite fitting the observations to a very similar level. Advantages of our equivalent source method include its local nature, the possibility for regional grid refinement and the production of local power spectra, the ability to implement constraints and regularization depending on geographical position, and the ease of transforming the equivalent source values into spherical harmonics.

Description: Previous works revealed a close relationship between magnetic susceptibility (MS) and heavy metal (HM) contents originating from industrial sources. However, despite general statements on the usefulness of magnetic mapping, the benefit of this procedure for geochemistry was not quantified yet. We present a study on fly ash pollution in soil around a coal-burning power plant complex and simulate a stepwise approach of magnetic pre-screening and subsequent targeted sampling for chemical analysis. The aim of this study is not to discuss correlations between MS and HM, but to show that a combined stepwise magnetic-chemical approach is the most efficient way for outlining HM contamination. In order to provide quantitative evidence, we explored map similarities of spatial HM distributions based on magnetochemical data and chemical data only. We determined 3-D triangular planes defined by categorized HM values at the sampling coordinates and calculated the average dihedral angle of the normal vectors as a similarity result. The study shows that the ‘Targeted’ HM map (selection of 30 sites based magnetic pre-screening) has a higher similarity with the ‘True’ Pollution HM map (85 sites) than HM maps resulting from site selections (30 sites) without using magnetic pre-screening information.

Description: Archaeological baked clay remains provide valuable information about the geomagnetic field in historical past, but determination of the geomagnetic field characteristics, especially intensity, is often a difficult task. This study was undertaken to elucidate the reasons for unsuccessful intensity determination experiments obtained from two different Bulgarian archaeological sites (Nessebar – Early Byzantine period and Malenovo – Early Iron Age). With this aim, artificial clay samples were formed in the laboratory and investigated. The clay used for the artificial samples preparation differs according to its initial state. Nessebar clay was baked in the antiquity, but Malenovo clay was raw, taken from the clay deposit near the site. The obtained artificial samples were repeatedly heated eight times in known magnetic field to 700 °C. X-ray diffraction analyses and rock-magnetic experiments were performed to obtain information about the mineralogical content and magnetic properties of the initial and laboratory heated clays. Two different protocols were applied for the intensity determination—Coe version of Thellier and Thellier method and multispecimen parallel differential pTRM protocol. Various combinations of laboratory fields and mutual positions of the directions of laboratory field and carried thermoremanence were used in the performed Coe experiment. The obtained results indicate that the failure of this experiment is probably related to unfavourable grain sizes of the prevailing magnetic carriers combined with the chosen experimental conditions. The multispecimen parallel differential pTRM protocol in its original form gives excellent results for the artificial samples, but failed for the real samples (samples coming from previously studied kilns of Nessebar and Malenovo sites). Obviously the strong dependence of this method on the homogeneity of the used subsamples hinders its implementation in its original form for archaeomaterials. The latter are often heterogeneous due to variable heating conditions in the different parts of the archaeological structures. The study draws attention to the importance of multiple heating for the stabilization of grain size distribution in baked clay materials and the need of elucidation of this question.

Description: Magnetic minerals acquire a remanent magnetization when cooled or warmed in an applied field H through a phase transition. This paper reports the first observations of transition cooling and warming remanent magnetizations (TrCRM and TrWRM) in pyrrhotite (Fe 7 S 8 ) and hematite (αFe 2 O 3 ), with phase transitions at 32 K (Besnus transition) and 250 K (Morin transition). TrWRM was produced by warming in a 2 mT field from 10 to 300 K and TrCRM by cooling in 2 mT from 300 to 10 K. In both experiments, magnetization M was measured at 5 K intervals. The largest changes in M occurred over a 20 K interval flanking the Besnus or Morin transition but substantial changes also occurred away from any transition, presumably due to continuous changes in anisotropy. The triclinic low-temperature phase of pyrrhotite has high remanence efficiency compared to the room-temperature monoclinic phase. When H was zeroed at 300 K, only 35–40% of M survived as TrWRM, but in TrCRM experiments when H was zeroed at 10 K, practically 100% of M was preserved. After zero-field cooling to 10 K, 50–60% of TrWRM was retained but 〈15% of TrCRM survived zero-field warming to 300 K. In hematite, both TrWRM of the room-temperature spin-canted phase and TrCRM of the low-temperature defect phase exist, as shown by net remanence production following complete warming–cooling or cooling–warming cycles between 10 and 300 K. Details of remanence acquisition/demagnetization in crossing the Morin transition could not be resolved, however. Symmetries were observed for both minerals. Single-domain size pyrrhotite crystals had mirror-image field-off and field-on warming curves below 32 K, during TrWRM production/TrCRM loss at the 32 K transition, and up to 150 K. Submicron single-domain hematite had almost complete symmetry between field-off and field-on curves during cooling as well as warming, both below and above the Morin transition but not in the transition region itself. This is the first time that mirror-image symmetries, observed previously for pseudo-single-domain and multidomain magnetite, have been documented for mineral grains in a single-domain state.

Description: A computational method is given for solving the forward modelling problem for transient electromagnetic exploration. Its key features are the discretization of the quasi-static Maxwell's equations in space using the first-kind family of curl-conforming Nédélec elements combined with time integration using rational Krylov methods. We show how rational Krylov methods can also be used to solve the same problem in the frequency domain followed by a synthesis of the transient solution using the fast Hankel transform, and we argue that the pure time-domain solution is more efficient. We also propose a new surrogate optimization approach for selecting the pole parameters of the rational Krylov method which leads to convergence within an a priori determined number of iterations independent of mesh size and conductivity structure. These poles are repeated in a cyclic fashion, which, in combination with direct solvers for the discrete problem, results in significantly faster solution times than previously proposed schemes.

Description: In 2010, a wide-angle seismic reflection/refraction profile was acquired along the Hecataeus Rise, an area of shallow seabed immediately south of Cyprus in the eastern Mediterranean. The profile crossed from the Hecataeus Rise, through the Cyprus Arc to the Levantine Basin beyond. Due to the short length of the profile and the corresponding lack of deep ray coverage, velocity modelling was complemented by gravity modelling to gain constraints on deep crustal structure. The resultant model reveals velocities for the Hecataeus Rise that show no evidence of shallow ophiolites like those seen on mainland Cyprus, and the velocities are not diagnostic of a unique crustal affinity. Low-velocity sediments make up at least 7 km of the upper structure of Hecataeus Rise and these sediments overlie a two-layered crust. From the gravity modelling, the combined sediments and crust of Hecataeus Rise appear to be thinner than the Eratosthenes Seamount block to the southwest. A high-velocity lower crustal block is modelled under the seaward edge of Hecataeus Rise and, based on the gravity modelling, is inferred to extend landwards beneath the Rise. Similar high-velocity blocks were identified on the southwestern edge of Hecataeus Rise along nearby refraction lines and were interpreted as remnant Tethyan oceanic crust, foundered in the Cyprus Arc, along which subduction has ceased in this area. Given the thin two-layered crust beneath a thick accumulation of sediments modelled for Hecataeus Rise, we interpret that Hecataeus Rise represents a collage of oceanic fragments, accreted together within the failed subduction zone. Outboard of the crust of Hecataeus Rise, a 5-km deep low-velocity basin, possibly an accretionary wedge, is imaged that appears to correspond with the Cyprus Arc deformation zone imaged on both coincident and along-strike seismic reflection lines. A similar and wider feature is observed on seismic refraction lines to the west and combined, these may be revealing an eastward tapering zone of crustal deformation. To the south of the profile, the Levantine Basin appears undisturbed by the collision to the north and exhibits a uniform and homogeneous velocity structure.

Description: The Rogaland Igneous Complex (RIC) in southern Norway intruded into post-Sveconorwegian granulite facies crust ~930 Ma. It includes three massif anorthosites, several small leuconorite bodies and the ~7 km thick norite-quartz mangerite layered Bjerkreim-Sokndal (BKS) intrusion. The intrusion consists of five rhythmic megaunits created by repeated magma influxes topped by a transition zone and more evolved mangerites and quartz mangerites. Over 70 palaeomagnetic sites have been collected in the BKS, sampling all the megacyclic subunits and overlying mangerites. Remanence within the BKS is held in hemo-ilmenite-only rocks (lower parts of the megacyclic units), mixed hemo-ilmenite and magnetite rocks (upper parts of the lower megacyclic units) and magnetite only rocks in the upper highest megacyclic unit and overlying mangerites. Due to the different oxides present magnetic susceptibility varies over four orders of magnitude with a bimodal distribution (mean susceptibility of 6.4 x 10 –3 SI for hemo-ilmenite rocks, and 8.7 x 10 –2 SI for magnetite rocks). NRM values do not show a strong bimodal distribution as many of the rocks lacking magnetite have hemo-ilmenite with strong lamellar magnetism; average NRM for the entire suite is 8.83 A m –1 . All sites within the cyclic part of the intrusion have stable remanence and produce well-clustered site means. Samples from the upper mangerite rocks, dominated by MD magnetite, are commonly unstable and not all sites provide acceptable data. Mean directions for 66 sites spanning the entire intrusion are I = –73.5°, D = 303.4°, with α 95 = 3.7° and k = 24. The resulting pole position is at 35.9°S and 217.9°E, with a palaeolatitude for this part of Baltica of –59.4°. Examination of the magnetic mineralogy combined with geochronology for RIC rocks and cooling rates for the region yields an age of magnetization of 916 Ma. Metamorphic country rocks yield similar directions at least 10 km from the contact, confirming the presence of a contact aureole around the intrusion. Comparison to other early Neoproterozoic palaeomagnetic poles from southern Scandinavia confirms the high southern latitude position of Baltica at this time, and combined with the few other ~900 Ma poles from Baltica affirm the presence of the Rodinia supercontinent at this time.

Description: The Internal Sierras (IS) in the southern margin of the Western and Central Axial Zone (Southern Pyrenees) are affected by a syn-orogenic remagnetization that provides information to reconstruct deformation geometries at the time of acquisition of magnetization. Furthermore, the IS structure changes strike along its structural trend, from ~N120 to 130°E in the western and eastern margins to ~N070–090°E in the central part. Palaeomagnetic techniques have been used to (i) accurately define the timing of remagnetization with regard to deformation and (ii) determine if the along-strike trend variation in the IS was induced by deformation and thrust emplacement during the Pyrenean compression or, on the contrary, was the result of a primary orientation controlled by structures inherited from pre-orogenic times. From 23 new palaeomagnetic sites, collected in Upper Cretaceous marls and marly limestones, two meaningful and stable palaeomagnetic components were resolved, principally carried by magnetite: (1) a lower-temperature component (B) that unblocks between 200 °C and 325–400 °C and (2) a higher-temperature component (C) that has been successfully isolated by means of combined thermal (up to 400 °C) and AF demagnetization (generally up to 50–100 mT). The B component is a late remagnetization that post-dates folding and emplacement of basement thrust sheets in the IS (mainly the Gavarnie thrust). It supports small but statistically significant clockwise rotations in the western part of the IS (from +18 to +26°). These rotations can be attributed to the westwards shortening decrease in the thrust system below the Gavarnie unit that results from its along-strike structural change, with a higher number of basement thrusts to the east. The C component has been interpreted as an early remagnetization, based on the results of conglomerate and fold tests. This component predates basement thrusting and is diachronous across the study area: reverse and normal polarities dominate in the eastern and western margins of the IS, respectively. New and previous palaeomagnetic data point out that curvature in the IS is probably a primary feature and the along-strike change in their trend could be interpreted as the result of basement geometrical features inherited from Variscan, Late Variscan or Mesozoic times. A complex, multi-episodic remagnetization probably related to burial and deformation processes occurred during Eocene times.

Description: A comprehensive magnetic field model named CM5 has been derived from CHAMP, Ørsted and SAC-C satellite and observatory hourly-means data from 2000 August to 2013 January using the Swarm Level-2 Comprehensive Inversion (CI) algorithm. Swarm is a recently launched constellation of three satellites to map the Earth's magnetic field. The CI technique includes several interesting features such as the bias mitigation scheme known as Selective Infinite Variance Weighting (SIVW), a new treatment for attitude error in satellite vector measurements, and the inclusion of 3-D conductivity for ionospheric induction. SIVW has allowed for a much improved lithospheric field recovery over CM4 by exploiting CHAMP along-track difference data yielding resolution levels up to spherical harmonic degree 107, and has allowed for the successful extraction of the oceanic M 2 tidal magnetic field from quiet, nightside data. The 3-D induction now captures anomalous Solar-quiet features in coastal observatory daily records. CM5 provides a satisfactory, continuous description of the major magnetic fields in the near-Earth region over this time span, and its lithospheric, ionospheric and oceanic M 2 tidal constituents may be used as validation tools for future Swarm Level-2 products coming from the CI algorithm and other dedicated product algorithms.

Description: The observation of electromagnetic signals by ocean bottom electromagnetometers in association with the Tohoku tsunami of 2011 March 11 has raised the opportunity to re-examine the physics of motional induction due to oceanic long waves in the framework of 1-D magnetotellurics (MT). Although a propagating tsunami has a complex structure, the induced electromagnetic field can be simply approximated as a plane wave (though a simple thin-sheet approximation is not valid at higher frequencies). We found that the MT impedance due to a surface gravity wave (or the ‘motional impedance’) is influenced largely by the dispersion of the wave if the period is sufficiently short or the electrical conductivity of the seabed is low. The tipper due to the motional induction (or the ‘motional tipper’) and motional impedance are essentially identical if the underneath structure is 1-D. It would be possible to estimate the motional impedance and tipper from the observed ocean bottom electromagnetic field at the time of passing of a tsunami. The wave amplitude must be much greater than several tens of centimetres for the motional impedance and tipper estimation to be free from the effects of external sources. However, the obtained motional impedance and tipper will mostly represent the property of the wave and use of them may not be suitable to discuss the subseafloor conductivity structure.

Description: The lacustrine deposits in Heqing basin provide an excellent archive for long-term high-resolution palaeoclimate studies in the monsoon-dominated southeastern Tibetan Plateau region. In this study, we investigate the climatic significance of magnetic parameters for analysing the variability of the past climate. For this, we performed comprehensive time-series and statistical analysis on previously published proxy data from a 168-m-long drill core (Core-HQ) that spans ~900–30 ka, comprising mainly magnetic parameters and carbonate content (CC). Moreover, we investigated magnetic properties of modern soil in the catchment that predominantly formed on limestone bedrock. Key findings are: (1) modern soils and sediments of Core-HQ both contain a mixture of magnetite (Mt), maghemite (Mgh) and hematite (Ht), but magnetic concentration of the soils is one order higher; (2) a superparamagnetic (SP) fraction of Mt/Mgh dominates in the soils whereas in Core-HQ the SP contribution is generally very low; (3) a larger grain-size fraction of Mt/Mgh and Ht is also present in the soils. We explain variations of magnetic concentration and CC in Core-HQ by an increased wind transport of soil and a decreased surface water transport of carbonate and soil in less humid periods. Low-temperature oxidation of magnetite in the catchment is as another crucial process that reflects weathering conditions and is likely sensitive to humidity; the degree of LTO can be semi-quantified by the magnetic parameters ARM/SIRM and S -ratio. Combining CC, ARM/SIRM and S -ratio values, we derive a weathering intensity (WI). The WI index variation along Core-HQ shows strong fluctuations on a 100-kyr eccentricity scale in the lower part, especially during ~630–380 ka, followed by a long period (~320–80 ka) of persistently weaker weathering (drier?) conditions with low variability, and a rapid return of much stronger weathering (wetter?) conditions at ~80 ka. We suggest that a reduced influence of the Indian summer monsoon accounts for less moisture supply to the region and lower climatic variability.

Description: Magnetic mineralogy of crustal rocks has important implications for understanding continental crustal evolution and origin of regional magnetic anomalies. However, magnetic properties of the deep continental crust are still poorly understood. In this paper, measurements of density ( ), mass-specific magnetic susceptibility ( ), natural remanent magnetization (NRM) and magnetic hysteresis loops, temperature-dependent magnetic susceptibility ( – T ), chemical and mineral analyses were conducted on Archean gneiss samples from the Jixian petrophysical section in the Precambrian terrain, northeastern North China Craton, with the aim of refining understanding of magnetic phase transformations in the deep crustal rocks. Results show that density and rock magnetic properties change distinctly with metamorphic facies. The dominant magnetic mineral is magnetite, while a little hematite is present in a few samples. Together with geochemical and mineralogical compositions, it is inferred that progressive increase in metamorphic grade from east to west is the major cause for magnetic enhancement of the lower crust in the studied section. Therefore, we conclude that study of magnetic phases of deep crustal rocks can offer important insights into the history of high metamorphic grade terranes.

Description: The C -response is a conventional transfer function in global electromagnetic induction research and is traditionally determined from observations of magnetic variations in the vertical and horizontal components. Its interpretation relies on the assumption that the source of the variations is well approximated by a large-scale symmetric (magnetospheric) ring current, described by a single spherical harmonic. However, there is growing evidence for a more complex structure of this source. In this paper, we investigate the variability of C -responses due to sources different from the dominating large-scale symmetric ring current. We show that the effect is significant and persists at all periods. Describing the magnetospheric source by a single spherical harmonic coefficient thus injects substantial errors into the estimated responses. To overcome the problem, we introduce arrays of alternative transfer functions that relate the components of the magnetic variation to different spherical harmonic coefficients. These transfer functions can handle a complex spatial structure of the magnetospheric source. Compared to C -responses, we observe a significant increase in the coherencies relating input and output quantities of the new transfer functions, especially at high latitudes. This increases the usability of observatory magnetic data for the recovery of global 3-D mantle conductivity structure.

Description: Magnetic mineralogy of crustal rocks has important implications for understanding continental crustal evolution and origin of regional magnetic anomalies. However, magnetic properties of the deep continental crust are still poorly understood. In this paper, measurements of density ( ), mass-specific magnetic susceptibility ( ), natural remanent magnetization (NRM) and magnetic hysteresis loops, temperature-dependent magnetic susceptibility ( – T ), chemical and mineral analyses were conducted on Archean gneiss samples from the Jixian petrophysical section in the Precambrian terrain, northeastern North China Craton, with the aim of refining understanding of magnetic phase transformations in the deep crustal rocks. Results show that density and rock magnetic properties change distinctly with metamorphic facies. The dominant magnetic mineral is magnetite, while a little hematite is present in a few samples. Together with geochemical and mineralogical compositions, it is inferred that progressive increase in metamorphic grade from east to west is the major cause for magnetic enhancement of the lower crust in the studied section. Therefore, we conclude that study of magnetic phases of deep crustal rocks can offer important insights into the history of high metamorphic grade terranes.