Spectral-finite-element approach to two-dimensional electromagnetic induction in a spherical earth

We present a spectral-finite-element approach to the 2-D forward problem for electromagnetic induction in a spherical earth. It represents an alternative to a variety of numerical methods for 2-D global electromagnetic modelling introduced recently (e.g. the perturbation expansion approach, the finite difference scheme). It may be used to estimate the effect of a possible axisymmetric structure of electrical conductivity of the mantle on surface observations, or it may serve as a tool for testing methods and codes for 3-D global electromagnetic modelling. The ultimate goal of these electromagnetic studies is to learn about the Earth's 3-D electrical structure.
Since the spectral-finite-element approach comes from the variational formulation, we formulate the 2-D electromagnetic induction problem in a variational sense. The boundary data used in this formulation consist of the horizontal components of the total magnetic intensity measured on the Earth's surface. In this the variational approach differs from other methods, which usually use spherical harmonic coefficients of external magnetic sources as input data. We verify the assumptions of the Lax-Milgram theorem and show that the variational solution exists and is unique. The spectral-finite-element approach then means that the problem is parametrized by spherical harmonics in the angular direction, whereas finite elements span the radial direction. The solution is searched for by the Galerkin method, which leads to the solving of a system of linear algebraic equations. The method and code have been tested for Everett & Schultz's (1995) model of two eccentrically nested spheres, and good agreement has been obtained.     

Slab low-velocity layer in the eastern Aleutian subduction zone

Local earthquakes in the vicinity of the Alaskan Peninsula's Shumagin Islands often produce arrivals between the main P and S arrivals not predicted by standard traveltime tables. Based on traveltime and polarization, these anomalous arrivals appear to be from P -to- S conversions at the surface of the subducted Pacific Plate beneath the recording stations. The P -to- S conversion occurs at the top of a low-velocity layer which extends to at least 150 km depth and is 8 ˜ 2 per cent slower than the overlying mantle. The slab is ˜ 7 per cent faster than the mantle. The low-velocity layer contains the foci of the earthquakes in the upper plane of the double seismic zone and confines PS ray paths to lie within it. These observations indicate that layered structures persist to positions well past the surface location of the volcanic front. Reactions forming high-pressure minerals do not yield slab-like velocities until beyond the point that subduction zone magma genesis occurs. If the subducted oceanic crust forms the layer, it is subducted essentially intact.     

Seismic hazard in central southern Africa

Seismic hazard maps of central-southern Africa where hazard has been expressed in terms of peak ground acceleration for an annual probability in excess of 10^-1 show relatively high values that distinguish the seismic hazard potential of the Deka fault zone, the mid-Zambezi basin-Luangwa rift and western central Mozambique. In areas such as central-southern Africa where little is known about the geology of the region and the fault systems have not been fully mapped, seismic hazard potential may be estimated from seismicity and broad-scale fault features. For this region, such potential is based on earthquake magnitude M_s ≥ 6. Events of such magnitude have recently occurred in the mid-Zambezi basin, southern Zimbabwe and western-central Mozambique. This paper follows the conventional probabilistic hazard analysis procedure, defining seismic source zones from seismicity based on instrumental records from a cataloque that spans a period of 83 years. Geological and geomorphological features in the region are described on the mesoscale and are correlated with the seismicity as broad fault zones. The scarcity of strong-motion accelerogram data necessitated the formulation of attenuation values based on random vibration theory (RVT).     

A quantitative methodology to extract regional magnetotelluric impedances and determine the dimension of the conductivity structure

A methodology to determine quantitatively the dimensionality of the dominant conducting structures and the resolution of the structural parameters in magnetotelluric data is presented. In addition, the method recovers the regional impedance responses when the regional structure can be characterized, at least approximately, as 1- or 2-D. The methodology is based upon three general models of the MT tensor, each of which has a distinct parameterization and physical interpretation. A weighted statistical residual describes quantitatively the fit of the model response to the data within the scatter of the measured data and hence permits: (1) tests of dimensionality, (2) determination of the appropriate strike angle, and (3) recovery of the regional responses.
The method has been tested extensively with synthetic data and proven to be successful. These synthetic studies give insight into the different physical parameterizations and the stability of the parameters determined. We describe and illustrate some of these synthetic studies. With field data, the methodology is not always as straightforward, but its application to a great many sites has proven valuable. Data from two closely spaced sites, which are both affected by strong but very different 3-D effects, are analysed to illustrate the geological significance of the results. The analyses reveal and recover regional responses within the data which indicate the presence of electrical anisotropy located deep in the crust and upper mantle. Analyses of the entire data set, of which these two form a part, confirm this finding.     

