Magnetic fabric and its significance in the Florianópolis dyke swarm, southern Brazil

Magnetic fabric was determined by applying the anisotropy from the low-field magnetic susceptibility (AMS) technique in 62 mafic dykes from the Mesozoic Florianópolis (Santa Catarina Island) dyke swarm, southern Brazil. These dykes cut the crystalline basement rocks, which are mainly Proterozoic. They are vertical or subvertical in dip and trend mainly NE, although NW-trending dykes are also found. Dykes are tholeiitic in composition and are geochemically similar to those from the Ponta Grossa swarm. Thicknesses vary from 0.3 to 60 m. Polished sections show that titanomagnetites carry the AMS in these dykes. Hysteresis parameters show that the magnetic minerals fall in the PSD range. Two types of magnetic fabric are recognized. Type I is characterized by K ₁- K ₂ parallel to the dyke wall, representing magma flow within the dykes; type II, with K ₁- K ₃ parallel to the dyke wall, was found in four dykes. Type I is found in 94 per cent of the dykes, and approximately 20 per cent of these have K ₁ inclinations of less than 30°, suggesting horizontal or subhorizontal flow. About 80 per cent have K ₁ inclinations of greater than 30°, due to inclined to vertical flow. The comparison of AMS studies from both the Florianópolis and the Ponta Grossa dykes suggests a source position closer to Santa Catarina Island than the Ponta Grossa arch.     

Magnetic susceptibility anisotropy measurements of coal from the South Wales coal field: a coalification process indicator

Summary. Anisotropy of magnetic susceptibility measurements have been performed on 154 specimens taken from 25 different seams in the South Wales Coal Field, up to 24 samples per seam. The magnetic fabric determined is represented, in general, by oblate spheroids with minimum ( K _c) AMS axes usually clustered at a near-vertical, or occasionally, near-horizontal pole which is normal to the magnetic foliation plane defined by great circle distributions of maximum ( K _a) and intermediate ( K _b) axes. The results suggest that the coalification process and variations in rank within the coal field are due to the combined effects of vertical compression and compaction resulting from variations in overburden pressure, and to horizontal compression, due to variations in tectonic thrusting during the Variscan Orogeny. Analysis of the AMS of closely spaced samples across a single anthracite seam shows no systematic variability in magnetic properties within the seam.     

Multimode surface waveform tomography of the Pacific Ocean: a closer look at the lithospheric cooling signature

We present a regional surface waveform tomography of the Pacific upper mantle, obtained using an automated multimode surface waveform inversion technique on fundamental and higher mode Rayleigh waves, to constrain the   V_SV   structure down to ∼400 km depth. We have improved on previous implementations of this technique by robustly accounting for the effects of uncertainties in earthquake source parameters in the tomographic inversion. We have furthermore improved path coverage in the South Pacific region by including Rayleigh wave observations from the French Polynesian Pacific Lithosphere and Upper Mantle Experiment deployment. This improvement has led to imaging of vertical low-velocity structures associated with hotspots within the South Pacific Super-Swell region. We have produced an age-dependent average cross-section for the Pacific Ocean lithosphere and found that the increase in   V_SV   with age is broadly compatible with a half-space cooling model of oceanic lithosphere formation. We cannot confirm evidence for a Pacific-wide reheating event. Our synthetic tests show that detailed interpretation of average   V_SV   trends across the Pacific Ocean may be misleading unless lateral resolution and amplitude recovery are uniform across the region, a condition that is difficult to achieve in such a large oceanic basin with current seismic stations.     

Further evidence for oceanic excitation of polar motion

