ALBERT

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Description: In this work, we carry out numerical simulations of the seismo-electromagnetic fields associated with a fault in a porous medium by considering the electrokinetic effect. In addition to porous materials, the adopted layered models comprise solid materials in which the electrokinetic effect is inoperative. First, sensitivity study is performed for the evanescent and direct radiation electromagnetic (EM) waves generated by a double couple point source embedded in a porous half-space below a solid half-space. Results suggest that both the evanescent and direct radiation EM waves are sensitive to some medium properties, for example porosity, salinity, fluid viscosity, and conductivity of solid layer. Then, adopting an eight-layer half-space model, we simulate the seismic and EM wavefields generated by the rupture process of a finite fault. It is shown that the electrokinetic effect is able to generate observable corupture EM signals, but the observability depends on some factors such as the epicentral distance, properties of the medium where the fault is located, and local activity levels of EM noise. Time synchronization coseismic EM signals are recorded when the receiver is close to the ground water level but located in a solid medium. In addition to the post-seismic electric field, our results also show the post-seismic magnetic field which has not been identified in previous simulation studies on the electrokinetic effect. The generation of the post-seismic magnetic field probably requires a sufficiently strong medium heterogeneity or fluid-pressure gradient.

Description: Fault ruptures in the Earth's crust generate both elastic and electromagnetic (EM) waves. If the corresponding EM signals can be observed, then earthquakes could be detected before the first seismic waves arrive. In this study, I consider the piezomagnetic effect as a mechanism that converts elastic waves to EM energy, and I derive analytical formulas for the conversion process. The situation considered in this study is a whole-space model, in which elastic and EM properties are uniform and isotropic. In this situation, the governing equations of the elastic and EM fields, combined with the piezomagnetic constitutive law, can be solved analytically in the time domain by ignoring the displacement current term. Using the derived formulas, numerical examples are investigated, and the corresponding characteristics of the expected magnetic signals are resolved. I show that temporal variations in the magnetic field depend strongly on the electrical conductivity of the medium, meaning that precise detection of signals generated by the piezomagnetic effect is generally difficult. Expected amplitudes of piezomagnetic signals are estimated to be no larger than 0.3 nT for earthquakes with a moment magnitude of ≥7.0 at a source distance of 25 km; however, this conclusion may not extend to the detection of real earthquakes, because piezomagnetic stress sensitivity is currently poorly constrained.

Description: This paper presents a methodology to sample equivalence domain (ED) in nonlinear partial differential equation (PDE)-constrained inverse problems. For this purpose, we first applied state-of-the-art stochastic optimization algorithm called Covariance Matrix Adaptation Evolution Strategy (CMAES) to identify low-misfit regions of the model space. These regions were then randomly sampled to create an ensemble of equivalent models and quantify uncertainty. CMAES is aimed at exploring model space globally and is robust on very ill-conditioned problems. We show that the number of iterations required to converge grows at a moderate rate with respect to number of unknowns and the algorithm is embarrassingly parallel. We formulated the problem by using the generalized Gaussian distribution. This enabled us to seamlessly use arbitrary norms for residual and regularization terms. We show that various regularization norms facilitate studying different classes of equivalent solutions. We further show how performance of the standard Metropolis–Hastings Markov chain Monte Carlo algorithm can be substantially improved by using information CMAES provides. This methodology was tested by using individual and joint inversions of magneotelluric, controlled-source electromagnetic (EM) and global EM induction data.

Description: Recent high-resolution aeromagnetic surveys in South Norway have revealed numerous remanent anomalies over Mesoproterozoic metamorphic rocks. Studies on the nature of the minerals that are the remanent carriers has led to discoveries of titanohematite samples with unusual magnetic properties caused by nanoscale exsolution lamellae with their related lamellar magnetism. Here we focus on a rock unit dominated by quartz-plagioclase-biotite granulite containing titanohematite grains with a strong lattice-preferred orientation parallel to regional foliation. When samples with their natural remanent magnetization (NRM), acquired nearly 1 billion years ago, are cooled to 10 K and hysteresis loops measured, these loops show bi-modal exchange bias caused by the magnetism induced within the ilmenite by antiferromagnetic coupling with the adjacent lamellar NRM. By contrast when the samples are cooled in a strong magnetic field (1.5 Tesla), this results in unimodal lamellar magnetism, and, below the T N of ilmenite it adopts a consistent negative orientation, giving rise to unimodal negative exchange bias of 〉500 mT. The results presented here cover the chemical and magnetic properties, Mossbauer results and transmission electron microscopy of the titanohematite and ilmenite lamellae. Initial magnetic experiments indicated the shifts found in the exchange-bias experiments were directly related to the orientation of the sample to the applied field and the initial state of the NRM. In most samples with these unusual magnetic properties, ilmenite lamellae could not be seen in an optical or a scanning electron microscope. However magnetic experiments gave proof of the presence of ilmenite, later confirmed by Mössbauer spectroscopy. Several attempts were made to identify ilmenite in TEM studies, finally successful in showing ilmenite lamellae parallel to (001) of hematite with thicknesses ~1.2 to 1.7 nm and aspect ratios 7–13. Here we compare new TEM images and the magnetic behaviour of these samples to the MOD2 samples that previously showed extraordinary exchange bias properties, and investigate further the nature of these magnetic minerals.

