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.     

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 .     

Teleseismic body wave radiation from a seismic source in a layered medium

Summary. A simple method is presented which combines the reciprocity theorem and the flat layer theory to yield teleseismic body wave radiation from seismic sources embedded in the Earth's crust. The source is represented by its equivalent body forces and can be quite general. The effect of Yucca Flat geology on explosion signals is studied in detail. In particular, the m _b— M _s relation is shown to be dependent on detonation medium and source depth. Application to shallow earthquake faults demonstrates the strong influence of free surface and layering on the shape of P- and S-wave teleseismic spectra.     

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.     

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:      

Odd zonal harmonics in the geopotential, from analysis of 28 satellite orbits

Summary. The geopotential is usually expressed as an infinite series of spherical harmonics, and the odd zonal harmonics are the terms independent of longitude and antisymmetric about the equator: they define the 'pear-shape' effect. The coefficients J ₃, J ₅, J ₇, … of these harmonics have been evaluated by analysing the variations in eccentricity of 28 satellite orbits from near-equatorial to polar. Most of the orbits from our previous determination in 1974 are used again, but three new orbits are added, including two at inclinations between 62° and 63°, which have been specially observed for more than five years by the Hewitt cameras. With the help of the new orbits and revised theory, we have obtained sets of J -coefficients with standard deviations about 40 per cent lower than before. A 9-coefficient set is chosen as representative, and is as follows (all × 10⁹): J ₃=– 2530 ± 4, J ₅=–245 ± 5, J ₇=–336 ± 6, J ₉=–90 ± 7, J ₁₁= 159 ± 9, J ₁₃=–158 ± 15, J ₁₅=– 20 ± 15, J ₁₇=– 236 ± 14, J ₁₉=– 27 ± 19. With this set of values, the pear-shape asymmetry of the geoid (north polar minus south polar radius) amounts to 45.1 m instead of the previous 44.7 m. The accuracy of the longitude-averaged geoid profile is estimated as 50 cm, except at latitudes above 86°. The geoid profile and predicted amplitude of the oscillation in eccentricity are compared with those from other sources.     

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.     

Anomalies of some tidal waves of UT1

Summary. The M _f and M _m waves of UT1 have been analysed from the BIH data during the period 1967.0 to 1984.0 in order to derive the Love number k .
These analyses performed during successive intervals of this period show some anomalies in the values of the Love number k derived from the M _f wave. Large variations with time of the amplitude and the phase appear for this wave while the values derived from the M _m wave present a good stability during the whole period.
The spectrum of the UT1-residuals (as obtained by removing the theoretical zonal tidal UT1 terms) shows the existence of a perturbation wave near the M _f period. This wave lies in the range 13.5–13.9 day according to the analysed interval; it could be related with the perturbation wave noticed by some authors near the M _f gravimetric wave.     

Finite-difference modelling of magnetotelluric fields in two-dimensional anisotropic media

An algorithm for the numerical modelling of magnetotelluric fields in 2-D generally anisotropic block structures is presented. Electrical properties of the individual homogeneous blocks are described by an arbitrary symmetric and positive-definite conductivity tensor. The problem leads to a coupled system of partial differential equations for the strike-parallel components of the electromagnetic field. E _x, and H _x These equations are numerically approximated by the finite-difference (FD) method, making use of the integro-interpolation approach. As the magnetic component H _x, is constant in the non-conductive air, only equations for the electric mode are approximated within the air layer. The system of linear difference equations, resulting from the FD approximation, can be arranged in such a way that its matrix is symmetric and band-limited, and can be solved, for not too large models, by Gaussian elimination. The algorithm is applied to model situations which demonstrate some non-trivial phenomena caused by electrical anisotropy. In particular, the effect of 2-D anisotropy on the relation between magnetotelluric impedances and induction arrows is studied in detail.     

A comment on depths of sources of anomalies in the Earth's gravity field

