B-polarization induction in two generalized thin sheets at the surface of a conducting half-space

Summary. A new closed-form solution is obtained analytically for a B- polarization induction problem of geophysical interest, in which a local region of the Earth is represented by a generalized thin sheet at the surface of and in electrical contact with a uniformly conducting half-space. The generalized sheet, first introduced by Ranganayaki & Madden, is a mathematical idealization of a double layer which consists, in this problem, of two adjacent half-planes with distinct conductances representing a surface conductivity discontinuity such as an ocean—coast boundary, underlain by a uniform sheet of finite integrated resistivity representing the lower crust. The resistive sheet exerts a considerable mathematical influence on the solution causing, under certain conditions, an additional pole to appear in one of the forms of contour integral by which the solution can be expressed; it also weakens or eliminates field singularities that would otherwise occur at the conductance discontinuity. A numerical calculation is made for model parameters typifying an ocean—coast boundary underlain by a highly resistive crust. It is found that the residue of the pole associated with the resistive sheet dominates the solution for this example, the main consequence of which is a huge increase in the horizontal range over which the induced currents adjust themselves between the different 'skin-effect' distributions at infinity on either side of the model. Moreover the solution shows that this 'adjustment distance' has a more complicated dependence on the conductance and integrated resistivity of the sheet than that given simply by the square root of their product which was the length parameter proposed by Ranganayaki & Madden.     

Two-dimensional induction in a thin sheet of variable integrated conductivity at the surface of a uniformly conducting earth

Summary. The forward solution of the general two-dimensional problem of induction in a model earth comprising a uniformly conducting half-space covered by a thin sheet of variable integrated conductivity is obtained. Unlike some previous treatments of similar problems, the method presented here does not require the field to be separated into its normal and anomalous parts. Both the E - and B -polarization modes of induction are considered and in each case the solution is expressed in terms of the horizontal component of the electric field satisfying, on the surface of the conductor, a singular integral equation whose kernel is a well-known analytic function. A recently published solution of the coast effect is included as a special case. The numerical procedure for solving the integral equations is described and some illustrative calculations are presented.     

Electric field at the seafloor due to a two-dimensional ionospheric current

The relation between the seafloor electric field and the surface magnetic field is studied. It is assumed that the fields are created by a 2-D ionospheric current distribution resulting in the E-polarization. The layered earth below the sea water is characterized by a surface impedance. The electric field at the seafloor can be expressed either as an inverse Fourier transform integral over the wavenumber or as a spatial convolution integral. In both integrals the surface magnetic field is multiplied by a function that depends on the depth and conductivity of the sea water and on the properties of the basement. The fact that surface magnetic data are usually available on land, not at the sea surface, is also considered. Test computations demonstrate that the numerical inaccuracies involved in the convolution method are negligible. The theoretical equations are applied to calculate the seafloor electric fields due to an ionospheric line current or associated with real magnetic data collected by the IMAGE magnetometer array in northern Europe. Two different sea depths are considered: 100 m (the continental shelf) and 5 km (the deep ocean). It is seen that the dependence of the electric field on the oscillation period is weaker in the 5 km case than for 100 m.     

Inverse scattering of surface waves: imaging of near-surface heterogeneities

An efficient inverse scattering method is developed for imaging near-surface heterogeneities using scattered surface waves. Three dimensional elastodynamic wave propagation and scattering in a laterally invariant embedding medium is considered. The Born Approximation is used and the scattered wavefield is expressed as a domain type integral representation. The computation time of Green's tensor elements is reduced by considering the radial symmetry of the medium. The method is validated by numerical tests. Ultrasonic laboratory data obtained from a scale model experiment are used for imaging the near-surface inhomogeneities caused by an epoxy-filled hole in the surface of an aluminum block. Both synthetic and the scale model tests show that the location, the actual density contrast and the depth of the inhomogeneities are reasonably well estimated.     

Love waves excited by discontinuous propagation of a rupture front

Summary. Analytical results are presented for Love waves generated by sudden changes of the rate of advance of a curved rupture front in an inclined fault plane that is embedded in an elastic half-space. The boundary condition at the surface of the half-space approximates the presence of an overlying layer. The calculation consists of two parts. First, ray theory is used to calculate far-field approximations to the horizontally polarized wavefields which are emitted when the speed of the rupture front suddenly changes. These fields can be expressed as products of emission coefficients (which govern the angular dependence) and propagation terms. Secondly, a representation integral for the Love wave over a surface enclosing the rupture front is constructed, using the emitted signal and an appropriate Green's function. This integral is evaluated asymptotically. The resulting approximate Love-wave spectrum shows an explicit dependence on the nature of the rupture process, on the rupture-front and fault-plane geometry, and on the magnitude of a sudden change in the rate of advance of the rupture front.     