Description: Rock-magnetic and geochemical characteristics of three Vertisol profiles with different degree of textural differentiation have been studied. Thermomagnetic analyses, thermal demagnetization of laboratory remanences and acquisition of isothermal remanence curves are applied for identification of iron oxide mineralogy. The main magnetic minerals in Vertisols are ferrihydrite, single-domain magnetite, maghemite and hematite. Variations in magnetic susceptibility, anhysteretic remanent magnetization, isothermal remanent magnetization, as well as different ratios (Xarm/X, ARM/SIRM, S-ratio) along depth are studied. Concentration of magnetic minerals in Vertisols is low, influenced by the intense reductomorphic processes. The lowest magnetic susceptibility is found in the most texturally differentiated soil. However, rock-magnetic data suggest the presence of small, but well defined fraction of single domain-like magnetite with relatively wide grain-size distribution found in those parts of the profiles, which are subjected to most intense and frequent seasonal changes in oxidation-reduction conditions. It is suggested that this fraction is formed as a result of transformations of ferrihydrite under repeated cycles of anaerobic/aerobic conditions. Based on geochemical data, CALMAG weathering index was calculated for the three Vertisols. Using the established relation between CALMAG and mean annual precipitation (MAP), palaeo-MAP was evaluated for the studied profiles. The obtained MAP estimations fall in the range 1000–1200 mm and are much higher compared to contemporary precipitation in the area (MAP in the interval 540–770 mm). This finding confirms the relict character of Vertisols on Bulgarian territory and gives more information about the palaeoclimate during the initial stages of Vertisol formation.

Description: This work presents MARE2DEM, a freely available code for 2-D anisotropic inversion of magnetotelluric (MT) data and frequency-domain controlled-source electromagnetic (CSEM) data from onshore and offshore surveys. MARE2DEM parametrizes the inverse model using a grid of arbitrarily shaped polygons, where unstructured triangular or quadrilateral grids are typically used due to their ease of construction. Unstructured grids provide significantly more geometric flexibility and parameter efficiency than the structured rectangular grids commonly used by most other inversion codes. Transmitter and receiver components located on topographic slopes can be tilted parallel to the boundary so that the simulated electromagnetic fields accurately reproduce the real survey geometry. The forward solution is implemented with a goal-oriented adaptive finite-element method that automatically generates and refines unstructured triangular element grids that conform to the inversion parameter grid, ensuring accurate responses as the model conductivity changes. This dual-grid approach is significantly more efficient than the conventional use of a single grid for both the forward and inverse meshes since the more detailed finite-element meshes required for accurate responses do not increase the memory requirements of the inverse problem. Forward solutions are computed in parallel with a highly efficient scaling by partitioning the data into smaller independent modeling tasks consisting of subsets of the input frequencies, transmitters and receivers. Non-linear inversion is carried out with a new Occam inversion approach that requires fewer forward calls. Dense matrix operations are optimized for memory and parallel scalability using the ScaLAPACK parallel library. Free parameters can be bounded using a new non-linear transformation that leaves the transformed parameters nearly the same as the original parameters within the bounds, thereby reducing non-linear smoothing effects. Data balancing normalization weights for the joint inversion of two or more data sets encourages the inversion to fit each data type equally well. A synthetic joint inversion of marine CSEM and MT data illustrates the algorithm's performance and parallel scaling on up to 480 processing cores. CSEM inversion of data from the Middle America Trench offshore Nicaragua demonstrates a real world application. The source code and MATLAB interface tools are freely available at http://mare2dem.ucsd.edu .

Description: The spectral complex conductivity of a water-bearing sand during interaction with carbon dioxide (CO 2 ) is influenced by multiple, simultaneous processes. These processes include partial saturation due to the replacement of conductive pore water with CO 2 and chemical interaction of the reactive CO 2 with the bulk fluid and the grain-water interface. We present a laboratory study on the spectral induced polarization of water-bearing sands during exposure to and flow-through by CO 2 . Conductivity spectra were measured successfully at pressures up to 30 MPa and 80 °C during active flow and at steady-state conditions concentrating on the frequency range between 0.0014 and 100 Hz. The frequency range between 0.1 and 100 Hz turned out to be most indicative for potential monitoring applications. The presented data show that the impact of CO 2 on the electrolytic conductivity may be covered by a model for pore-water conductivity, which depends on salinity, pressure and temperature and has been derived from earlier investigations of the pore-water phase. The new data covering the three-phase system CO 2 -brine-sand further show that chemical interaction causes a reduction of surface conductivity by almost 20 per cent, which could be related to the low pH-value in the acidic environment due to CO 2 dissolution and the dissociation of carbonic acid. The quantification of the total CO 2 effect may be used as a correction during monitoring of a sequestration in terms of saturation. We show that this leads to a correct reconstruction of fluid saturation from electrical measurements. In addition, an indicator for changes of the inner surface area, which is related to mineral dissolution or precipitation processes, can be computed from the imaginary part of conductivity. The low frequency range between 0.0014 and 0.1 Hz shows additional characteristics, which deviate from the behaviour at higher frequencies. A Debye decomposition approach is applied to isolate the feature dominating the data at low frequencies. We conclude from our study that electrical conductivity is not only a highly sensitive indicator for CO 2 saturation in pore space. When it is measured in its full spectral and complex form it contains additional information on the chemical state of the system, which holds the potential of getting access to both saturation and interface properties with one monitoring method.

