Detection of Local Inhomogeneities by Means of Cross-borehole Seismic Tomography

—This numerical study is devoted to distortions of local anomalies which are revealed by 2-D and 3-D cross-borehole seismic tomography based on first arrival travel times. The fact that prolonged subvertical anomalies may disappear in the final velocity model is well known. But distortions are also inherent to images of local inhomogeneities. These distortions are accompanied with false anomalies of the opposite sign located just above and below true ones. I offer a semi-empirical explanation of their existence, proving that they are an inherent consequence of observation geometry incompleteness. To improve the quality of the model obtained with cross-borehole tomography, a 3-D acquisition geometry may be used when additional boreholes are positioned around the target object. Results of 3-D tomographic experiments however may be considerably distorted for another reason nonaccounted bending of seismic rays.¶To generate the travel times, a ray-tracing procedure making use of the graph theory was applied. Tomographic inversion was performed by an algorithm based on the assumption of the unknown function smoothness.     

Synthetic Seismograms for a Layered Earth Geological Model Using the Absorption and Dispersion Phenomena

—An algorithm has been developed to compute the dispersive and dissipative seismic response using FUTTERMAN’S (1962) third attenuation-dispersion relationship. In the computation, frequency-dependent velocity and quality factor Q have been used but in the case of the nondispersive synthetic seismogram, frequency-independent velocity has been used. The model’s parameters are density, phase velocity, quality factors and thicknesses of the layers. Dispersive and nondispersive synthetic seismograms have been computed with and without absorption for a layered earth geological model. Fast Fourier transform (FFT) technique has been adopted for converting the frequency domain response into the time domain. The frequency spacing, Δf = 0.976?Hz, has been considered to avoid the aliasing effect. The results have revealed changes in the reflected waveforms in the frequency domain as well as in the time domain for absorption and dispersion cases. It is also concluded that dispersion reduces the arrival time and this effect is increasing with the travel time. The effect of constant Q on the seismic response has also been studied.     

Afternoon mid-latitude current system and low-latitude geomagnetic field asymmetry during geomagnetic storms

For four geomagnetic storms of middle intensity the relationship between the low-latitude magnetic field asymmetry using ASY indices and the intensity of the auroral eastward and westward electro-jet was considered. It was asked whether there exists a connection between ASY and the eastward electrojet. To answer this question equivalent current systems were estimated in mid-latitudes. It was found that the observations obviously show no correlative relationship between the low-latitude magnetic-field asymmetry and the eastward electrojet, whereas one exists between ASY and the westward electrojet. To explain the generally accepted common three-dimensional current system between the partial ring current and the eastward electrojet, a condensor model of the three-dimensional current system was developed. It could be shown that the short periodic variations of the partial ring current are shielded by the condensor and cannot influence the eastward-electrojet current.     

Influence of the Load Wavelength on the Permeability of a Viscosity Interface in the Mantle

Assuming a radially stratified Newtonian mantle in a steady-state approximation, we demonstrate that the permeability of a viscosity interface at 660-km depth strongly depends on the wavelength of buoyancy forces driving the flow. The flow induced by long-wavelength loads penetrates through the boundary freely even if the viscosity increases by two orders. In contrast, the boundary is practically impermeable for short-wavelength loads located in the upper mantle. Thus, a stepwise increase of viscosity is a significant obstacle for small descending features in the upper mantle, but huge upper mantle downwellings, or upwellings formed in the-lower mantle can overcome it easily. This indicates that certain care is necessary in interpreting the seismic structure of the mantle by means of flow models. The global tomographic image includes only the first few degrees of the harmonic series and, consequently, its interpretation in terms of a present-day flow field results in a predominantly whole-mantle circulation even for extreme viscosity contrasts.     

Statistics of the largest geomagnetic storms per solar cycle (1844–1993)