Induction in a thin sheet of variable conductance at the surface of a stratified earth – II. Three-dimensional theory

Summary. The algorithm of Dawson & Weaver for modelling electromagnetic induction effects in a thin sheet at the surface of a uniform earth is modified to permit the use of a layered earth model. The theory is developed in Fourier space in terms of the toroidal and poloidal transfer functions instead of with the Green's function approach which was used by Dawson & Weaver. The integral equation for the surface electric field and most of the integral formulae for the derived field components are the same as before, except for the inclusion of additional integral the kernel of which has to be calculated numerically with the aid of fast Hankel transforms. The accuracy of the results is tested by comparing solutions with those obtained from a related 2-D algorithm and finally an example of 3-D modelling is presented.     

Exact synthetic seismograms for an inhomogeneity in a layered elastic half-space

Summary. The propagation of a pulsed elastic wave in the following geometry is considered. An elastic half-space has a surface layer of a different material and the layer furthermore contains a bounded 3-D inhomogeneity. The exciting source is an explosion, modelled as an isotropic pressure point source with Gaussian behaviour in time.
The time-harmonic problem is solved using the null field approach (the T matrix method), and a frequency integral then gives the time-domain response. The main tools of the null field approach are integral representations containing the free space Green's dyadic, expansions in plane and spherical vector wave functions, and transformations between plane and spherical vector wave functions. It should be noted that the null field approach gives the solution to the full elastodynamic equations with, in principle, an arbitrarily high accuracy. Thus no ray approximations or the like are used. The main numerical limitation is that only low and intermediate frequencies, in the sense that the diameter of the inhomogeneity can only be a few wavelengths, can be considered.
The numerical examples show synthetic seismograms consisting of data from 15 observation points at increasing distances from the source. The normal component of the velocity field is computed and the anomalous field due to the inhomogeneity is sometimes shown separately. The shape of the inhomogeneity, the location and depth of the source, and the material parameters are all varied to illustrate the relative importance of the various parameters. Several specific wave types can be identified in the seismograms: Rayleigh waves, direct and reflected P -waves, and head waves.     

Bias estimation and correction for triangle-based surface area calculations

The calculation of surface area is meaningful for a variety of space-filling phenomena, e.g., the packing of plants or animals within an area of land. With Digital Elevation Model (DEM) data we can calculate the surface area by using a continuous surface model, such as by the Triangulated Irregular Network (TIN). However, just as the triangle-based surface area discussed in this paper, the surface area is generally biased because it is a nonlinear mapping about the DEM data which contain measurement errors. To reduce the bias in the surface area, we propose a second-order bias correction by applying nonlinear error propagation to the triangle-based surface area. This process reveals that the random errors in the DEM data result in a bias in the triangle-based surface area while the systematic errors in the DEM data can be reduced by using the height differences. The bias is theoretically given by a probability integral which can be approximated by numerical approaches including the numerical integral and the Monte Carlo method; but these approaches need a theoretical distribution assumption about the DEM measurement errors, and have a very high computational cost. In most cases, we only have variance information on the measurement errors; thus, a bias estimation based on nonlinear error propagation is proposed. Based on the second-order bias estimation proposed, the variance of the surface area can be improved immediately by removing the bias from the original variance estimation. The main results are verified by the Monte Carlo method and by the numerical integral. They show that an unbiased surface area can be obtained by removing the proposed bias estimation from the triangle-based surface area originally calculated from the DEM data.     

Induction in a thin sheet of variable conductance at the surface of a stratified earth - I. Two-dimensional theory

Summary. An existing 2-D integral equation method for modelling electromagnetic induction in a thin sheet at the surface of a uniform half-space can be generalized to deal with a layered half-space by the inclusion of an extra term in the integral equation. The results obtained are shown to be in excellent agreement with finite difference solutions to the same modelling problem.     

Electromagnetic response of an ocean-coast model to E-polarization induction

summary . An ocean-coast model which consists of a uniformly conducting half-space screened by a perfectly conducting half-plane (the model ocean) is studied. On the land the electric field decreases continuously to zero as the coast is approached. The horizontal magnetic field component is found to vary rapidly, but remains finite; the vertical component on the other hand, increases to infinity at the coast. On the surface of the model ocean as well as on the sea floor, electric field and vertical magnetic field are both nil, but the horizontal magnetic field becomes singular as the seashore is approached. This horizontal magnetic field however, is different on the sea floor and at the ocean surface, because the integrated ocean current is finite, even growing to infinity as the shore is approached. The very large ocean currents near the shore act as an extremely long line antenna, which radiates far afield. This antenna feature explains the very long range of the ocean-coast effects observed under E -polarization induction, compared to the corresponding H -polarization effects where no such antenna-like feature occurs. A similarly large difference of ranges can be expected for all shallow structures with large lateral conductivity contrasts. The present study may therefore be of some interest in relation to geomagnetic depth soundings by the inductive and magnetotelluric methods, as well as in understanding the ocean-coast effect known for some time from records of coastal observatories.     