Description: Our study attempts to identify a characteristic magnetic signature of overbank sediments exhibiting anthropogenically induced magnetic enhancement and thereby to distinguish them from unenhanced sediments with weak magnetic background values, using a novel approach based on data mining methods, thus providing a mean of rapid pollution determination. Data were obtained from 539 bulk samples from vertical profiles through overbank sediment, collected on seven rivers in the eastern Czech Republic and three rivers in northwest England. k -Means clustering and hierarchical clustering methods, paired group (UPGMA) and Ward's method, were used to divide the samples to natural groups according to their attributes. Interparametric ratios: SIRM/; SIRM/ARM; and S –0.1T were chosen as attributes for analyses making the resultant model more widely applicable as magnetic concentration values can differ by two orders. Division into three clusters appeared to be optimal and corresponded to inherent clusters in the data scatter. Clustering managed to separate samples with relatively weak anthropogenically induced enhancement, relatively strong anthropogenically induced enhancement and samples lacking enhancement. To describe the clusters explicitly and thus obtain a discrete magnetic signature, classification rules (JRip method) and decision trees (J4.8 and Simple Cart methods) were used. Samples lacking anthropogenic enhancement typically exhibited an S –0.1T 〈 c. 0.5, SIRM/ARM 〈 c. 150 and SIRM/ 〈 c. 6000 A m –1 . Samples with magnetic enhancement all exhibited an S –0.1T 〉 0.5. Samples with relatively stronger anthropogenic enhancement were unequivocally distinguished from the samples with weaker enhancement by an SIRM/ARM 〉 c. 150. Samples with SIRM/ARM in a range c. 126–150 were classified as relatively strongly enhanced when their SIRM/ 〉 18 000 A m –1 and relatively less enhanced when their SIRM/ 〈 18 000 A m –1 . An additional rule was arbitrary added to exclude samples with fd% 〉 6 per cent from anthropogenically enhanced clusters as samples with natural magnetic enhancement. The characteristics of the clusters resulted mainly from the relationship between SIRM/ARM and the S –0.1T , and SIRM/ and the S –0.1T . Both SIRM/ARM and SIRM/ increase with increasing S –0.1T values reflecting a greater level of anthropogenic magnetic particles. Overall, data mining methods demonstrated good potential for utilization in environmental magnetism.

Description: The construction of geomagnetic, archaeomagnetic or palaeomagnetic field models requires some prior knowledge about the actual field, which can be gathered from the statistical properties of the field over a variety of length-scales and timescales. However, available geomagnetic data on centennial to millennial periods are too sparse to infer directly these statistical properties. We thus use high-resolution numerical simulations of the geodynamo to test a method for estimating the temporal power spectra (or equivalently the autocovariance functions) of the individual Gauss coefficients that describe the geomagnetic field outside the Earth's fluid outer core. Based on the spectral analysis of our simulations, we argue that a prior for the observational geomagnetic field over decennial to millennial periods can be constructed from the statistics of the field during the short satellite era. The method rests on the assumption that time-series of spherical harmonic coefficients can be considered as realizations of stationary and differentiable stochastic processes, namely order 2 autoregressive (AR2) processes. In the framework of these processes, the statistics of Gauss coefficients are well constrained by their variance and one or two timescales. We find that the time spectra in the dynamo simulations of all Gauss coefficients but the axial dipole are well approximated by the spectra of AR2 processes characterized by only one timescale. The process parameters can simply be deduced from instantaneous estimates of the spatial power spectra of the magnetic field and of its first time derivative. Some deviations of the Gauss coefficients statistics from this minimal model are also discussed. Characterizing the axial dipole clearly requires a more sophisticated AR2 process, with a second distinct timescale.

Description: The Luonan Basin, located in the transitional zone between temperate and subtropical China, is an important locality for human evolution during the early to middle Pleistocene. The loess deposits in the Luonan Basin contain numerous in situ lithic artefacts; the deposits also constitute suitable material for dating the artefacts and are potentially useful for reconstructing the climatic fluctuations which is important for studying the adaptation and occupation of the area by early humans. We carried out a combined rock magnetic and geochemical investigation of a loess sequence from the Liuwan Palaeolithic site in the Luonan Basin. The results indicate a mixture of magnetic minerals, including magnetite/maghemite and hematite/goethite. Magnetic susceptibility was used as a palaeoclimate proxy on the Chinese Loess Plateau; however, its application to the Luonan Basin may be problematic because the provenance of the loess parent material, as well as the depositional environment, differs from that of the Chinese Loess Plateau. We found that rock magnetic parameters related to the grain size of magnetic minerals, such as SIRM/ and ARM /SIRM, are better palaeoclimatic indicators than magnetic susceptibility. Overall, the magnetic results, together with the results of bulk grain-size and chemical index of alteration, indicate that the interglacial environment of early humans in Luonan Basin was warmer and more humid than the coeval environment of the Chinese Loess Plateau.

