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.     

Three-dimensional induction in a non-uniform thin sheet at the surface of a uniformly conducting earth

Summary. A new method for solving problems in three-dimensional electromagnetic induction in which the Earth is represented by a uniformly conducting half-space overlain by a surface layer of variable conductance is presented. Unlike previous treatments of this type of problem the method does not require the fields to be separated into their normal and anomalous parts, nor is it necessary to assume that the anomalous region is surrounded by a uniform structure; the model may approach either an E- or a B -polarization configuration at infinity. The solution is expressed as a vector integral equation in the horizontal electric field at the surface. The kernel of the integral is a Green's tensor which is expressed in terms of elementary functions that are independent of the conductance. The method is applied to an illustrative model representing an island near a bent coastline which extends to infinity in perpendicular directions.     

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.     

Interaction of a compressional impulse with a slot normal to the surface of an elastic half space – II

Summary. An improved finite difference scheme has been used to simulate the propagation of a plane P -impulse in an elastic half space with a slot normal to its surface. Various angles of incidence and dimensions of slot are considered. The numerical results are presented in several visualizations; each emphasized a different type of wave and all representations help in understanding the scattered and diffracted wave pattern. Experiments were carried out using 0.5—6 MHz ultrasonic pulses on a duralumin semidisc with a surface-breaking slot and the results are compared with those given by the numerical models.
The scattered wavefield includes compressional and Rayleigh pulses whose amplitude increases at the front of the slot and decreases behind it, as the angle of incidence is reduced. A diffracted compressional pulse is generated with a semicircular wavefront centred at the mid-point of the bottom of the slot. Also, two elliptical eddies are excited at the lower corners. These ellipses propagate into the medium and eventually spread out to form arc-shaped shear pulses.
In the shadow zone, behind the slot, the two components of displacement show independent behaviour. The horizontal component decreases either with decreasing angle of incidence or as the slot is made deeper. For acute angles, a reduction of displacement amplitude of about 50 per cent is obtained when the depth of the slot is made a half pulse width. On the other hand there is no diminution of the vertical displacement behind the slot, and, near the upper right corner, it is even amplified.     

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.     

Kirchhoff–Helmholtz reflection seismograms in a laterally inhomogeneous multi-layered elastic medium – II. Computations

Summary. High-frequency reflection and refraction seismograms for laterally variable multi-layered elastic media are computed by using the frequency domain elastic Kirchhoff–Helmholtz (KH) theory of Frazer and Sen. Both source and receiver wavefields are expanded in series of generalized rays and then elastic (KH) theory is applied to determine the coupling between each source ray and each receiver ray at each interface. The motion at the receiver is given as a series of integrals, one for each generalized ray. We use geometrical optics and plane wave reflection and transmission coefficients for rapid evaluation of the integrand. When the source or the receiver ray field has caustics on the surface of integration geometrical ray theory breaks down and this gives rise to singularities in the KH integrand. We repair this using methods suggested by Frazer and Sen.
Examples of reflection seismograms for 2-D structures computed by elastic KH theory are shown. Those for a vertical fault scarp structure are compared with the seismograms obtained by physical modelling. Then OBS data obtained from the mid-America trench offshore Guatemala area are analysed by computing KH synthetics for a velocity model that has been proposed for that area. Our analysis indicates the existence of a small low-velocity zone off the trench axis.
No head wave arrivals are obtained in our KH synthetics since we do not consider multiple interactions of a ray with an interface. The nearly discontinuous behaviour of elastic R/T coefficients near the critical angle causes small spurious phases which arrive later than the correct arrivals.     

Three-dimensional massively parallel electromagnetic inversion—II. Analysis of a crosswell electromagnetic experiment

A 3-D electromagnetic inversion scheme has been successfully applied to crosswell electromagnetic data collected at the Richmond Field Station near Berkeley, California. By comparing images of data collected before and after the injection of 50000 gallons of salt water, a 3-D image of the plume has been developed, which shows the location of zones of maximum permeability surrounding the injection well through which the salt water has migrated. A resolution analysis has determined that the location of the plume is fairly accurate. However, the image of the geology will be distorted due to incomplete data coverage. This latter problem is further complicated by the fact that data residuals indicate that the wells deviate from the vertical. These deviations have been shown to cause artefacts within the images and thus further reduce the accuracy of the images with respect to the geology.     

Crustal structure across the Grand Banks–Newfoundland Basin Continental Margin – II. Results from a seismic reflection profile

New multichannel seismic reflection data were collected over a 565 km transect covering the non-volcanic rifted margin of the central eastern Grand Banks and the Newfoundland Basin in the northwestern Atlantic. Three major crustal zones are interpreted from west to east over the seaward 350 km of the profile: (1) continental crust; (2) transitional basement and (3) oceanic crust. Continental crust thins over a wide zone (∼160 km) by forming a large rift basin (Carson Basin) and seaward fault block, together with a series of smaller fault blocks eastwards beneath the Salar and Newfoundland basins. Analysis of selected previous reflection profiles (Lithoprobe 85-4, 85-2 and Conrad NB-1) indicates that prominent landward-dipping reflections observed under the continental slope are a regional phenomenon. They define the landward edge of a deep serpentinized mantle layer, which underlies both extended continental crust and transitional basement. The 80-km-wide transitional basement is defined landwards by a basement high that may consist of serpentinized peridotite and seawards by a pair of basement highs of unknown crustal origin. Flat and unreflective transitional basement most likely is exhumed, serpentinized mantle, although our results do not exclude the possibility of anomalously thinned oceanic crust. A Moho reflection below interpreted oceanic crust is first observed landwards of magnetic anomaly M4, 230 km from the shelf break. Extrapolation of ages from chron M0 to the edge of interpreted oceanic crust suggests that the onset of seafloor spreading was ∼138 Ma (Valanginian) in the south (southern Newfoundland Basin) to ∼125 Ma (Barremian–Aptian boundary) in the north (Flemish Cap), comparable to those proposed for the conjugate margins.     

Shear-wave polarizations near the North Anatolian Fault – II. Interpretation in terms of crack-induced anisotropy

Summary. The polarizations of shear waves recorded by networks of digital three-component seismometers immediately above small earthquakes near the North Anatolian Fault in Turkey display shear-wave splitting on almost all shear-wave seismograms recorded within the shear-wave window. This splitting is incompatible with source radiation-patterns propagating through simple isotropic structures but is compatible with effective anisotropy of the internal structure of the rock along the ray paths. This paper interprets the phenomena in terms of widespread crack-induced anisotropy. Distributions of stress-induced cracks model many features of the observations, and synthetic polarization diagrams calculated for propagation through simulated cracked rock are similar to the observed patterns. This evidence for widespread crack-induced anisotropy lends strong support to the hypothesis of extensive-dilatancy anisotropy (EDA) suggested by laboratory experiments in subcritical crack-growth. The crucial evidence confirming some form of EDA would be observations of temporal changes in shear-wave splitting as the stress field alters the crack density and crack geometry. There is some weak evidence for such temporal changes at one site, but further analysis of suitable digital three-component seismometer networks in seismic areas is required to confirm EDA.