Application of the boundary element method to a terrestrial heat flow problem

b
In this paper, we apply the boundary element method (BEM) to the 2-D steady state heat flow problem of what would be the perturbation to the regional temperature gradient, and hence heat flow density, determined from temperatures measured in a borehole that passes close to, but does not penetrate, a body of anomalous thermal properties. This type of problem with an infinity boundary is particularly well suited to the BEM.
The results have been compared with those obtained from analytical solutions for bodies of simple shape; it is found that for the worst case of a close approach to a boundary of small radius of curvature, a numerical modelling error of less than 1 per cent can still be obtained provided the length of each element is less than the shortest distance between the calculation point and the object.     

The theory of finite frequency body wave synthetic seismograms in inhomogeneous elastic media

b
A new method is presented by means of which one can compute finite frequency synthetic seismograms for media whose velocity and density are continuous functions of two or three spatial variables. Basically, the method is a generalization of the familiar phase integral method, to which it reduces in a stratified medium. For a given source location the travel-time and distance functions needed to compute synthetics are obtained by numerically tracing rays through the model. This information is then used to evaluate a double integral over frequency and take-off angle at the source. The solution obtained reduces to the geometrical optics solution wherever that is valid but it also works in shadows and at caustics without knowing explicitly where these may be located. The method can be used as a spectral method, in which the integral over take-off angle is evaluated first, or as a slowness method, in which the frequency integral is evaluated first.     

Scattering of surface waves modelled by the integral equation method

The integral equation method is used to model the propagation of surface waves in 3-D structures. The wavefield is represented by the Fredholm integral equation, and the scattered surface waves are calculated by solving the integral equation numerically. The integration of the Green's function elements is given analytically by treating the singularity of the Hankel function at   R = 0  , based on the proper expression of the Green's function and the addition theorem of the Hankel function. No far-field and Born approximation is made. We investigate the scattering of surface waves propagating in layered reference models imbedding a heterogeneity with different density, as well as Lamé constant contrasts, both in frequency and time domains, for incident plane waves and point sources.     

Bimodal electromagnetic induction in non-uniform thin sheets with an application to the northern Pyrenean induction anomaly

In order to handle the distortion of large-scale induced electric currents by local conductivity anomalies, the problem of electromagnetic induction in non-uniform thin sheets has been reformulated in terms of an integral equation over the anomalous domain. This formulation considers in the layered substratum in addition to toroidal currents also the poloidal current mode (vertical current loops), at the expense that two scalar functions have to be determined. Simple formulas for the required kernels are derived. The algorithm is applied to model the gross features of the northern Pyrenean induction anomaly. It is suggested that this pronounced anomaly results from a conductive channel between the Atlantic Ocean and the Mediterranean Sea.     

The complex-image method for calculating the magnetic and electric fields produced at the surface of the Earth by the auroral electrojet

For studying the auroral electrojet and for examining the effects it can produce in power systems on the ground, it is useful to be able to calculate the magnetic and electric fields that the electrojet produces at the surface of the Earth. Including the effects of currents induced in the Earth leads to a set of integral expressions, the numerical computation of which is complicated and demanding of computer resources. An approximate solution can be achieved by representing the induced currents by an image current at a complex depth. We present a simple derivation of the complex-image expressions and use them to calculate the fields produced by the auroral electrojet at the surface of an earth represented by layered conductivity models. Comparison of these results with ones obtained using the exact integral solution show that the errors introduced are insignificant compared to the uncertainties in the parameters used. The complex-image method thus provides a simple, fast and accurate means of calculating the magnetic and electric fields.     

Developmetn stage threshold of watershed Landforms in Loess Plateau and separation of erosion mechanism

Threshold is a limit and marginal point when a qualitative change takes place.Since threshold principle was firstly used in landform research by S.A.Schumm in the 1970s,the quantitative research on watershed landform development stage has come true.Davis proposed a three-stage model on landform evolution:young ,mature and old stages.Thereafter Stahler quantified this model by hypsometric analysis method.The authors thought that the material movement stage cannot be expressed by hypsometric method in watershed landform at development stage,because of the uncertainty on stage delimitaion.To meet this shortcoming,this paper presents an integral erosion value method.A clear delimitation on landform development stage in the Loess Plateau region has been tested by this method.The result shows that gullied loessial hilly area is at the mature stage,and gullied loessial tableland area is at the young stage.It is estimated that from the point of erosion related sediment yield,natural erosion accounts for 70% of the total erosion amount,and artificial accelerating erosion is 30%.Therefore soil and water conservation is very crucial for the Loess Plateau.     