While the role of the atmosphere in driving variations in polar motion is well established, the importance of the oceans has been recognized only recently. Further evidence for the role of the oceans in the excitation of polar motion is presented. To estimate the equatorial excitation functions, χ ₁ and χ ₂ , for the ocean, we use velocity and mass fields from a constant-density ocean model, driven by observed surface wind stresses and atmospheric pressure, for the period 1993–1995; comparison with similar functions derived from a more complex density-stratified ocean model indicates the effectiveness of the simple constant-density modelling approach. Corresponding atmospheric excitation functions are computed from NCEP/NCAR re-analyses. Results indicate significant improvements in the agreement with the observed polar motion excitation when the simulated oceanic effects are added to atmospheric excitation. Correlations between the polar motion and the geophysical signals at periods of 15–150 days increase from 0.53 to 0.80 and from 0.75 to 0.88 for χ ₁ and χ ₂ , respectively. The oceanic signals are particularly important for seasonal variations in χ ₁ (correlation increases from 0.28 to 0.85 when oceanic excitation is included). A positive impact of the oceans on more rapid polar motion is also observed, up to periods as short as 5 days. The sensitivity of the results to different forcing fields and different amounts of friction in the oceans is also discussed.     

Palaeomagnetism of alkaline lava flows from El Salto-Almafuerte, Córdoba Province, Argentina

Summary. Stable components of magnetization have been isolated in 15 lava flows (mean K-Ar age 123 ± 4 Myr) from the alkaline sequence outcropping at El Salto-Almafuerte, Province of Cordoba, Argentina. Magnetic and geologic stratigraphy, as well as K-Ar ages indicate that this sequence was probably extruded in the Lower Cretaceous during the first volcanic cycle of the Sierra de los Cóndores Group (Vulcanitas Cerro Colorado Formation).
The palaeomagnetic pole-position for El Salto-Almafuerte lava flows, computed from the mean of 15 virtual geomagnetic poles and denoted SAK₇, is: 25° E, 72° S ( k = 35, α₉₅= 6.5°); it is fairly close to other Lower Cretaceous palaeomagnetic poles for South America. The elongated distribution of Cretaceous palaeomagnetic poles suggest recurrent drift for South America in early Cretaceous time.
The palaeomagnetic and radiometric data for the igneous rocks from El Salto-Almafuerte support the magnetic reversal time-scale for the early Cretaceous suggested by oceanic magnetic lineations.     

Observation of Tidal Tilt On Kilauea Volcano, Hawaii

Summary. We have analysed the east-west tilt components, O₁, K₁, N₂, M₂ and S₂ from a continuously recording tiltmeter located in Uwekahuna Vault on Kilauea Volcano, Hawaii, for the period 1971—79. Detailed analysis of the M₂ component gives values of 30.9 ± 2.0 (95 per cent) nrad and 116.0 ± 2.0° for the amplitude and phase, respectively, compared to values of 48.5 nrad and 139.4° for the equilibrium tide. the total theoretical tide, found by summing the equilibrium and load tides, amounts to 37.2 nrad at a phase of 121.7°. the 20 per cent discrepancy with that observed may be due to an inaccurate cotical chart, cavity effects in the vault, strain—tilt coupling or an inappropriate solid earth model. In the vicinity of Hawaii (≤ 3°) two independent cotidal charts give almost identical results for the near field ocean load. At greater distances, we use the Schwiderski (1978) cotidal chart. We estimate that local cavity and strain—tilt coupling effects are less than 12 per cent owing to the agreement between geodetically determined and instrumental tilt but we can not rule out regional effects. Assuming these are small and the cotical charts correct, we find that the M₂ results are brought into satisfactory agreement if, instead of using an average oceanic earth model in the (< 75 km) vicinity of Hawaii, we use an earth model with nearly one-half the oceanic rigidity. Such a low upper mantle and crustal rigidity is consistent with Kilauea's position above the thermal upwelling associated with the Hawaiian hotspot.     

The Canary Islands swell: a coherence analysis of bathymetry and gravity

