An algebraic formulation of geophysical inverse problems

Many geophysical inverse problems derive from governing partial differential equations with unknown coefficients. Alternatively, inverse problems often arise from integral equations associated with a Green's function solution to a governing differential equation. In their discrete form such equations reduce to systems of polynomial equations, known as algebraic equations. Using techniques from computational algebra one can address questions of the existence of solutions to such equations as well as the uniqueness of the solutions. The techniques are enumerative and exhaustive, requiring a finite number of computer operations. For example, calculating a bound to the total number of solutions reduces to computing the dimension of a linear vector space. The solution set itself may be constructed through the solution of an eigenvalue problem. The techniques are applied to a set of synthetic magnetotelluric values generated by conductivity variations within a layer. We find that the estimation of the conductivity and the electric field in the subsurface, based upon single-frequency magnetotelluric field values, is equivalent to a linear inverse problem. The techniques are also illustrated by an application to a magnetotelluric data set gathered at Battle Mountain, Nevada. Surface observations of the electric ( E _y ) and magnetic ( H _x ) fields are used to construct a model of subsurface electrical structure. Using techniques for algebraic equations it is shown that solutions exist, and that the set of solutions is finite. The total number of solutions is bounded above at 134 217 728. A numerical solution of the algebraic equations generates a conductivity structure in accordance with the current geological model for the area.     

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.     

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.     

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.     

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.     

Observation of electric currents in the Alaska oil pipeline resulting from auroral electrojet current sources

Summary. Currents in the 1.28 × 10³ km Alaska oil pipeline, induced from the ionospheric, auroral electrojet, were measured at three sites, near Fairbanks, Paxson and Valdez, Alaska, using a gradient configuration of two fluxgate magnetometers. The observed pipeline current magnitudes, which reached 50 A during times of mild geomagnetic activity, displayed a linear relationship with the electric earth potential. Using the induction relationship between the electric and magnetic fields and the typical spectral composition of the geomagnetic field at high latitudes, I obtained a spectral appearance of the current that shows a maximum in the range of 4.5- to 10-min period. Near Fairbanks the pipeline current amplitudes, I (Amperes), could be represented, approximately, by I = 0.65 B _x T ^−0.5, where B _x(nT) is the north—south geomagnetic field variation amplitude and T (min) is its apparent period. There is much less pipelines current at the sites south of Fairbanks. A previously established relationship between the local electric field and the planetary geomagnetic activity index, Ap , permitted a prediction of the pipeline current surge amplitudes in the Fairbanks region as approximately I = 5.0 Ap . Current surges larger than 500 A may be expected rather often in the Alaska pipeline during large geomagnetic storms.     

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.     

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.     

Review of the low-temperature magnetic properties of magnetite from a rock magnetic perspective

The low-temperature magnetic properties of magnetite are reviewed, and implications for rock magnetism considered. The behaviour of fundamental properties of magnetite at low temperatures near the Verwey transition ( T _v ) are documented, and attention is given to various Verwey transition theories. The low-temperature behaviour of the magnetic energies that control domain structure is reviewed in detail. For the first time in rock magnetic literature, the low-temperature anomaly in spontaneous magnetization ( M _s ) is documented and the differences between the saturation magnetization and M _s near the Verwey transition are discussed. It is argued that the low-temperature behaviour of the magnetocrystalline anisotropy, and in particular the anomaly at T _v , is most likely to affect multidomain remanence during low-temperature cycling. For multidomain crystals it is calculated that the large increase in magnetocrystalline anisotropy intensity and reduction in symmetry on cooling through T _v is likely to reduce the stability of closure domains.     

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.     

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.     

On the connection between the critical value of the "vertical" electrical anisotropy and earthquake prediction

Summary. It is shown that at any orientation of a dipole of direct current at the surface of a uniform anisotropic halfspace there exists a unique critical value η_cr of the "vertical" electric anisotropy η. For η > η_cr not all field lines are closed, and three straight field lines pass through the centre of the dipole. Also there exists a unique value η_eq > η_cr, where equatorial η is equal to zero. If η > η_eq then equatorial ρ_app is negative since the electric vector is turned more than through 90° at the equator of dipole, and it is shown that η_cr= 1.79, η_eq= 1.93 if the axis of dipole is directed symmetrically with respect to the axes of anisotropy. Routine observations of variations of ρ_app with the aim of earthquake prediction should be carried out taking into account the possibility of great fluctuations on η far from the epicentre of an earthquake.     