Description: We develop the theory for the magnetic to magnetic and kinetic to magnetic energy transfer between different spherical harmonic degrees due to the interaction of fluid flow and radial magnetic field at the top of the Earth's core. We show that non-zero secular variation of the total magnetic energy could be significant and may provide evidence for the existence of stretching secular variation, which suggests the existence of radial motions at the top of the Earth's core—whole core convection or MAC waves. However, the uncertainties of the small scales of the geomagnetic field prevent a definite conclusion. Combining core field and flow models we calculate the detailed magnetic to magnetic and kinetic to magnetic energy transfer matrices. The magnetic to magnetic energy transfer shows a complex behaviour with local and non-local transfers. The spectra of magnetic to magnetic energy transfers show clear maxima and minima, suggesting an energy cascade. The kinetic to magnetic energy transfers, which are much weaker due to the weak poloidal flow, are either local or non-local between degree one and higher degrees. The patterns observed in the matrices resemble energy transfer patterns that are typically found in 3-D MHD numerical simulations.

Description: The magnetic susceptibility measured in alternating field can in general be resolved into a component that is in-phase with the applied field and a component that is out-of-phase. While in non-conductive diamagnetic, paramagnetic and many ferromagnetic materials the phase is effectively zero, in some ferromagnetic minerals, such as pyrrhotite, hematite, titanomagnetite or small magnetically viscous grains of magnetite, it is clearly non-zero. The anisotropy of out-of-phase susceptibility (opAMS) can then be used as a tool for the direct determination of the magnetic subfabrics of the minerals with non-zero phase. The error in determination of out-of-phase susceptibility non-linearly increases with decreasing phase angle. This may result in imprecise determination of the opAMS in specimens with very low phase angle. The degree of opAMS is higher than that of ipAMS, which may in contrast result in slightly increasing precision n the opAMS determination. It is highly recommended to inspect the results of the statistical tests of each specimen and to exclude the specimens whose opAMS is determined with insufficient precision from further processing. In rocks, whose magnetism is dominated by the mineral with non-zero out-of-phase susceptibility, the principal directions of the opAMS and ipAMS are virtually coaxial, while the degree of opAMS is higher than that of ipAMS. In some specific cases, the opAMS provides us with similar data to those provided by anisotropies of low-field dependent susceptibility and frequency-dependent susceptibility. The advantage of the opAMS compared to the other two anisotropies is its simultaneous measurement with the ipAMS during one measuring process, while the other two anisotropies require the AMS measurements in several fields or at least at two operating frequencies.

Description: The temperature dependence of magnetic susceptibility ( –T curves) and magnetization ( M–T curves) has been used as a routine rock magnetic tool to characterize the magnetic mineralogy and magnetic granulometry of Chinese loess/palaeosols. However, palaeoclimatic interpretation of these thermomagnetic analyses remains controversial. In the present study, total organic carbon (TOC), thermomagnetic and low-temperature magnetic experiments on Mid-Late Pleistocene loess/palaeosols in Central Asia and the Chinese Loess Plateau (CLP) have been conducted. We found that the M ( T ) cooling curves at room temperature were mostly lower than the corresponding heating curves, whereas for the ( T ) analyses the cooling curves at room temperature were always much higher than the heating curves. Low-temperature magnetic measurements demonstrated that a large amount of superparamagnetic ferrimagnetic particles were produced during the thermal treatment and resulted in the aforementioned differences. This finding further indicated that the use of the M–T curves to estimate the relative content of maghemite in the loess/palaeosols from the CLP was problematic. In addition, a positive correlation exists between the TOC and the frequency-dependent susceptibility ( FD ) in the CLP, suggesting that stronger pedogenesis would result in the simultaneous increase in the content of both maghemite and organic matter. Consequently, the parameters 1 (representing the relative content of pedogenic maghemite), 2 ([ ph – ] + 1 ) and ph (related to the organic matter concentration), which can be calculated from the –T analyses, can potentially be used as new indicators of pedogenesis and palaeoclimate in Central Asia and the CLP.