The Canary Archipelago is an intraplate volcanic chain, located near the West African continental margin, emplaced on old oceanic lithosphere of Jurassic age, with an extended volcanic activity since Middle Miocene. The adjacent seafloor does not show the broad oceanic swell usually observed in hotspot-generated oceanic islands. However, the observation of a noticeable depth anomaly in the basement west of the Canaries might indicate that the swell is masked by a thick sedimentary cover and the influence of the Canarian volcanism. We use a spectral approach, based on coherence analysis, to determine the swell and its compensation mechanism. The coherence between gravity and topography indicates that the swell is caused by a subsurface load correlated with the surface volcanic load. The residual gravity/geoid anomaly indicates that the subsurface load extends 600 km SSW and 800 km N and NNE of the islands. We used computed depth anomalies from available deep seismic profiles to constrain the extent and amplitude of the basement uplift caused by a relatively low-density anomaly within the lithospheric mantle, and coherence analysis to constrain the elastic thickness of the lithosphere ( T_e ) and the compensation depth of the swell. Depth anomalies and coherence are well simulated with T_e =28–36 km, compensation depth of 40–65 km, and a negative density contrast within the lithosphere of ∼33 kg m⁻³. The density contrast corresponds to a temperature increment of ∼325°C, which we interpret to be partially maintained by a low-viscosity convective layer in the lowermost lithosphere, and which probably involves the shallower parts of the asthenosphere. This interpretation does not require a significant rejuvenation of the mechanical properties of the lithosphere.     

Coercive force of single crystals of magnetite at low temperatures

The temperature dependence of coercive force H _c was studied on well-characterized and stoichiometric millimetre-sized single crystals of magnetite at a series of 16 temperatures from 300 to 10 K using a SQUID magnetometer. H _c decreases gradually with cooling to the isotropic temperature, T _i = 130 K, where the first magnetocrystalline anisotropy constant K ₁ becomes zero. H _c exhibits a sharp increase at the Verwey transition, T _v = 120 K, where the structure changes from cubic to monoclinic. In crossing the Verwey transition, H _c increases by more than two orders of magnitude, from 20 μT to 2.4 mT, and the shape of the hysteresis loops becomes wasp-waisted.
Observed coercivity between 300 K and 170 K varies with temperature as λ _s / M _s , where λ _s is the magnetostriction constant and M _s is the saturation magnetization, indicating that the coercivity in MD magnetite is controlled mainly by internal stress associated with dislocations or other crystal defects. It seems likely that the stable single-domain-like magnetic memory observed in large MD magnetite crystals is due to magnetoelastically pinned domain walls. The discontinuous change in H _c at the Verwey transition is controlled by abrupt changes in magnetocrystalline and magnetostriction constants due to crystal deformation from cubic to monoclinic structure.     

Magnetic properties of metal-substituted haematite

Mineral and isothermal magnetic properties of Al-, Mn- and Ni-substituted haematites were characterized and their relationships evaluated in order to interpret better the results of magnetic analyses of soils and recent sediments. Aluminium, manganese and nickel haematites generally behaved as single-domain (SD) particles. The influence of incorporated Al on the magnetic behaviour of haematite was consistent with Al acting as a paramagnetic dilutent. Mass magnetic susceptibility ( χ ) and SIRM₈₀₀ decreased as the level of Al substitution increased. Incorporation of Mn and Ni increased χ , which could be associated with enhancement of the spin canting effect of haematite. The stability of SIRM₈₀₀ to demagnetization for Al-haematite appears to be related to a defect mechanism associated with the development of smaller crystallites arising from Al substitution. Magnetic domain rotation or flipping was probably inhibited, being blocked by structural defects during magnetization and demagnetization, and resulted in a low but stable partial SIRM (SIRM₈₀₀ ). %IRM/SIRM₈₀₀ demagnetization curves and estimated ( B _o )_CR values of ≤100  mT for Mn-haematite indicate pseudo-single-domain/multidomain-like behaviour despite Mn-haematite having particle and crystallite dimensions similar to Ni-haematite, which did not show this behaviour. Data indicate that parameters involving unsaturated, partial SIRM should be used with caution in magnetic studies of soils and sediments.     