A comparison of the stability of stationary and non-stationary magnetotelluric analysis methods

A Magnetotelluric (MT) sounding was carried out at a site in south-east Queensland, in the Clarence-Moreton Basin. The synoptic recordings were taken over a period of four months at sampling frequencies from 500 Hz to 5 × 10-5 Hz. The resulting data was analysed by the stationary cross-frequency and the Cone kernel time-frequency distribution (TFD) methods of MT analysis. The results were compared as apparent resistivities on a daily basis for frequencies above 1 Hz, as well as over all the available data. The TFD MT apparent-resistivity results were more stable and less noisy on an daily basis than the cross-frequency results. Similarly the TFD analysis gave less noisy results than the cross-frequency analysis when all available data was processed. Application of these new non-stationary analysis techniques to MT processing should decrease the bias error problem of the MT methods and so increase reliability and repeatability of MT soundings.     

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.     

Two efficient algorithms for iterative linearized inversion of seismic waveform data

We present two equivalent algorithms for iterative linearized waveform inversion for 3-D Earth structure with respect to an arbitrary 3-D starting model; one is a matrix formulation, and the second is a wavefield formulation. Both algorithms require the computation of accurate synthetic seismograms, but neither requires that any particular method be used to compute the synthetics. The matrix formulation is equivalent to our previously published algorithm (Hara, Tsuboi & Geller 1991), but requires less than 10 per cent of the CPU time of the previous algorithm. The wavefield algorithm is equivalent to that of Tarantola (1986) and Mora (1987), but appears to be substantially more efficient.     

Microearthquake seismicity and fault-plane solutions in the southern Aegean and its geodynanic implications

Seismicity ( Ml <3.5) in the southern Aegean, located using data collected during seven weeks of recording by a temporary network of seismological stations, largely follows the Hellenic arc; the Sea of Crete is nearly aseismic, and only little activity is located south of the Hellenic trench, within the African plate. Focal mechanisms exhibit reverse faulting in the external part of the arc and normal faulting inside it. This normal faulting indicates N-S extension in the northern Aegean, the Gulf of Corinth, the Cyclades and Dodecanese Islands, but NW-SE extension in southern Peloponnese and western Crete and E-W extension in eastern Crete. This non-uniform strain pattern suggests that the Aegean region not only extends in a N-S sense, with the Hellenic arc moving south-westward relative to the Eurasian plate, but also by E-W extension of its southern margin, so that there is a net divergence of material.     

Global mapping of upper mantle reflectors from long-period SS precursors

Long-period precursors to SS resulting from underside reflections off upper mantle discontinuities ( SdS where d is the discontinuity depth) can be used to map the global distribution and depth of these reflectors. We analyse 5,884 long-period seismograms from the Global Digital Seismograph Network (1976-1987, shallow sources, transverse component) in order to identify SdS arrivals. Corrections for velocity dispersion, topography and crustal thickness at the SS bounce point, and lateral variation in mantle velocity are critical for obtaining accurate estimates of discontinuity depths. The 410 and 660 km discontinuities are observed at average depths of 413 and 653 km, and exhibit large-scale coherent patterns of topography with depth variations up to 40 km. These patterns are roughly correlated with recent tomographic models, with fast anomalies in the transition zone associated with highs in the 410 km discontinuity and lows in the 660 km discontinuity, a result consistent with laboratory measurements of Clapeyron slopes for the appropriate phase changes. The best resolved feature in these maps is a trough in the 660 km discontinuity in the northwest Pacific, which appears to be associated with the subduction zones in this region. Amplitude variations in SdS arrivals are not correlated with discontinuity depths and probably result from focusing and defocusing effects along the ray paths. The SdS arrivals suggest the presence of regional reflectors in the upper mantle above 400 km. However, only the strongest of these features are above probable noise levels due to sampling inadequacies.