Relationship between median destructive field and remanent coercive forces for dispersed natural magnetite, titanomagnetite and hematite

Summary. Remanent acquisition curves, remanent hysteresis curves and alternating field demagnetization curves were determined for a number of artificial rock specimens containing well-defined grain-size fractions between 5 and 250 μm of natural magnetite, titanomagnetite and hematite. From these curves, the remanent acquisition coercive force H '_cr, the remanent coercive force H _cr and the median destructive field of IRM H _½I were determined. Theoretically these parameters should be the same for an assembly of non-interacting, homogeneously distributed, randomly oriented single-domain grains. For a given hematite specimen H '_cr, H _cr and H _½I have about the same value in spite of the strong grain-size dependence of these parameters. For each specimen of magnetite and titanomagnetite the value of H '_cr is larger than H _cr which again is larger than H _½I. However, the ratios H '_cr/ H _cr and H _½I/ H _cr appear to have a (different) constant value. An interesting relationship which appears to hold for dispersed magnetite, titanomagnetite or hematite grains between 5 and 250 μm, independently of grain-size, quantity and packing density of the magnetic material, is:      

Multiple scattering and energy transfer of seismic waves – separation of scattering effect from intrinsic attenuation – I. Theoretical modelling

Summary. In order to separate the scattering effect from the intrinsic attenuation, we need a multiple scattering model for seismic wave propagation in random heterogeneous media. In this paper, we apply radiative transfer theory to seismic wave propagation and formulate in the frequency domain the energy density distribution in space for a point source. We consider the cases of isotropic scattering and strong forward scattering. Some numerical examples are shown. It is seen that the energy density–distance curves have quite different shapes depending on the values of medium seismic albedo B ₀=η_s/(η_s+η_a) where η_s is the scattering coefficient and η_a is the absorption coefficient of the medium. For a high albedo ( B > 0.5) medium, the energy–distance curve is of arch shape and the position of the peak is a function of the extinction coefficient of the medium η_e=η_s+η_a. Therefore it is possible to separate the scattering effect and the absorption based on the measured energy density distribution curves.     

Magnetotelluric response of a uniformly stratified earth containing a magnetized layer

Summary The magnetotelluric (MT) response is studied of a uniformly stratified earth which contains a magnetized layer. The impedance as a function of the layer parameters (resistivity, ρ permeability, μ and thickness h ) is discussed. The MT response from a layer (μ, ρ, h ) is equivalent to that from a layer (μ/μ_r=μ₀, μ_rρ, μ_r h ) where μ_r is the relative permeability of the layer. Thus the effect of a magnetized layer is to make it apppear μ_r times more resistive and μ_r times thicker than an unmagnetized equivalent layer. Master curves of apparent resistivity and phase are computed for three-layer models with varying permeability associated with varying resistivity in each layer. An example of MT field data is presented in which the most reasonable interpretation is that a magnetized layer exists beneath the observatory site.     

The solution of the inverse gravity potential problem for cylindrical bodies

Summary. The inverse gravity potential problem consists in the determination of the form and the density of the body by its exterior gravity potential. We describe two similar classes of bodies for which this problem has a unique constructive solution.
(1) The first class contains the cylindrical bodies with finite length, arbitrary form of section and ρ( R , ø, z) =ρ₁( z )ρ₂( R , ø) density distribution, where z is the cylindrical coordinate; R , ø are the polar coordinates in a section plane. This class is important for prospecting geophysics in that it allows us to determine in a unique and constructive way, the function ρ₁( R , ø), the length, form and orientation of the cylinder if we know the function ρ₁( z ) and the exterior potential. The classical moment problem of functions is the basis for the solution of this problem.
(2) The analogous problem for the class of the spherical cylinders, or bodies bounded by arbitrary similar sections of two different concentric spheres and the radial lateral surface, appears when bodies of planetary size are studied. (An example of these bodies would be the Moon mascons.) The density distribution of these cylinders is ρ(τ, θ, ø) =ρ₁(τ)ρ₂(θ, ø) where τ, θ, ø are the spherical coordinates. The function ρ₁(θ, ø), length and form of spherical sections can be uniquely determined by exterior potential if we know the function ρ₁(τ). We propose a new constructive method for harmonic continuation of the gravity potential into the region containing the perturbing masses for the solution of the problem.     

The Influence of Surface Inhomogeneities On Deep Electromagnetic Soundings of the Earth