TRM deviations in anisotropic assemblages of multidomain magnetite

Anisotropic assemblages of multidomain magnetite particles develop an anisotropy of magnetic susceptibility (AMS), which in turn induces deviations of thermo-remanent magnetization (TRM) from the field direction. From the theories of multidomain TRM acquisition, it is shown that the TRM anisotropy tensor has its eigenvalue ratios ( T _i) related to the principal weak-field susceptibility ratios ( P _i) by the order of magnitude T _i≃ P ²_i. This relation has been experimentally verified on two sets of highly anisotropic rock samples. The exponent has been determined to be 1.94 in the samples from a Peruvian gabbro, and 1.81 in those from the granite of Flamanville (NW France). Accounting for experimental difficulties in determining the TRM anisotropy tensors, these exponents are judged to agree well with the expected one. It is therefore stressed that AMS measurements provide a good means of evaluating the magnetic field direction from deviated TRM directions, providing magnetic carriers are mainly multidomain magnetites.     

Electromagnetic induction fields in the deep ocean north-east of Hawaii: implications for mantle conductivity and source fields

Summary. We have analysed a thirty-six day recording of the natural electric and magnetic field variations obtained on the deep ocean floor north-east of Hawaii. The electromagnetic fields are dominated by tides which have an appreciable oceanic component, especially in the east electric and north magnetic components. The techniques of data analysis included singular value decomposition (SVD) to remove uncorrelated noise. There are three degrees of freedom in the data set for periods longer than five hours, indicating a correlation of the vertical magnetic field and the horizontal components, suggesting source field inhomogeneity. Tensor response functions were calculated using spectral band averaging with both SVD and least squares techniques and rotated to the principal direction. One diagonal component, determined mainly by the north electric and east magnetic fields, is not interpretable as a one-dimensional induction phenomenon. The other diagonal term of the response function indicates a rapid rise in conductivity to 0.05 mho m⁻¹ near 160 km. No decrease in conductivity below this depth is resolvable. Polarization analysis of the magnetic field indicates moving source fields with a wavelength near 5000 km. Model studies suggest that the two dimensionality in the response function may be caused by motion in the ionospheric current system.     

A three-dimensional radially anisotropic model of shear velocity in the whole mantle

We present a 3-D radially anisotropic S velocity model of the whole mantle (SAW642AN), obtained using a large three component surface and body waveform data set and an iterative inversion for structure and source parameters based on Non-linear Asymptotic Coupling Theory (NACT). The model is parametrized in level 4 spherical splines, which have a spacing of ∼ 8°. The model shows a link between mantle flow and anisotropy in a variety of depth ranges. In the uppermost mantle, we confirm observations of regions with   V_SH > V_SV   starting at ∼80 km under oceanic regions and ∼200 km under stable continental lithosphere, suggesting horizontal flow beneath the lithosphere. We also observe a   V_SV > V_SH   signature at ∼150–300 km depth beneath major ridge systems with amplitude correlated with spreading rate for fast-spreading segments. In the transition zone (400–700 km depth), regions of subducted slab material are associated with   V_SV > V_SH   , while the ridge signal decreases. While the mid-mantle has lower amplitude anisotropy (<1 per cent), we also confirm the observation of radially symmetric   V_SH > V_SV   in the lowermost 300 km, which appears to be a robust conclusion, despite an error in our previous paper which has been corrected here. The 3-D deviations from this signature are associated with the large-scale low-velocity superplumes under the central Pacific and Africa, suggesting that   V_SH > V_SV   is generated in the predominant horizontal flow of a mechanical boundary layer, with a change in signature related to transition to upwelling at the superplumes.     

Rapid computation of magnetic anomalies with demagnetization included, for arbitrarily shaped magnetic bodies

