The magnetic structure of convection-driven numerical dynamos

The generation of a magnetic field in numerical simulations of the geodynamo is an intrinsically 3-D and time-dependent phenomenon. The concept of magnetic field lines and the frozen-flux approximation can provide insight into such systems, but a suitable visualization method is required. This paper presents results obtained using the Dynamical Magnetic Field line Imaging (DMFI) technique, which is a representation of magnetic field lines accounting for their local magnetic energy, together with an algorithm for the time evolution of their anchor points. The DMFI illustrations are consistent with previously published dynamo mechanisms, and allow further investigation of spatially and temporally complex systems. We highlight three types of magnetic structures: (i) magnetic cyclones and (ii) magnetic anticyclones are expelled by, but corotate with axial flow vortices; (iii) magnetic upwellings are amplified by stretching and advection within flow upwellings, and show structural similarity with helical plumes found in rotating hydrodynamic experiments. While magnetic anticyclones are responsible for the regeneration of a stable axial dipole, here we show that excursions and reversals of the dipole axis are caused by the emergence of magnetic upwellings, which amplify and transport a generally multipolar magnetic field from the inner to the outer boundary of the models. Geomagnetic observations suggest the presence of magnetic structures similar to those found in our models; thus, we discuss how our results may pertain to Earth's core dynamo processes. In order to make DMFI a standard tool for numerical dynamo studies, a public software package is available upon request to the authors (supplementary material is available at: http://www.ipgp.jussieu.fr/~aubert/DMFI.html ).     

The nature of particle motion in regional seismograms and its utilization for phase identification

The particle motion of regional arrivals is frequently treated in automatic phase-recognition schemes as that appropriate to simple P or S waves incident on an elastic, laterally homogeneous half-space. This model implies that the motion in ' P -type' phases can be described in terms of a single, generalized signal process and ' S -type' phases in terms of two independent processes ( SV and SH ) and thus, all regional arrivals could be fully characterized by three components of motion. In this paper, we present anlyses of the particle-motion patterns of various regional arrivals recorded at the ARCESS array from closely spaced events in the Kola Peninsula. We have found that only Pn -particle motion, described in terms of two independent signal processes, can be reliably characterized by three-component recordings. On the other hand, the various regional arrivals following Pn , such as Pg and Sn and Lg , can only be poorly characterized on the basis of three-component recordings alone. The reason is that these arrivals must be described in terms of more than two independent generalized signal processes, at least three for Pg and Sn , and possibly up to five for Lg . Recognition of these phases will thus require the use of more sensors than signal processes in the observing sensor configuration, such as three-component sensors combined with a small tripartite array. We have investigated the feasibility of adaptive, automatic recognition of regional arrivials by a wavefield extrapolation scheme utilizing such a mini-array. The process, which appear to be promising, adaptively learns the particle-motion patterns of individual arrivals, including complex site-response functions, from examples of closely located regional events.     

On the frequency contents of local events: source or path effect?

Local events which occurred in the Timanfaya volcanic field are analysed and their spectral contents are interpreted as being due to the effect of a finely layered medium with large impedance contrast, the presence of which is well known from well-logging data. By considering the layered medium as a filter, its spectral response (amplitude and phase) is computed from different randomly generated values of its parameters (total thickness, mean thickness of the individual layers and impedance contrast), showing that this stratigraphic filter acts as a low-pass filter. This effect can be correctly modelled by means of an exponential correlation function. However, in some cases resonant peaks appear in the spectrum superimposed on the low-pass filter effect, which can only be attributed to resonant scattering. An analysis of the phase and group delay associated with the amplitudes is made and strongly suggests that it is necessary to incorporate the theory of resonant scattering for the retrieval of medium properties. From the point of view of the amplitudes, the superposition of both effects, low-pass filter and resonant peaks, allows us to reproduce the gross features that are usually attributed to source effects. The need for a careful analysis of the local site response previous to any interpretation of the spectral features in terms of source characteristics is thus emphasized.     

The Aegean region: deep structures and seismological properties

3-D P -wave velocity anomaly patterns, the geometry of the Hellenic Wadati-Benioff zone and of the seismically active fracture zones in the overlying continental wedge, and the depth distribution of seismogenic properties (seismic activity, energy and b value) for the Aegean region are examined in this study. The western and eastern parts of the area studied differ significantly in terms of the depth distribution of the considered seismological parameters and the velocity structures. The results indicate different physical conditions in the western and eastern parts of the region in question. Data for the surface heat flow and the evolution of the volcanism in the Aegean region since the early Eocene are employed to interpret the results of the present study. Possible geodynamic processes at the plate boundaries between Africa and Eurasia are discussed.     

Late Quaternary kinematics of the Pallatanga strike-slip fault (Central Ecuador) from topographic measurements of displaced morphological features