Description: This is the second of three papers investigating properties of titanohematite-bearing quartzo-feldspathic rocks that create a significant remanent magnetic anomaly in the Modum District, South Norway. The first paper provided initial magnetic results, mineralogical characterization and evidence for the presence of lamellar magnetism. In this paper, knowledge of lamellar magnetic properties is explored through experiments where ilmenite lamellae were magnetized below 57 K, and interact magnetically along interfaces with the titanohematite host. Samples with known NRM directions were placed in specific orientations in an MPMS then cooled in zero field to 5 K, where hysteresis loops were measured in fields up to 5 Tesla. This assured that results were ultimately related to the natural lamellar magnetism produced during cooling ~1 billion years ago. In a second set of experiments the same oriented samples, were subjected to a +5 Tesla (T) field then field cooled to 5 K before hysteresis experiments. The first experiments consistently produced asymmetric shifted hysteresis loops with two loop separations, one in a positive field and one in a negative field. Without exception, when the NRM was oriented toward the negative field end of the MPMS, the bimodal loop showed a dominant loop separation in a positive field. By contrast, when the NRM was oriented toward the positive field end of the MPMS, the bimodal loop showed a dominant loop separation in a negative field. Both observations are consistent with antiferromagnetic coupling between the hard magnetization of ilmenite and the more easily shifted lamellar magnetism of the hematite. The bimodal nature of the loops indicates that the NRMs are vector sums of natural lamellar moments, which are oriented both positively and negatively, and that these opposite moments control the orientations of ilmenite magnetizations when cooling through 57 K. Here, extreme exchange biases up to 1.68 T were measured. The second set of experiments produced asymmetric shifted hysteresis loops with one opening always in the negative field. These observations indicate that the +5 T field applied at room temperature rotated the hematite lamellar magnetism in a positive direction, so that upon cooling all the ilmenite lamellae acquired negative magnetic moments, thus causing unimodal negatively shifted loops. Here, the largest exchange bias among the unimodal loops was only 0.7 T. These results will be used in paper III to build a better understanding of lamellar magnetism at the atomic layer scale.

Description: Non-heating palaeointensity methods are a vital tool to explore magnetic field strength variations recorded by thermally sensitive materials of both terrestrial and extraterrestrial origin. One such method is the calibrated pseudo-Thellier method in which a specimen's natural remanent magnetization is alternating field demagnetized and replaced with a laboratory induced anhysteretic remanent magnetization (as an analogue of a thermoremanent magnetization, TRM). Using a set of 56 volcanic specimens given laboratory TRMs in fields of 10–130 μT, we refine the calibration of the pseudo-Thellier method and better define the uncertainty associated with its palaeointensity estimates. Our new calibration, obtained from 32 selected specimens, resolves the issue of non-zero intercept, which is theoretically predicted, but not satisfied by any previous calibration. The range of individual specimen calibration factors, however, is relatively large, but consistent with the variability expected for SD magnetite. We explore a number of rock magnetic parameters in an attempt to identify selection thresholds for reducing the calibration scatter, but fail to find a suitable choice. We infer that our careful selection process, which incorporates more statistics then previous studies, may be largely screening out any strong rock magnetic dependence. Some subtle grain size or mineralogical dependencies, however, remain after selection, but cannot be discerned from the scatter expected for grain size variability of SD magnetite. As a consequence of the variability in the calibration factor, the uncertainty associated with pseudo-Thellier results is much larger than previously indicated. The scatter of the calibration is ~25 per cent of the mean value, which implies that, when combined with the scatter of results typically obtained from a single site, the uncertainty of averaged pseudo-Thellier results will always be 〉25 per cent. As such, pseudo-Thellier results should be complementary to, and cross-validated with results from other methods. Nevertheless, the pseudo-Thellier method remains a valuable tool for obtaining palaeointensity estimates from thermally sensitive terrestrial and extraterrestrial materials and with careful data selection and analysis can yield results that are accurate to within a factor of 4 or better.

Description: As a consequence of measuring time variations of the electric and the magnetic field, which are related to current flow and charge distribution, magnetotelluric (MT) data in 2-D and 3-D environments are not only sensitive to the geoelectrical structures below the measuring points but also to any lateral anomalies surrounding the acquisition site. This behaviour complicates the characterization of the electrical resistivity distribution of the subsurface, particularly in complex areas. In this manuscript we assess the main advantages of complementing the standard MT impedance tensor ( Z ) data with interstation horizontal magnetic tensor ( H ) and geomagnetic transfer function ( T ) data in constraining the subsurface in a 3-D environment beneath a MT profile. Our analysis was performed using synthetic responses with added normally distributed and scattered random noise. The sensitivity of each type of data to different resistivity anomalies was evaluated, showing that the degree to which each site and each period is affected by the same anomaly depends on the type of data. A dimensionality analysis, using Z , H and T data, identified the presence of the 3-D anomalies close to the profile, suggesting a 3-D approach for recovering the electrical resistivity values of the subsurface. Finally, the capacity for recovering the geoelectrical structures of the subsurface was evaluated by performing joint inversion using different data combinations, quantifying the differences between the true synthetic model and the models from inversion process. Four main improvements were observed when performing joint inversion of Z , H and T data: (1) superior precision and accuracy at characterizing the electrical resistivity values of the anomalies below and outside the profile; (2) the potential to recover high electrical resistivity anomalies that are poorly recovered using Z data alone; (3) improvement in the characterization of the bottom and lateral boundaries of the anomalies with low electrical resistivity; and (4) superior imaging of the horizontal continuity of structures with low electrical resistivity. These advantages offer new opportunities for the MT method by making the results from a MT profile in a 3-D environment more convincing, supporting the possibility of high-resolution studies in 3-D areas without expending a large amount of economical and computational resources, and also offering better resolution of targets with high electrical resistivity.