Modelling of electromagnetic fields in thin heterogeneous layers with application to field generation by volcanoes—theory and example

When interpreting electromagnetic fields observed at the Earth's surface in a realistic geophysical environment it is often necessary to pay special attention to the effects caused by inhomogeneities of the subsurface sedimentary and/or water layer and by inhomogeneities of the Earth's crust. The inhomogeneities of the Earth's crust are expected to be especially important when the electromagnetic field is generated by a source located in a magma chamber of a volcano. The simulation of such effects can be carried out using generalized thin-sheet models, which were independently introduced by Dmitriev (1969 ) and Ranganayaki & Madden (1980 ). In the first part of the paper, a system of integral equations is derived for the horizontal current that flows in the subsurface inhomogeneous conductive layer and for the vertical current crossing the inhomogeneous resistive layer representing the Earth's mantle. The terms relating to the finite thickness of the laterally inhomogeneous part of the model are retained in the equations. This only marginally complicates the equations, whilst allowing for a significant expansion of the approximation limits.
  The system of integral equations is solved using the iterative dissipative method developed by the authors in the period from 1978 to 1988. The method can be applied to the simulation of the electromagnetic field in an arbitrary inhomogeneous medium that dissipates the electromagnetic energy. When considered on a finite numerical grid, the integral equations are reduced to a system of linear equations that possess the same contraction properties as the original equations. As a result, the rate at which the iterative-perturbation sequence converges to the solution remains independent of the numerical grid used for the calculations. In contrast to previous publications on the method, aspects of the algorithm implementation that guarantee its effectiveness and robustness are discussed here.     

An Integral Equation and its Solution for some Two-and Three-Dimensional Problems in Resistivity and Induced Polarization

Summary. The potential, U , about a point electrode, at the surface of a layered ground in which there is an heterogeneity embedded, satisfies the integral equation:
Here, U * and σ* are the corresponding quantities for the potential and conductivity without the heterogeneity. The integral is taken over the surface of the heterogeneity, ∂ U /∂ n is the normal derivative (in the direction of the outward normal) of U , and G is a Green's function.
Solutions to this equation can readily be found by using the Galerkin method of solving integral equations. The solutions of this equation when the heterogeneity is a sphere or a cylinder in a uniform ground or beneath a conductive overburden are the most readily found.
When the solution of the integral has been found for the potential it is a simple matter to calculate the apparent resistivity or chargeability for any electrode configuration.     

Calculation of electromagnetic sensitivities in the time domain

The speed of calculating sensitivities for 3-D conductivity structures for time- domain electromagnetic methods is significantly improved by applying the reciprocity theorem directly in the time domain. The sensitivities are obtained by convolving the electric field in the subsurface due to a transmitter at the surface with the electric field impulse response due to another transmitter, which replaces the original receiver. The acceleration compared to the classical perturbation method is approximately P/R , where P is the number of model parameters and R is the number of receiver positions. If the sensitivity has to be calculated very close to the receiver, approximate sensitivities can be obtained using an integral condition. Comparisons with the classical perturbation approach show that the method gives accurate results. Examples using transmitter–receiver configurations from a long-offset transient electromagnetics survey demonstrate the usefulness of sensitivities for the evaluation of resolution properties.     

Development stage threshold of watershed landforms in Loess Plateau and separation of erosion mechanism

Threshold is a limit and marginal point when a qualitative change takes place. Since threshold principle was firstly used in landform research by S.A.Schumm in the 1970s, the quantitative research on watershed landform development stage has come true. Davis proposed a three-stage model on landform evolution: young, mature and old stages. Thereafter Strahler quantified this model by hypsometric analysis method. The authors thought that the material movement stage cannot be expressed by hypsometric method in watershed landform at development stage, because of the uncertainty on stage delimitation. To meet this shortcoming, this paper presents an integral erosion value method. A clear delimitation on landform development stage in the Loess Plateau region has been tested by this method. The result shows that gullied loessial hilly area is at the mature stage, and gullied loessial tableland area is at the young stage. It is estimated that from the point of erosion related sediment yield, natural erosion accounts for 70% of the total erosion amount, and artificial accelerating erosion is 30%. Therefore soil and water conservation is very crucial for the Loess Plateau.