The northeast-trending Pallatanga right-lateral strike-slip fault runs across the Western Cordillera connecting N50E-N70E trending normal faults in the Gulf of Guayaquil with N-S reverse faults in the Interandean Depression. Over most of its length, the fault trace has been partly obscured by erosional processes and can be inferred in the topography only at the large scale. Only the northern fault segment, which follows the upper Rio Pangor valley at elevations above 3600 m, is prominent in the morphology. Valleys and ridges cut and offset by the fault provide an outstanding record of right-lateral cumulative fault displacement. The fault geometry and kinematics of this particular fault segment can be determined from detailed topographic levellings. The fault strikes N30E and dips 75 to the NW. Depending on their size and nature, transverse morphological features such as tributaries of the Rio Pangor and intervening ridges, reveal right-lateral offsets which cluster around 27 ± 11m, 41.5 ± 4 m, 590 ± 65 m and 960 ± 70 m. The slip vector deduced from the short-term offsets shows a slight reverse component with a pitch of about 11.5 SW. The 41.5 ± 4 m displacements are assumed to be coeval with the last glacial termination, yielding a mean Holocene slip-rate of 2.9- 4.6 mm yr⁻¹. Assuming a uniform slip rate on the fault in the long term, the 27 m offset appears to correlate with an identified middle Holocene morphoclimatic event, and the long term offsets of 590 m and 960 m coincide with the glacial terminations at the beginning of the last two interglacial periods.     

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.     

Observations of geomagnetic pulsations and variations with a new borehole magnetometer down to depths of 3000 m

b
A triaxial magnetometer has been developed for investigating the in situ skin effect of horizontal geomagnetic pulsations and variations in deep boreholes. The observations were carried out in the pilot borehole of the Continental Deep Drilling Program of Germany (KTB) down to depths of 3000 m and up to temperature of 90 C. A weak skin effect, due to the known very low conductivity of the penetrated crystalline rocks, of 90 to 95 per cent in amplitude and -5 to -10 rotation in phase has been observed at periods of 10 s and magnetometer depth of 2400 m.
An essential prerequisite for all calculations is the accurate determination of the orientation of the downhole magnetometer. It is demonstrated how oriented record samples of temporal variations at depth correlate precisely with those from the surface.
Results from surface magnetotelluric (MT) investigations show strong local distortions of the telluric field. The distortion of the MT tensor response has been determined by means of newly introduced skin-effect transfer functions, which are assumed to be undistorted.     

Glacial rebound of the British Isles—II. A high-resolution, high-precision model

Observations of ice movements across the British Isles and of sea-level changes around the shorelines during Late Devensian time (after about 25 000 yr BP) have been used to establish a high spatial and temporal resolution model for the rebound of Great Britain and associated sea-level change. The sea-level observations include sites within the margins of the former ice sheet as well as observations outside the glaciated regions such that it has been possible to separate unknown earth model parameters from some ice-sheet model parameters in the inversion of the glacio-hydro-isostatic equations. The mantle viscosity profile is approximated by a number of radially symmetric layers representing the lithosphere, the upper mantle as two layers from the base of the lithosphere to the phase transition boundary at 400 km, the transition zone down to 670 km depth, and the lower mantle. No evidence is found to support a strong layering in viscosity above 670 km other than the high-viscosity lithospheric layer. Models with a low-viscosity zone in the upper mantle or models with a marked higher viscosity in the transition zone are less satisfactory than models in which the viscosity is constant from the base of the lithosphere to the 670 km boundary. In contrast, a marked increase in viscosity is required across this latter boundary. The optimum effective parameters for the mantle beneath Great Britain are: a lithospheric thickness of about 65 km, a mantle viscosity above 670 km of about (4-5) 10²⁰ Pa s, and a viscosity below 670 km greater than 4 × 10²¹ Pa s.     

Reciprocity theorem to compute the static deformation due to a point dislocation buried in a spherically symmetric earth

Intriguing reciprocity relations exist between the static deformation excited by a point dislocation in a SNREI earth and those generated by external forces, such as tidal force, surface loading and surface shear forces. Coseismic deformations can be rewritten as follows: (1) potential change in terms of the tide deformation field, (2) radial displacement in terms of the load and tidal deformation fields, and (3) tangential displacement in terms of shear and torsional deformation fields. The relations greatly reduce the effort to compute the coseismic crustal deformation in a spherically symmetric earth.     

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.     

Wavepath traveltime tomography

The elastic-wave equation is used to construct sensitivity kernels relating perturbations in elastic parameters to traveltime deviations. Computation of the functions requires a correlation of the forward-propagating seismic wavefield with a backward propagation of the residual wavefield. The computation of the wavefields is accomplished using a finite difference algorithm and is efficiently executed on a CM-2 parallel processor. The source and receiver locations have maximum sensitivity to velocity structure. The sensitivity kernels or wavepaths are well suited for transmission traveltime inversion such as cross-borehole tomography and vertical seismic profiling. Conventional ray tomography and wavepath tomography are applied to a set of P -wave arrival times, from a cross-borehole experiment at Kesterson, California. Because the wavepaths have increased sensitivity near the source and receiver there are differences in resolution of the velocity structure. Both techniques recover the same relative variations in velocity where the coverage is adequate. The wavepath solution is more laterally continuous and the dominant variation is vertical, as is expected for the layered sediments in this region.     

Geometrical structure of shear wave surfaces near singularity directions in anisotropic media

There are three types of surfaces which are used for studying wave propagation in anisotropic media: normal surfaces, slowness surfaces and wave surfaces. Normal surfaces and slowness surfaces have been researched in detail. Wave surfaces are the most complicated and comparatively poorly known compared with the other two. Areas of complicated geometrical structure of the wave surfaces are located in the vicinity of conical acoustic axes. There is an elliptical hole on the quick shear wave surface and complicated folds and cusps on the slow shear wave surface. Decomposition of the slow shear wave surface into smooth sheets is used for the study of its geometrical structure. Complexity of shear wave surfaces can be expressed by the number of waves corresponding to a fixed ray. An original approach to the calculation of wave normals depending on ray direction is presented.