Description: Inference of globally averaged eustatic sea level (ESL) rise since the Last Glacial Maximum (LGM) highly depends on the interpretation of relative sea level (RSL) observations at Barbados and Bonaparte Gulf, Australia, which are sensitive to the viscosity structure of Earth's mantle. Here we examine the RSL changes at the LGM for Barbados and Bonaparte Gulf ( ${\rm{RSL}}_{\rm{L}}^{{\rm{Bar}}}$ and ${\rm{RSL}}_{\rm{L}}^{{\rm{Bon}}}$ ), differential RSL for both sites ( $\Delta {\rm{RSL}}_{\rm{L}}^{{\rm{Bar}},{\rm{Bon}}})$ and rate of change of degree-two harmonics of Earth's geopotential due to glacial isostatic adjustment (GIA) process (GIA-induced 2 ) to infer the ESL component and viscosity structure of Earth's mantle. Differential RSL, $\Delta {\rm{RSL}}_{\rm{L}}^{{\rm{Bar}},{\rm{Bon}}}$ and GIA-induced 2 are dominantly sensitive to the lower-mantle viscosity, and nearly insensitive to the upper-mantle rheological structure and GIA ice models with an ESL component of about (120–130) m. The comparison between the predicted and observationally derived $\Delta {\rm{RSL}}_{\rm{L}}^{{\rm{Bar}},{\rm{Bon}}}$ indicates the lower-mantle viscosity higher than ~2  x  10 22 Pa s, and the observationally derived GIA-induced 2 of –(6.0–6.5)  x  10 –11 yr –1 indicates two permissible solutions for the lower mantle, ~10 22 and (5–10)  x  10 22 Pa s. That is, the effective lower-mantle viscosity inferred from these two observational constraints is (5–10)  x  10 22 Pa s. The LGM RSL changes at both sites, ${\rm{RSL}}_{\rm{L}}^{{\rm{Bar}}}$ and ${\rm{RSL}}_{\rm{L}}^{{\rm{Bon}}}$ , are also sensitive to the ESL component and upper-mantle viscosity as well as the lower-mantle viscosity. The permissible upper-mantle viscosity increases with decreasing ESL component due to the sensitivity of the LGM sea level at Bonaparte Gulf ( ${\rm{RSL}}_{\rm{L}}^{{\rm{Bon}}}$ ) to the upper-mantle viscosity, and inferred upper-mantle viscosity for adopted lithospheric thicknesses of 65 and 100 km is (1–3)  x  10 20 Pa s for ESL~130 m and (4–10)  x  10 20 Pa s for ESL~125 m. The former solution of (1–3)  x  10 20 Pa s is consistent with the inferences from the postglacial differential RSL changes in the Australian region and also inversion study of far-field sea-level data. The inference of the viscosity structure based on these four observational constraints is, however, relatively insensitive to the viscosity structure of D'' layer.

Description: We present the theory, algorithms and implementation of a parallel finite-volume algorithm for the solution of the incompressible magnetohydrodynamic (MHD) equations using unstructured grids that are applicable for a wide variety of geometries. Our method implements a mixed Adams-Bashforth/Crank-Nicolson scheme for the nonlinear terms in the MHD equations and we prove that it is stable independent of the time step. To ensure that the solenoidal condition is met for the magnetic field, we use a method whereby a pseudo-pressure is introduced into the induction equation; since we are concerned with incompressible flows, the resulting Poisson equation for the pseudo-pressure is solved alongside the equivalent Poisson problem for the velocity field. We validate our code in a variety of geometries including periodic boxes, spheres, spherical shells, spheroids and ellipsoids; for the finite geometries we implement the so-called ferromagnetic or pseudo-vacuum boundary conditions appropriate for a surrounding medium with infinite magnetic permeability. This implies that the magnetic field must be purely perpendicular to the boundary. We present a number of comparisons against previous results and against analytical solutions, which verify the code's accuracy. This documents the code's reliability as a prelude to its use in more difficult problems. We finally present a new simple drifting solution for thermal convection in a spherical shell that successfully sustains a magnetic field of simple geometry. By dint of its rapid stabilization from the given initial conditions, we deem it suitable as a benchmark against which other self-consistent dynamo codes can be tested.

Description: We consider an inviscid and perfectly conducting fluid sphere in rapid rotation and permeated by a background magnetic field. Such a system admits normal modes in the form of torsional oscillations, namely azimuthal motions of cylinders coaxial with the rotation axis. We analyse this system for a particular background magnetic field that provides a new closed form normal mode solution. We derive Wentzel–Kramers–Brillouin–Jeffreys (WKBJ) approximations to the normal modes, and focus particularly on the reflections that take place on the rotation axis and at the equator. We propose a procedure to calculate the reflection coefficients and we discuss the analogy of our findings with well-known seismological results. Our analytical results are tested against numerical calculations and show good agreement.