Summary. Asymptotic expressions for components of the electromagnetic field of a grounded electric dipole are considered for the model consisting of a thin surface-layer overlapping a stratified medium with a highly resistive screen on the roof. It is shown that the method of spatial derivatives makes it possible to obtain proper estimates of the impedance at distances of r ≥|λ₀| from the nearest edge of the surface anomaly (|λ₀| being the effective depth of the field penetration in the underlying section). the magnetotelluric methods allow one to obtain the true values of impedance, provided r ≥ max {|λ₀|, |/( S ⁻¹+ Z ⁻₀|^1/2} where S is the integrated conductivity of the surface layer, is the transverse resistance of the screen, and Z ⁻₀ is the Tikhonov—Cagniard impedance for the medium underlying the surface layer.     

A study of geomagnetic storms

Summary. An attempt is made to find interplanetary magnetic field and solar-wind parameters which control the development of geomagnetic storms. For this purpose, the interplanetary energy flux is estimated in terms of the Poynting flux ( E × B /4π), and its time variations are compared with the rate of energy dissipation in terms of the ring-current particle injection u _i( t ), Joule dissipation in the ionosphere u_j ( t ) and auroral particle injection u_p ( t ) for 15 major geomagnetic storms.
It is shown that the growth of geomagnetic storms, namely the time variations of the rate of the total energy dissipation, u ( t ) = u _i( t ) + u _j( t ) + u _p( t ), is closely related to the Poynting flux by the following relation:
where l ₀≅ 7 R _E and θ' is a measure of the angle between the interplanetary magnetic field vector and the magnetospheric field vector at the front of the magnetosphere in the equatorial plane. Further, it is shown that within a factor of 2 for each storm period.
A large increase of u ( t ) is associated with substorm activity. Thus, the energy flux ɛ( t ) entering the magnetosphere is dissipated through magneto-spheric substorm processes within the magnetosphere, and their accumulated effects can be understood as geomagnetic storm phenomena.     

Experimental study of the anisotropy of longitudinal and transverse waves from local earthquake records

Summary. The paper gives the results of a study of the anisotropy of seismic wave velocities within the Ashkhabad test field in Central Asia. The anisotropy was studied by analysing variations in the values of apparent velocities of first arrivals for epicentral distances ranging from 30 to 130 km and by analysing the delays (Δ t_s1-_s2 ) between the arrival times of shear waves with different polarizations.
The velocities of P -waves vary with azimuth from 5.3 to 6.27 km s^-1 and the velocities of S -waves vary from 3.15 to 3.5 km s^-1.
The delay times Δ t_S1 - _S2 depend on the direction of the propagation. The character of the variation of the propagation velocity of the longitudinal wave, the presence of two differently polarized shear waves S 1 and S 2 propagating at different velocities, and the character of the distribution of Δ t_S1 - _S2 on the stereogram suggest that the symmetry of the anisotropic medium is close to hexagonal with a nearly horizontal symmetry axis coinciding with the direction of maximal velocity. The azimuth of the symmetry axis of the medium is 140° and coincides with the direction of geological faults.     

On the consistency of earthquake moment release and space geodetic strain rates: Europe

In this paper, approximately 100 VLBI/SLR/GPS velocities map European strain rates from <0.09 × 10⁻⁸ to >9.0 × 10⁻⁸ yr⁻¹ with regional uncertainties of 20 to 40 per cent. Kostrov's formula translates these strain-rate values into regional geodetic moment rates M¯˙ _geodetic . Two other moment rates, M¯˙ _seismic , extracted from a 100-year historical catalogue and M¯˙ _plate , taken from plate-tectonic models, contrast the geodetic rates. In Mediterranean Europe, the ratios of M¯˙ _seismic to M¯˙ _geodetic are between 0.50 and 0.71. In Turkey the ratio falls to 0.22. Although aseismic deformation may contribute to the earthquake deficit ( M¯˙ _seismic values less than M¯˙ _geodetic ), the evidence is not compelling because the magnitudes of the observed shortfalls coincide with the random variations expected in a 100-year catalogue. If the lack of aseismic deformation inferred from the 100-year catalogue holds true for longer periods, then much of Europe's strain budget would have to be accommodated by more frequent or larger earthquakes than have been experienced this century to raise the ratios of M¯˙ _seismic to M¯˙ _geodetic to unity. Improved geological fault data bases, longer historical earthquake catalogues, and densification of the continent's space geodetic network will clarify the roles of aseismic deformation versus statistical quiescence.     

Nature and origin of magnetic minerals within the Middle Jurassic shallow-water carbonate rocks of the Paris Basin, France: implications for magnetostratigraphic dating

The Middle and Upper Jurassic Bathonian-Oxfordian shallow-water carbonate rocks from the Paris Basin, France, consist mainly of oolitic and bioclastic limestones that are hydrocarbon reservoirs in the subsurface. Despite a preliminary positive study, these deposits have been considered to be largely remagnetized (Rochette, private communication), and hence not amenable to palaeomagnetic dating. To establish their magnetic mineralogy and test this remagnetization hypothesis, we have used an integrated investigation combining petrographic, geochemical, rock-magnetic and palaeomagnetic measurements on samples extracted from five cores from the Paris Basin and from outcrops in Burgundy. Magnetic minerals in the Bathonian-Oxfordian carbonates include: (1) primary biogenic single-domain magnetite and detrital multi-domain Ti-magnetite and their oxidized form, maghemite; (2) authigenic spheres of magnetite probably related to hydrocarbons; and (3) goethite, either restricted to ferruginous ooid layers or resulting from surficial alteration, notably replacement of pyrite framboids. Rock-magnetic experiments carried out on 68 samples reveal H _cr/ H _c and M _rs/ M _s ratios ranging from 1.88 to 5.58 and 0.017 to 0.314, respectively. These values are clearly distinct from diagnostic values for a chemical remagnetization. Pyrrhotite was not identified within these sediments. Moreover, the average H _cr/ H _c ratio of 3.14 is significantly different from the value of 1.333 for natural pyrrhotite (Dekkers 1988). These results have a direct implication for the preservation of the primary magnetization; consequently, these deposits are selectively amenable for magnetostratigraphic dating and possible regional correlations.     

Rayleigh-wave amplitudes from earthquakes in the range 0–150°

Summary . Vertical component Rayleigh-wave amplitudes from 1461 shallow earthquakes recorded in the distance range 0–150° are analysed to separate the effects of earthquake size, epicentral distance (Δ) and recording station.
The estimated decay of amplitude with distance has the form of a theoretical curve for the decay of Rayleigh waves with distance if the assumption is made that the decay due to dispersion for the data analysed is that of an Airy phase. Writing the decay due to anelastic attenuation as exp (- k Δ), k is estimated to be 0.676/rad over the whole range of distance. If the distance effects are represented by a straight line of the form h log Δ+ constant, h is estimated to be 1.15. The calibration function for computing M _s derived from the estimated distance effects is very similar to that of Marshall & Basham.
Station effects on Rayleigh-wave amplitudes though statistically significant are small, and can probably be ignored in the computation of M _s.
Comparing the estimated surface-wave magnitudes (earthquake size) obtained in this study with the long and short period body-wave magnitudes ( m ^LP_b and m ^SP_b respectively) obtained by Booth, Marshall & Young for the same earthquake shows that m ^LP_b is about equal to M _s over the magnitude range of interest (˜4.0–7.0). The m ^LP_b and M_s relationship shows that the greater the long-period energy radiated by an earthquake the smaller proportionately is the short-period energy.     

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.     

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.     

Two-dimensional non-linear inversion of VLF-R data using simulated annealing

The VLF-R (very low frequency-resistivity) data, i.e. the apparent resistivity ( ρ _a ) and phase ( φ ) data, were inverted individually and jointly using the VFSA (very fast simulated annealing) global inversion approach. Global inversion results for synthetic data without and with various amounts of random and normally distributed Gaussian noise reveal that the inversion of neither the ρ _a nor φ data alone yields the true parameters of the structures. However, the joint inversion of the ρ _a and φ data yields very good estimates of the model parameters. Five models, representing typical subsurface structures in the shield areas, are studied here. Various models achieved after 10 VFSA runs were used to compute the mean model and the corresponding covariance and correlation matrices, which were used to estimate the uncertainties in the mean model parameters and correlations between the model parameters. We observe that these correlations follow the physics associated with the problem. VLF-R field data due to a nearly vertical contact structure and a very thick dyke-like structure were also inverted to demonstrate the efficacy of the approach in the delineation of the parameters of 2-D structures.     

Detailed gravity-tide spectrum between one and four cycles per day

Summary. The results of previous work by the authors is used to remove most of the effects of ocean and atmospheric loading from an 18-month Earth gravity-tide record. The remaining signal is examined for additional influence of ocean and atmosphere and for evidence of the frequency-dependence of the response of the solid earth. Variations in time of the measured tides are shown to result from the atmospheric tide at S ₂ and appear to result from variations in ocean tides at other frequencies. The frequency-dependence of the solid earth response near 1 cycle per siderial day is found to be consistent with the nearly diurnal free wobble. However, the influence of the ocean on the small but crucial Ψ₁ tide is uncertain. Anomalous responses are observed at several other frequencies but except for the case of ρ₁ it is argued that anomalous ocean tides are plausible and could therefore explain the observations.     

Focal mechanisms of recent earthquakes in the Southern Korean Peninsula

We evaluate the stress field in and around the southern Korean Peninsula with focal mechanism solutions, using the data collected from 71 earthquakes ( M_L = 1.9–5.2) between 1999 and 2004. For this, the hypocentres were relocated and well-constrained fault plane solutions were obtained from the data set of 1270 clear P -wave polarities and 46 SH / P amplitude ratios. The focal mechanism solutions indicate that the prevailing faulting types in South Korea are strike-slip-dominant-oblique-slip faultings with minor reverse-slip component. The maximum principal stresses (σ₁) estimated from fault-slip inversion analysis of the focal mechanism solutions show a similar orientation with E–W trend (269°–275°) and low-angle plunge (10°–25°) for all tectonic provinces in South Korea, consistent with the E–W trending maximum horizontal stress (σ_Hmax) of the Amurian microplate reported from in situ stress measurements and earthquake focal mechanisms. The directions of the intermediate (σ₂) and minimum (σ₃) principal stresses of the Gyeongsang Basin are, however, about 90 deg off from those of the other tectonic provinces on a common σ₂–σ₃ plane, suggesting a permutation of σ₂ and σ₃. Our results incorporated with those from the kinematic studies of the Quaternary faults imply that NNW- to NE-striking faults (dextral strike-slip or oblique-slip with a reverse-slip component) are highly likely to generate earthquakes in South Korea.     

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.     

On the linear relation between mb and Ms for discrimination between explosions and earthquakes

Summary. The statistical capability of the m _b: M _s discriminant for the discrimination of earthquake and explosion populations is examined by application of discriminant functions to a group of 83 explosions and 72 earthquakes in Eurasia. Equations are derived for the probability that an event is an earthquake or an explosion. The positive sign of DIS in the decision index equation, DIS _i = 34.3383 – 11.9569 mb _t + 7.1161 M _si , indicates that the event i is an earthquake. Its negative sign indicates that event i is an explosion. The probability of correct classification for an event, P _i , is related to its DIS _i value, by P _i = [1-exp (DIS _i )]⁻¹, where a large, positive DIS indicates a high probability that an event is an earthquake and a large, negative DIS indicates a high probability that an event is an explosion. The discrimination line M _s = 1.680 m _b– 4.825, or m _b= 0.595 M _s+ 2.872 very successfully separates the explosion population from the earthquake population. The points on this line have an equal chance of being an earthquake or an explosion; moreover, for any event, the distance parallel to the M _s-axis from the point representing that event in the m _b: M _s plane to this line is a measure of the probability for the correct classification of that event.     

A New Magnitude-Frequency Relation for Earthquakes and a Classification of Relation Types

Summary. In this paper some aspects of the deviations of frequency distributions of earthquakes from a linear trend are considered. For recurrence curves with inflexions, the formula
is proposed as a magnitude-frequency relation, with critical magnitude M _cr. Comparison with actual data shows better agreement than other formulae. The slope of the frequency graph is considered in the most general case to be a function of M and some parameters p ₁, p ₂, …. A classification of the types of magnitude-frequency relation is given and all known types of relations N : M and N: K ( K = log E , E = seismic energy) are summarized.     

The traveltime perturbations for seismic body waves in factorized anisotropic inhomogeneous media

The traveltime perturbation equations for the quasi-compressional and the two quasi-shear waves propagating in a factorized anisotropic inhomogeneous (FAI) media are derived. The concept of FAI media simplifies considerably these equations. In the FAI medium, the density normalized elastic parameters a _ijkl ( X _i ) can be described by the relation a _ijkl ( X _i) = f ²( x _i ) A _ijkl, where A _ijkl are constants, independent of coordinates x _i and f ²( x _i) is a continuous smooth function of x _i . The types of anisotropy ( A _ijkl ) and inhomogeneity [ f ( x _i)] are not restricted. The traveltime perturbations of individual seismic body waves ( q ^P , qS 1 and qS 2) propagating in the FAI medium depend, of course, both on the structural pertubations [δ f ²( x _i)] and on the anisotropy perturbations (δ A _ijkl ), but both these effects are fully separated. The perturbation equations for the time delay between the two qS -waves propagating in the FAI medium are simplified even more. If the unperturbed (background) medium is isotropic, the perturbation of the time delay does not depend on the structural perturbations (δ f ²( x _i) at all. This striking result, valid of course only in the framework of first-order perturbation theory, will simplify considerably the interpretation of the time delay between the two split qS -waves in inhomogeneous anisotropic media. Numerical examples are presented.     

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.