Summary. The potential function ø for a magnetic body of susceptibility μ in a medium of susceptibility μ* satisfies the integral equation
Here Φ* is the potential function for the region without the heterogeneity and R is the distance from the point of observation to the point on the surface, s , of the body. δΦ /δn is the normal derivative, in the direction of the outward normal. The equation allows for the effects of demagnetization. For numerical purposes the surfaces can be divided into N facets over which δΦ/δ n is a constant. The unknown quantities δΦ/δn_j can be found from the system of equations defined by:
The prime on the summation sign denotes that the summation does not include the i th element. The magnetic field in the direction of the unit vector P( P ₁, P ₂, P₃ ) is given by:      

Palaeoclimatic record from a loess–soil profile in northeastern Bulgaria—I. Rock magnetic properties

A detailed record of mineral magnetic properties of a loess–palaeosol profile comprising seven loess horizons, six interbedded palaeosols and recent soil at the top in NE Bulgaria is analysed. A strong contrast between the soil and loess susceptibilities as well as other concentration-dependent hysteresis parameters is present, similar to the well-documented magnetic characteristics of the Chinese loess ( Hus & Han 1992 ; Maher & Thompson 1992 ; Heller & Evans 1995 ; Hunt et al. 1995 ). The magnetic enhancement of the palaeosol units is caused by very fine-grained pedogenic magnetite with superparamagnetic behaviour. Thermomagnetic analyses on bulk material suggest magnetite and maghemite as the main ferrimagnetic carriers in both soil and loess horizons. Their relative proportions are shown to reflect different palaeoclimatic conditions. Chernozem soils, which include recent soil S₀ and first and second palaeosols S₁ and S₂ developed under steppe vegetation, show a high degree of low-temperature oxidation of the pedogenic magnetite to maghemite. This material is characterized by coercive force H _c showing even higher values than those of the parent loess material. The older palaeosols (S₄ to S₆ ) were formed during more humid climatic conditions and therefore probably developed as forest types. Rock magnetic data suggest the existence here of only partly oxidized magnetite grains. The behaviour of the thermomagnetic curves, characterized by a kink at 200 °C, may be due to either a release of internal stress (built up as a result of partial low-temperature oxidation) or interactions between two phases.     

Three-dimensional Rayleigh hysteresis of oriented core samples from the German Continental Deep Drilling Program: susceptibility tensor, Rayleigh tensor, three-dimensional Rayleigh law

Rayleigh hysteresis, as defined by the well-known Rayleigh relations, has been observed not only when magnetization of pyrrhotite-bearing KTB-samples is measured in parallel to a weak dc magnetic field, but also in experiments where field and measuring directions have been adjusted strictly perpendicularly to each other. Nine-tupels of independent Rayleigh hysteresis loops could thus be compiled. Their characteristic coefficients X ^ijk of initial susceptibility together with the Rayleigh loss coefficients α^jk have been proved to determine completely the samples' weak-field magnetic anisotropy. Interpreting the coefficient matrices ( X ^ijk) and (α^jk) as the tensor of initial susceptibility and the Rayleigh tensor, respectively, generalization of the isotropic Rayleigh relations in terms of corresponding tensor relationships has been suggested for the anisotropic case. Application to the KTB samples showed 3-D Rayleigh hysteresis measurements to be an excellent tool for rock magnetic analysis in terms of ore content and crystalline texture. In particular, a magnetocrystalline double texture of the basal planes of pyrrhotite precipitates and their [1120] directions of easy magnetization have been clearly detected. Surprisingly, the welt-known theorem α= const. X ²_I, formulated by Néel (1942) for the isotropic case, has been found to hold true even in tensor generalization (α_jk) = const ( X ²_jk). To reach sufficient resolution for the measurements performed, a sensitive vibrating coil magnetometer (VCM) has been developed.     

Significance of the sign changing of the imaginary arrows in geomagnetic induction investigation