Description: The 69 Ma Prospector Mountain stock is located in southwestern Yukon in the northern Canadian Cordillera. This massive monzonite-syenogranite stock is thought to be the intrusive volcanic centre for the surrounding coeval Carmacks Group volcanics. Anomalous palaeomagnetic data from these volcanics have provided the only evidence for the commonly posited hypothesis that the Yukon-Tanana terrane (YTT) was part of the far-travelled (1950 ± 600 km northward) ‘Baja BC’ terrane from 70 to 50 Ma. All other geologic evidence and averaged palaeomagnetic data support a northward displacement of ~415 ± 15 km. This study provides a direct test of the Carmacks volcanics’ estimate and examines the possible causes of its anomalous results. Both the stock and volcanics are unmetamorphosed and rest unconformably on metamorphosed basement rocks of the YTT. Palaeomagnetic and mineral magnetic results from 17 of 19 tested sites (218 specimens) in the Prospector stock and its peripheral skarn isolated a stable thermoremanent magnetization (TRM) in magnetite or low-Ti titanomagnetite that was mostly determined on demagnetization between temperatures of 500 and 580 °C. The TRM has a direction of Decl. = 8.3°, Incl. = 82.4° ( N = 17, k = 71.9, α 95 = 4.2°), providing a non-significant northwards translation estimate of 70 ± 880 km for the YTT. The normal-polarity TRM direction at Prospector Mountain provides a highly significant palaeomagnetic reversals test with the reversed-polarity TRM of the 70 Ma Swede Dome stock, another volcanic centre of the Carmacks Group about 190 km to the north. The test affirms to a high probability that both stocks carry primary TRMs and have not been tectonically tilted significantly since emplacement. Combining the palaeopoles for these two stocks with that for the 75 Ma Mount Lorne volcanic centre stock about 210 km south of Prospector Mountain yields a combined northward translation estimate of 330 ± 400 km for the YTT since ~71 Ma. This estimate agrees closely with the 415 ± 15 km estimate from geological constraints on strike-slip motion along the inboard Tintina fault zone since ~70 Ma, and with the palaeomagnetic estimate of 650 ± 450 km northward since ~108 Ma from four mid-Cretaceous batholiths. In contrast, the palaeomagnetic directions from four areas of the ~70 Ma Carmacks volcanics are poorly clustered. Nonetheless their estimates of ~1900 km northward have been used for a quarter of a century to support a far-travelled ‘Baja BC’ tectonic model for the YTT. This paper discusses and concludes that the anomalous far-travelled estimates from the Carmacks volcanics are attributable to the unfortunate additive effects of inadequate averaging of secular variation, dipole offset error, unrecognized primary dip and other possible causes.

Description: Previous magnetotelluric (MT) studies of the high-temperature Coso geothermal system in California identified a subvertical feature of low resistivity (2–5 Ohm m) and appreciable lateral extent (〉1 km) in the producing zone of the East Flank field. However, these models could not reproduce gross 3-D effects in the recorded data. We perform 3-D full-tensor inversion and retrieve a resistivity model that out-performs previous 2-D and 3-D off-diagonal models in terms of its fit to the complete 3-D MT data set as well as the degree of modelling bias. Inclusion of secondary Z xx and Z yy data components leads to a robust east-dip (60) to the previously identified conductive East Flank reservoir feature, which correlates strongly with recently mapped surface faults, downhole well temperatures, 3-D seismic reflection data, and local microseismicity. We perform synthetic forward modelling to test the best-fit dip of this conductor using the response at a nearby MT station. We interpret the dipping conductor as a fractured and fluidized compartment, which is structurally controlled by an unmapped blind East Flank fault zone.

Description: The island of Terceira (38.5°N, 27.2°W) in the Azores Archipelago boasts a wealth of well-dated recent volcanic products; it is therefore the ideal location to improve the resolution of the full-vector geomagnetic field record for the Mid-Atlantic Ocean. We investigated 21 (sub-)sites sampled from 10 different cooling units; six of these cooling units are from the Holocene, the other four are approximately coeval with the age of known geomagnetic excursions. After successfully determining the palaeodirections from all but one site, we applied the multimethod palaeointensity approach to all sites. Here, we used three different palaeointensity methods: the IZZI-Thellier protocol, the Multispecimen technique, and the calibrated pseudo-Thellier method. For the Holocene sites four robust palaeointensity estimates were obtained, they mostly agree with data from the nearby island of São Miguel and data from the Iberian Peninsula and West Africa. Two older sites, dated at 25.7 and 46.3 ka, produced low palaeointensities: 24.4 and 24.1 μT, respectively. As both the IZZI-Thellier and pseudo-Thellier methods yielded technically acceptable results for these sites we used them to improve the pseudo-Thellier calibration relation and establish its applicability to lower palaeointensities. Sixty per cent of all sampled cooling units/ages in our study yielded a technically successful palaeointensity estimate; this unusually high success rate is attributed both to our sampling strategy and the implementation of the multimethod approach.