Summary. One method to determine the depths of sources of anomalies in the Earth's gravity field is to plot log [(2 n + 1) σ _n ] versus n where σ _n is the n th term in the amplitude spectrum of the Earth's gravitational potential. This procedure assumes that the amplitude spectrum of the anomalous density variations does not vary with n. Such an assumption may not apply to the Earth.     

Range of validity of seismic ray and beam methods in general inhomogeneous media – I. General theory

Summary. The limitations of asymptotic wave theory and its geometrical manifestations are newly formalized and scrutinized. Necessary and sufficient conditions for the existence of acoustic and seismic rays and beams in general inhomogeneous media are expressed in terms of new physical parameters: the threshold frequency ω₀ associated with the P/S decoupling condition, the cut-off frequency ω_c associated with the radiation-zone condition, the total curvature of the wavefront and the Fresnel-zone radius.
The analysis is facilitated with the introduction of a new ancillary functional – the hypereikonal which is capable of representing ordinary as well as evanescent waves. The hypereikonal is the natural extension of the eikonal theory.
With the aid of the above new parameters, simple conditions are obtained for the decoupled far field, the decoupled near field, two point dynamic ray tracing, paraxial wavefields and Gaussian beams.     

A comparative study of the ARM and Isr of some natural magnetites of MD and PSD grain size

Summary. ARM and I _sr stability with respect to alternating field and low-temperature treatment were measured for a set of artificial rock specimens. Each of the specimens contained a known amount of a natural magnetic mineral of a known grain size. The following minerals were used: a homogeneous natural magnetite, a natural titanomagnetite with submicroscopic Fe-rich inclusions, and a natural titanomagnetite and a natural titano-maghemite that both show exsolution phenomena. The grain size of the fractions varied from < 5 to 250 μm. It is found that the stability of ARM and I _rs of multi-domain grains with respect to alternating fields and low-temperature treatment is similar, whereas for the smallest grains studied (< 5 μm) ARM is more stable than I _sr. Multi-domain ARM data from this study are in agreement with current multi-domain ARM theory. Low-temperature behaviour of ARM and I _sr can be explained either in terms of a model in which stable moments in the bulk of the grains are screened by the surrounding matrix or in terms of a model in which regeneration of the bulk remanence plays a role and in which stable moments are situated on the surfaces of the grains.     

Carbon dioxide and methane dynamics in a sub-Arctic peatland in northern Finland