Summary. In this investigation, we carry out a two-dimensional study of the dependence of the imaginary Parkinson arrows on the frequency of the inducing geomagnetic field. Our results demonstrate that the imaginary arrows reverse direction as the inducing period varies. Therefore, we consider that there is no way to fix a consistent sign convention for the imaginary arrows even when the time factor is taken into account. We find that in the twodimensional case there exists a characteristic period T _c at which the phase difference between the vertical and horizontal magnetic components is zero. It is anticipated that T _c is related to the parameters of the conductivity anomaly and the status of the half-space host.     

Scattering of waves by a quasi-thin body of arbitrary shape

Summary. The problem of scattering of a wave with a front of arbitrary shape by a curved quasi-thin body bounded by contour L of arbitrary shape is presented. The problem is solved using the Kirchhoff approximation (KA). The reflecting (scattering) properties of the body are given in terms of the reflection (transmission) coefficients. The routine technique of computation in the KA is complemented by significant improvements. The procedure of smooth continuation of the body surface S beyond its boundary contour L is suggested. The invariant presentation of the expressions for the scattered field u _s obtained in the paper is extremely important.
The leading term in the formula for the field u _s is obtained as a product u _s= u _ray W , where u _ray is the reflected wavefield computed, according to the formulae of ray theory, as if the body surface is boundless; W is the so-called weakening function, which takes account of the modification of the ray field caused by diffraction at the edges (the contour L ).
The arguments of the weakening function are dimensionless quantities depending on the relation between the body size and the radii (half-axes) r _f^{( i )} ( i = 1, 2) of the Fresnel zones on the surface S. The values of r _f^{( i )} are obtained in the paper. The method of computation of the function W turns out to depend essentially on the reciprocal position of the Fresnel zone and body contour. When the contour L is of a complex form, the surface S is divided into simple surfaces S_i , each of which is limited by a contour L , which permits the carrying out of a single valued (in the sense of the leading part of the field) smooth continuation of S_i .     

Three-dimensional VP and VP /VS models of the upper crust in the Friuli area (northeastern Italy)

3-D images of P velocity and P - to S -velocity ratio have been produced for the upper crust of the Friuli area (northeastern Italy) using local earthquake tomography. The data consist of 2565 P and 930 S arrival times of high quality. The best-fitting V _P and V _P / V _S 1-D models were computed before the 3-D inversion. V _P was measured on two rock samples representative of the investigated upper layers of the Friuli crust. The tomographic V _P model was used for modelling the gravity anomalies, by converting the velocity values into densities along three vertical cross-sections. The computed gravity anomalies were optimized with respect to the observed gravity anomalies. The crust investigated is characterized by sharp lateral and deep V _P and V _P / V _S anomalies that are associated with the complex geological structure. High V _P / V _S values are associated with highly fractured zones related to the main faulting pattern. The relocated seismicity is generally associated with sharp variations in the V _P / V _S anomalies. The V _P images show a high-velocity body below 6 km depth in the central part of the Friuli area, marked also by strong V _P / V _S heterogeneities, and this is interpreted as a tectonic wedge. Comparison with the distribution of earthquakes supports the hypothesis that the tectonic wedge controls most of the seismicity and can be considered to be the main seismogenic zone in the Friuli area.     

A diurnal resonance in the ocean tide and in the Earth's load response due to the resonant free 'core nutation'

Summary. The luni-solar forced nutations and body tide are believed to be resonant at frequencies near (1 + 1/460) cycle sidereal day⁻¹ as seen from the rotating Earth. This resonance is due to the Earth's rotating, elliptical fluid core. We show here that tides in the open ocean and the Earth's response to those tides must also be resonant at (1 + 1/460) cycle day⁻¹. We examine these resonant oceanic effects on the Earth's nutational motion and on the body tide. Effects on the forced nutations might be as large as 0.002 arcsec at 18.6 yr. The effects on the observed resonance in the body tide are more important. For tidal gravity, for example, the difference between K ₁ and 0 ¹ which is usually used to determine the resonance, can be perturbed by 30 per cent or more due to the oceanic resonance effects.