Description: Two lava flows from the Ceboruco volcano in west-central Mexico were sampled for palaeomagnetic dating. The younger one was emitted in 1870 and used to validate the method, while the older one known as Ceboruco flow is of unknown age but probably younger than ~1005 AD and older than 1528 AD. Each flow was sampled in at least four sites, in order to unravel between site variations. For the 1870 flow, between site differences were notable and additionally post-cooling block movements were important; therefore, two sites had to be rejected. Three sites from the vent area and one at the tip of the 1870 flow provided well-constrained directions. This is also true for Ceboruco lava flow, and overall mean directions and palaeointensities were then used for palaeomagnetic dating applying the Matlab tool archaeo_dating and the global palaeosecular variation model SHA.DIF.14k. For the 1870 lava flow, the dating resulted in an age ranging between 1755 and 1871 AD (95 per cent probability level), which includes the real emplacement age. In addition, the Ceboruco lava flow was dated between 1000 and 1134 AD, which is close to the large plinian Jala eruption producing the crater of Ceboruco volcano around 1005 AD. This age is older than previously assumed and suggests an emplacement only shortly after the Jala eruption. As this lava flow is considered to be the youngest one of seven post-Jala lava flows, the age also defines a period of inactivity of Ceboruco volcano of about 730–860 yr before the historic 1870 eruption. Future volcanic hazard analysis will have to take into account this result. Our work also shows that multiple sampling of single lava flows is important to obtain a reliable mean direction. Sampling sites have to be carefully selected so that they represent un-tilted parts of the flows. We interpret this to be the case for the Ceboruco lava flow, while three of the six sites of the 1870 lava flow may have been partly or completely affected by movements after thermoremanent magnetization acquisition. Unfortunately, no better sites were found for this flow.

Description: Archaeomagnetic observations are key to recovering the behaviour of the geomagnetic field over the past few millennia. The corresponding data set presents a highly heterogeneous distribution in both space and time. Furthermore, the data are affected by substantial age and experimental uncertainties. In order to mitigate these detrimental properties, time-dependent global archaeomagnetic field models are usually constructed under spatial and temporal regularization constraints, with the use of bootstrap techniques to account for data uncertainties. The models so obtained are the product of an adjustable trade-off between goodness-of-fit and model complexity. The spatial complexity is penalized by means of a norm reflecting the minimization of Ohmic dissipation within the core. We propose in this study to resort to alternative spatial constraints relying on the statistics of a numerical dynamo simulation with Earth-like features. To that end, we introduce a dynamo norm in an ensemble least-squares iterative framework, the goal of which is to produce single-epoch models of the archaeomagnetic field. We first validate this approach using synthetic data. We next construct a redistributed archaeomagnetic data set between 1200 BC and 2000 AD by binning the data in windows of 40-yr width. Since the dynamo norm is not adjustable, we can legitimately calculate a resolution matrix to quantify the resolving power of the available archaeomagnetic data set. Gauss coefficients are resolved up to spherical harmonic degree 3 for the first thousand years of the interval, to degree 4 for the next thousand years and to degree 5 during the last millennium. These conclusions are based on the distribution and uncertainties that characterize the data set, and do not take into account the possible presence of outliers. Comparison between our model, called AmR, and previously published archaeomagnetic field models confirms the archaeomagnetic resolution analysis: it highlights the dichotomy between data-driven coefficients for which model predictions coincide (within their respective uncertainties), and prior-driven coefficients. This study opens the way to physics-based models of the archaeomagnetic field; future work will be devoted to integrating the framework here introduced into a time-dependent ensemble assimilation scheme.

Description: We carried out an archaeomagnetic directional study on 38 oriented samples (bricks and baked clays) collected from four archaeological locations at three provinces in China. The ages of our samples, spanning from ~3000 BCE to ~1300 CE, were constrained using a combination of archaeological context, radiocarbon dating and stratigraphic information. Rock magnetic results demonstrate that the main magnetic minerals of the studied samples are magnetite and/or hematite in single domain and superparamagnetic states. A total of 20 new reliable archaeodirectional data from 12 independent sites are obtained after thermal demagnetization experiments. These are the first set of archaeodirectional data in China produced since the 1990s. The published data are largely from the past 2 kyr and data from older time periods are rare. Our new data, especially those from period older than 3 ka, fill many gaps of the presently published dataset and will provide strong constraints on palaeosecular variation of the geomagnetic field in Eastern Asia and on the improvement of global models. Quite a few inflection points in the direction of the geomagnetic field are recorded in Eastern Asia over the past 10 kyr and some of them synchronize with the maximums or minimums of the palaeointensity. The palaeosecular variation rates are very low (based on present data distribution) before 2000 BCE and then start to increase and fluctuate afterward, which is generally consistent with the pattern of palaeointensity variations in this area.