A previous application of extreme-value statistics to the first, second and third largest geomagnetic storms per solar cycle for nine solar cycles is extended to fourteen solar cycles (1844–1993). The intensity of a geomagnetic storm is measured by the magnitude of the daily aa index, rather than the half-daily aa index used previously. Values of the conventional aa index (1868– 1993), supplemented by the Helsinki Ak index (1844–1880), provide an almost continuous, and largely homogeneous, daily measure of geomagnetic activity over an interval of 150 years. As in the earlier investigation, analytic expressions giving the probabilities of the three greatest storms (extreme values) per solar cycle, as continuous functions of storm magnitude (ad), are obtained by least-squares fitting of the observations to the appropriate theoretical extreme-value probability functions. These expressions are used to obtain the statistical characteristics of the extreme values; namely, the mode, median, mean, standard deviation and relative dispersion. Since the Ak index may not provide an entirely homogeneous extension of the aa index, the statistical analysis is performed separately for twelve solar cycles (1868–1993), as well as nine solar cycles (1868–1967). The results are utilized to determine the expected ranges of the extreme values as a function of the number of solar cycles. For fourteen solar cycles, the expected ranges of the daily aa index for the first, second and third largest geomagnetic storms per solar cycle decrease monotonically in magnitude, contrary to the situation for the half-daily aa index over nine solar cycles. The observed range of the first extreme daily aa index for fourteen solar cycles is 159–352 nT and for twelve solar cycles is 215–352 nT. In a group of 100 solar cycles the expected ranges are expanded to 137–539 and 177–511 nT, which represent increases of 108% and 144% in the respective ranges. Thus there is at least a 99% probability that the daily aa index willAlso Visiting Reader in Physics, University of Sussex, Palmer, Brighton, BN1 9QH, UK     

Study of sporadic-E clouds by backscatter radar

It is shown that swept-frequency backscatter ionograms covering a range of azimuths can be used to study the dynamics of sporadic-E clouds. A simple technique based on analytic ray tracing can be used to simulate the observed narrow traces associated with E_s patches. This enables the location and extent of the sporadic-E clouds to be determined. The motion of clouds can then be determined from a time sequence of records. In order to demonstrate the method, results are presented from an initial study of 5 days of backscatter ionograms from the Jindalee Stage B data base obtained during March-April 1990. Usually 2–3 clouds were observed each day, mainly during the evening and up to midnight. The clouds lasted from 1–4 h and extended between 30°–80° in azimuth and 150–800 km in range. The clouds were mostly stationary or drifted generally westward with velocities of up to 80 ms^–1. Only one cloud was observed moving eastward.     

Artificial periodic irregularities in the auroral ionosphere

Artificial periodic irregularities (API) are produced in the ionospheric plasma by a powerful standing electromagnetic wave reflected off the F region. The resulting electron-density irregularities can scatter other high-frequency waves if the Bragg scattering condition is met. Such measurements have been performed at mid-latitudes for two decades and have been developed into a useful ionospheric diagnostic technique. We report here the first measurements from a high-latitude station, using the EISCAT heating facility near Troms0, Norway. Both F-region and lower-altitude ionospheric echoes have been obtained, but the bulk of the data has been in the E and D regions with echoes extending down to 52-km altitude. Examples of API are shown, mainly from the D region, together with simultaneous VHP incoherent-scatter-radar (ISR) data. Vertical velocities derived from the rate of phase change during the irregularity decay are shown and compared with velocities derived from the ISR. Some of the API-derived velocities in the 75–115-km height range appear consistent with vertical neutral winds as shown by their magnitudes and by evidence of gravity waves, while other data in the 50–70-km range show an unrealistically large bias. For a comparison with ISR data it has proved difficult to get good quality data sets overlapping in height and time. The initial comparisons show some agreement, but discrepancies of several metres per second do not yet allow us to conclude that the two techniques are measuring the same quantity. The irregularity decay time-constants between about 53 and 70 km are compared with the results of an advanced ion-chemistry model, and height profiles of recorded signal power are compared with model estimates in the same altitude range. The calculated amplitude shows good agreement with the data in that the maximum occurs at about the same height as that of the measured amplitude. The calculated time-constant agrees very well with the data below 60 km but is larger above 60 km by a factor of up to 2 at 64 km. The comparisons with the model are considered to be a good basis for more refined comparisons.