We studied carbon dynamics on various surface parts of a highly patterned fen, typical in northern Finland, to examine the importance of different microsites to the areal carbon fluxes. The studies were carried out in June-September 1995 on a mesotrophic flark fen (an aapa mire) in Kaamanen (69°08'N, 27° 17'E). Wet flarks, moist lawns and dry strings accounted for 60%, 10% and 30% of the surface area, respectively. A static chamber technique was applied to measure the CH₄ exchange, the instantaneous net ecosystem exchange (NEE, transparent chamber) and the ecosystem respiration (R_tot' opaque chamber) in several microsites. The static chamber results were compared with those obtained by the eddy covariance technique. The mean daytime areal net ecosystem CO₂ exchange rate measurement in conditions where photosynthesis was light saturated (PAR>400 μmol m_-2 s_-1) varied during the measurement period from −59 mg CO₂-C m⁻²h⁻¹ (release) to 250 (uptake). The mean CH₄ emission during the measuring period was 78 mg CH_4-Cm⁻² d⁻¹ on the flarks, 68 mg on the lawn and 6.0 mg on the strings. The strings without shrubs (mainly Betula nana ) were in general net sources of CO₂, even during the middle of the growing season, whereas the lawns, flarks and also strings growing B. nana showed a daytime net uptake of CO₂. Areally integrated chamber results showed lower CO₂ and higher CH₄ fluxes than predicted from the eddy covariance measurements.     

Estimation of the Bingham distribution function on nearly two-dimensional data sets

Summary. In cases where directional data, such as palaeomagnetic directions, lie nearly along a great circle, a good approximation to the maximum likelihood estimate of the intermediate concentration parameter k ₂ in the Bingham probability distribution is given by: 2( t ₂/ N ) – 1 = I ₁(1/2 k ₂)/ I ₀(1/2 k ₂), where t ₂ is the intermediate eigenvalue, N is the number of samples, and the I_i are the appropriate modified Bessel functions of the first kind. This estimate, the asymptotic limit as the smallest eigenvalue t ₁→ 0, corresponds to restricting all points to lie on a great circle. The limit is also useful as an endpoint for interpolation, especially since numerical calculation in this region is difficult.     

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:      

Finite element calculations of very high Rayleigh number thermal convection

Summary. A finite element method with uniform and variable resolution meshes is used to model very high Rayleigh number Ra thermal convection in a square box of infinite Prandtl number, Boussinesq fluid with constant viscosity and thermodynamic properties. Heating is either entirely from below or mostly from within and the boundaries are stress free. The variable mesh is coarse in the interior of the convection cell and it is fine in the very thin boundary layers and plumes surrounding the core. The highest resolution variable mesh has a dimensionless grid spacing of 0.027 in the core and 0.0017 in the boundary layers. The boundary layers contain about 10 mesh points even at the highest values of Ra considered and are thus highly resolved. The variable mesh approach is shown to yield reliable simulations of convection as long as the aspect ratio of the most elongated boundary layer elements is not too large; values of about 4 to 6 work well. This aspect ratio also measures the increase in resolution in the boundary layers as compared with the central core. Steady single-cell rolls are computed for bottom heating and Ra up to 5 × 10⁵ times the marginal instability value of the Rayleigh number Ra_cr. One and two-cell roll solutions are calculated for f = 1, 0.8 and 0.6, where f is the fraction of the heat escaping through the top of the box that is generated internally. The values of Ra_cr for f = 1, 0.8 and 0.6 are 1296, 1024 and 864, respectively. The largest of Ra/Ra_cr at which unicellular convection is stable (steady) are approximately 390, 610 and 970, for f = 1, 0.8 and 0.6.     

Tidal gravity in Britain: tidal loading and the spatial distribution of the marine tide

Summary. Tidal gravity measurements have been made at six sites in Britain with two nulled LaCoste and Romberg Earth tide gravitymeters. The M ₂ observations from these and two further sites are compared with calculations of the tidal loading from the seas around the British Isles and the major oceans. Models of the M ₂ marine tides are convolved with Green's functions for appropriate radially stratified Earth models. The differences between the M ₂ observations and the theoretical calculations are less than 0.6 μ gals and it is shown that these differences contain further information concerning the errors in the marine tide models. The M ₂ marine tides on the north-west European continental shelf are reasonably well known and this allows a useful test of the feasibility of using tidal gravity measurements for the inverse ocean tide problem in areas where the ocean tides are less well known. The differential gravity loading signal between pairs of gravity stations is shown to be important for considerations of the uniqueness and accuracy of the inverse problem. M ₂ tidal gravity loading maps for the British Isles and Europe have been produced which are of use in making corrections to various geodetic measurements.     

Some comments on the descriptions of the polarization states of waves

Summary. Procedures are developed for specifying the polarization characteristics of n -dimensional waves, and in particular three-dimensional waves of geophysical interest. We show that when a wave is in a pure state or is totally polarized, all the polarization information can be represented by a single vector u in an n -dimensional unitary space. Simple measures of the degree of polarization of the wave are constructed from the characteristic equation of the spectral matrix S . These measures are functions only of the scalar invariants of S and consequently S need not be diagonalized. If S represents a purely polarized wave, the unitary vector which contains the polarization information about the wave can be obtained directly from S using any 2 n – 1 equations of n ² possible equations. By multiplying by a phase-factor this unitary vector can be written in the form u = r ₁+ i r ₂ where r ₁ and r ₂ are orthogonal vectors in a real space. For an elliptically polarized wave, r ₁ and r ₂ locate the major and minor axes of the ellipse, and the ellipticity is given by the ratio of their magnitudes. The polarization parameters of ULF magnetic waves at the Earth's surface are computed from one set of five equations ( n = 3) and compared with parameters calculated using established techniques.