3-D seismic velocities calculated from lattice-preferred orientation and reflectivity of a lower crustal section: examples of the Val Sesia section (Ivrea zone, northern Italy)

To quantify the seismic properties of lower crustal rocks and to better constrain the origin of the lower crustal seismic reflectivity, we determined the complete 3-D seismic properties of a lower crustal section. Eight representative samples of the main lithologic and structural units outcropping in the Val Sesia (Ivrea zone) were studied in detail. The seismic velocities were calculated using the single crystal stiffness coefficients and the lattice preferred orientation (LPO) of each mineral in all samples. The 21 stiffness coefficients characterizing the elastic behaviour of each rock are determined. Mafic and ultramafic rocks such as pyroxenite and pyroxene-bearing gabbros display complex shear wave properties. These rocks are weakly birefringent (maximum 0.1 kms⁻¹) and it is difficult to find consistent relationships between the seismic properties and the rock structure. On the other hand, seismic properties of deformed felsic rocks are essentially controlled by mica. They display strong S -wave birefringence (0.3 km s⁻¹) and relatively high V _p anisotropy (7.6 per cent). Amphibole also strongly influences the rock birefringence patterns. For both kind of rocks, the foliation is highly birefringent and the fast polarized shear wave is systematically oriented parallel to the foliation. We show that the number of mineral phases in the rock strongly controls the anisotropy. The seismic anisotropy has a complex role in the P -wave reflectivity. Compared to the isotropic case, anisotropy enhances the reflection coefficient for about 60 per cent of the possible lithological interfaces. For 40 per cent of the interfaces, the reflection coefficient is much lower when one considers the medium as anisotropic.     

Rapid relaxation inversion of CSAMT data

In this paper an inversion algorithm for controlled-source audio frequency magnetotelluric data is presented. This algorithm combines 2.5-D finite element forward modelling with the concepts of rapid relaxation inversion of magnetotelluric data. The inversion uses the same technique to compute sensitivities as the rapid relaxation inversion, and these approximate sensitivities are validated by comparison with exact 2.5-D sensitivities. The comparison shows that the approximate sensitivities are similar in shape to the exact sensitivities when transmitter–receiver offsets are greater than one skin depth in the Earth. The magnitudes of the two sensitivities differ but the variations with depth are similar. Tests of the algorithm on synthetic data and field data provide promising results.     

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:      

Anisotropic reflection coefficients for a weak-contrast interface

In this article the interaction of plane waves with a weak-contrast interface between two weakly anisotropic half-spaces is investigated. The anisotropy dealt with is of a general type. The stress–displacement vectors of the plane waves are calculated by perturbation theory. By assuming that the jump in elastic parameters and density across the interface is small, one can derive a simple expression for the R _qPqP coefficient. In cases in which the wave motion is restricted to a symmetry plane of an anisotropic medium, simple expressions for the R _qSVqSV and R _SHSH coefficients are also derived.     

Gravity prediction using collocation and taking known mass density anomalies into account

Summary. The anomalous (gravitational) potential of the Earth, T , is split in two parts, T= T ^C + T _M. Here T _M is a harmonic function generated by known mass density anomalies and T ^C =T-T _M. This function will also be a harmonic function, which therefore may be approximated using the method of collocation, based on known gravity anomalies or altimeter derived geoid undulations, for example. Gravity anomalies can then be predicted using the known linearized relationship between T and Δ g . This procedure may give a 40–50 per cent increase in the precision of the prediction results as compared to a procedure where mass density anomalies are not taken into account.     

Error estimates in the measurement of anisotropic magnetic susceptibility

Hext's (1963 ) first-order analysis for error estimation in measurements of second-rank tensors is applied to the torque meter and the orthogonal-coil induction instrument, both representatives of the class of instruments that measures differences in magnetic susceptibility. The conventional measurement design for the torque meter is rotatable; an efficient and nearly rotatable design is proposed for the orthogonal-coil instrument. Error estimates are derived for a set of anisotropy parameters, K _mean , H and μ . The applicability of a first-order analysis is discussed and illustrated by simulation studies.     

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.     

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.