New techniques for auroral irregularity studies with COSCAT

The COSCAT system enables the detection of E-region auroral backscatter with the EISCAT remote receivers at magnetic aspect angles close to 90°. This is achieved by utilising a low-power transmitter stationed in Oulu, Finland. Many important observations of E-region irregularities have been achieved with this simple experiment. Recent studies have attempted to push the COSCAT system to its experimental limits. Firstly, the CW signal has been phase-modulated with 13-bit Barker codes with baud lengths of 40, 70 and 100 s. Interpretation of the received power allows the spatial distribution of the auroral scatterers to be determined. The second advance is in the use of a sophisticated correlator program which allows data to be buffered within the correlator at very high time resolution. This enables the coherent backscatter power to be sampled every 12.5 ms and the full auto-correlation function to be measured every 100 ms. These measurements allow the COSCAT system to be employed for the first time in an investigation of the growth and decay of the auroral irregularities.     

Multi-instrument observations of the electric and magnetic field structure of omega bands

High time resolution data from the CUTLASS Finland radar during the interval 01:30–03:30 UT on 11 May, 1998, are employed to characterise the ionospheric electric field due to a series of omega bands extending 5° in latitude at a resolution of 45 km in the meridional direction and 50 km in the azimuthal direction. E-region observations from the STARE Norway VHF radar operating at a resolution of 15 km over a comparable region are also incorporated. These data are combined with ground magnetometer observations from several stations. This allows the study of the ionospheric equivalent current signatures and height integrated ionospheric conductances associated with omega bands as they propagate through the field-of-view of the CUTLASS and STARE radars. The high-time resolution and multi-point nature of the observations leads to a refinement of the previous models of omega band structure. The omega bands observed during this interval have scale sizes 500 km and an eastward propagation velocity 0.75 km s^–1. They occur in the morning sector (05 MLT), simultaneously with the onset/intensification of a substorm to the west during the recovery phase of a previous substorm in the Scandinavian sector. A possible mechanism for omega band formation and their relationship to the substorm phase is discussed.     

Non-magnetic aspect sensitive auroral echoes from the lower E region observed at 50 MHz

Backscatter from E-region irregularities was observed at aspect angles close to 90° (almost parallel to the direction of the magnetic field) using the ALOMAR SOUSY radar at Andoya/Norway. Strong electric fields and increased E-region electron temperatures simultaneously measured with the incoherent scatter facility EISCAT proved that the Farley-Buneman plasma instability was excited. In addition, strong particle precipitation was present as inferred from EISCAT electron densities indicating that the gradient drift instability may have been active, too. Backscatter at such large aspect angles was not expected and has not been observed before. The characteristics of the observed echoes, however, are in many aspects completely different from usual auroral radar results: the Doppler velocities are only of the order of 10 m/s, the half-width of the spectra is around 5 m/s, the echoes originate at altitudes well below 100 km, and they seem to be not aspect-sensitive with respect to the magnetic field direction. We, therefore, conclude that the corresponding irregularities are not caused by the mentioned instabilities and that other mechanism have to be invoked.     

Simultaneous observations at different altitudes of ionospheric backscatter in the eastward electrojet

A common feature of evening near-range ionospheric backscatter in the CUTLASS Iceland radar field of view is two parallel, approximately L-shell-aligned regions of westward flow which are attributed to irregularities in the auroral eastward electrojet region of the ionosphere. These backscatter channels are separated by approximately 100–200 km in range. The orientation of the CUTLASS Iceland radar beams and the zonally aligned nature of the flow allows an approximate determination of flow angle to be made without the necessity of bistatic measurements. The two flow channels have different azimuthal variations in flow velocity and spectral width. The nearer of the two regions has two distinct spectral signatures. The eastern beams detect spectra with velocities which saturate at or near the ion-acoustic speed, and have low spectral widths (less than 100ms^–1), while the western beams detect lower velocities and higher spectral widths (above 200ms^–1). The more distant of the two channels has only one spectral signature with velocities above the ionacoustic speed and high spectral widths. The spectral characteristics of the backscatter are consistent with E-region scatter in the nearer channel and upper-E-region or F-region scatter in the further channel. Temporal variations in the characteristics of both channels support current theories of E-region turbulent heating and previous observations of velocity-dependent backscatter cross-section. In future, observations of this nature will provide a powerful tool for the investigation of simultaneous E- and F-region irregularity generation under similar (nearly co-located or magnetically conjugate) electric field conditions.     

基于三亚VHF雷达的场向不规则体观测研究:1.电离层E区连续性回波

基于三亚(109.6°E,18.4°N)VHF电离层相干散射雷达观测,分析了我国低纬电离层E区场向不规则结构连续性回波的发生特征.研究结果表明:白天,E区连续性回波的多普勒速度范围为-50至25m/s,多普勒宽度主要分布在20至70m/s;连续性回波的高度大约以1km/h的速度缓慢下降,与偶发E层(Es)底部所在高度(hbEs)有很好的相关性,表明在背景电场影响下,Es经梯度漂移不稳定性产生场向不规则结构,引起E区连续性回波.夜间,E区连续性回波的多普勒速度范围为-50至50m/s,多普勒宽度为20至110m/s,回波在时间-高度-强度图上常呈现多层结构,可能与潮汐引起的多个离子层相关;而E区连续性回波的短暂中断,以及120km以上高E区连续性回波的发生,则可能归因于赤道扩展F极化电场的影响.此外,对E区连续性回波多普勒速度与全天空流星雷达风场观测的比较发现,在100km以下,多普勒速度与子午风场有很好的相关性.     

A statistical study of underestimates of wind speeds by VHF radar

Comparisons are made between horizontal wind measurements carried out using a VHF-radar system at Aberystwyth (52.4°N, 4.1°W) and radiosondes launched from Aberporth, some 50 km to the southwest. The radar wind results are derived from Doppler wind measurements at zenith angles of 6° in two orthogonal planes and in the vertical direction. Measurements on a total of 398 days over a 2-year period are considered, but the major part of the study involves a statistical analysis of data collected during 75 radiosonde flights selected to minimise the spatial separation of the two sets of measurements. Whereas good agreement is found between the two sets of wind direction, radar-derived wind speeds show underestimates of 4–6% compared with radiosonde values over the height range 4–14 km. Studies of the characteristics of this discrepancy in wind speeds have concentrated on its directional dependence, the effects of the spatial separation of the two sets of measurements, and the influence of any uncertainty in the radar measurements of vertical velocities. The aspect sensitivity of radar echoes has previously been suggested as a cause of underestimates of wind speeds by VHF radar. The present statistical treatment and case-studies show that an appropriate correction can be applied using estimates of the effective radar beam angle derived from a comparison of echo powers at zenith angles of 4.2° and 8.5°.     

CUTLASS/IMAGE observations of high-latitude convection features during substorms

The CUTLASS Finland HF radar has been operational since February 1995. The radar frequently observes backscatter during the midnight sector from a latitude range 70–75° geographic, latitudes often associated with the polar cap. These intervals of backscatter occur during intervals of substorm activity, predominantly in periods of relatively quiet magnetospheric activity, with Kp during the interval under study being 2-and K_P for the day being only 8-. During August 1995 the radar ran in a high time resolution mode, allowing measurements of line-of-sight convection velocities along a single beam with a temporal resolution of 14s, and measurement of a full spatial scan of line-of-sight convection velocities every four minutes. Data from such scans reveal the radar to be measuring return flow convection during the interval of substorm activity. For three intervals during the period under study, a reduction in the spatial extent of radar backscatter occurred. This is a consequence of D region HF absorption and its limited extent in the present study is probably a consequence of the high latitude of the substorm activity, with the electrojet centre lying between 67° and 71° geomagnetic latitude. The high time resolution beam of the radar additionally demonstrates that the convection is highly time dependent. Pulses of equatorward flow exceeding 600 m s^–1 are observed with a duration of 5 min and a repetition period of 8 min. Their spatial extent in the CUTLASS field of view was 400–500 km in longitude, and 300–400 km in latitude. Each pulse of enhanced equatorward flow was preceded by an interval of suppressed flow and enhanced ionospheric Hall conductance. The transient features are interpreted as being due to ionospheric current vortices associated with field aligned current pairs. The relationship between these observations and substorm phenomena in the magnetotail is discussed.     

P-wave velocity anomalies in the earth's mantle from the Uppsala array observations

Summary Distribution of compressional-wave velocities in the mantle is determined fromdT/d measurements using the Uppsala seismograph array station (UPSAS). Short-period vertical-component seismograms from 181 events in the epicentral distance range 16°–100° have been used. The velocity distribution shows anomalous variations at depths of 750, 1500, 1800, 2300 and 2550 km. Evidence of lateral heterogeneity beneath the northern part of the Asian continent, in the depth range 1700–2300 km, is discussed. Computed travel times, based on this velocity-depth relation, are tested by an examination of travel-time residuals from the Long Shot and Milrow explosions on Amchitka, Aleutian Islands.     

Crustal and upper mantle structures of the brazilian coast

The Rayleigh wave phase and group velocities in the period range of 24–39 sec, obtained from two earthquakes which occurred in northeastern brazil and which were recorded by the Brazilian seismological station RDJ (Rio de Janeiro), have been used to study crustal and upper mantle structures of the Brazilian coastal region. Three crustal and upper mantle models have been tried out to explain crustal and upper mantle structures of the region. The upper crust has not been resolved, due basically to the narrow period range of the phase and group velocities data. The phase velocity inversions have exhibited good resolutions for both lower crust and upper mantle, with shear wave velocities characteristic of these regions. The group velocity data inversions for these models have showed good results only for the lower crust. The shear wave velocities of the lower crust (3.86 and 3.89 km/sec), obtained with phase velocity inversions, are similar to that (=3.89 km/sec) found byHwang (1985) to the eastern South American region, while group velocity inversions have presented shear velocity (=3.75 km/sec) similar to that (=3.78 km/sec) found byLazcano (1972) to the Brazilian shield. It was not possible to define sharply the crust-mantle transition, but an analysis of the phase and group velocity inversions results has indicated that the total thickness of the crust should be between 30 and 39 km. The crustal and upper mantle model, obtained with phase velocity inversion, can be used as a preliminary model for the Brazilian coast.     

The generation of non aspect sensitive plasma density irregularities by field aligned drifts in the lower ionosphere

A theory of the generation of plasma density irregularities with virtually no aspect sensitivity, in the lower ionosphere at high latitudes, by electron drifts aligned with the geomagnetic field, is presented. The theory is developed through fluid equations in which the destabilising mechanism involves positive feedback from electron collisional heating. When field aligned electron drift speeds exceed a few km s^–1, this effect destabilises waves with wavelengths in excess of a few tens of metres in the lower E-region, where collisional effects are sufficiently large. Furthermore, the threshold conditions are almost independent of the wave propagation direction and the unstable waves propagate at speeds well below the ion acoustic speed. The role that this new instability may play in recent radar backscatter observations of short scale irregularities propagating in directions close to that of the geomagnetic field, in the lower E-region is also considered.     

Structure of the crust in the Black Sea and adjoining regions from surface wave data

Group velocities of Rayleigh and Love waves along the paths across the Black Sea and partly Asia Minor and the Balkan Peninsula are used to estimate lateral variations of the crustal structure in the region. As a first step, lateral variations of group velocities for periods in the range 10–20 s are determined using a 2D tomography method. Since the paths are oriented predominantly in NE–SW or N–S direction, the resolution is estimated as a function of azimuth. The local dispersion curves are actually averaged over the extended areas stretched in the predominant direction of the paths. The size of the averaging area in the direction of the best resolution is approximately 200 km. As a second step, the local averaged dispersion curves are inverted to vertical sections of S-wave velocities. Since the dispersion curves in the 10–20 s period range are mostly affected by the upper crustal structure, the velocities are estimated to a depth of approximately 25 km. Velocity sections along 43° N latitude are determined separately from Rayleigh and Love wave data. It is shown that the crust under the sea contains a low-velocity sedimentary layer of 2–3 km thickness, localized in the eastern and western deeps, as found earlier from DSS data. Beneath the sedimentary layer, two layers are present with velocity values lying between those of granite and consolidated sediments. Velocities in these layers are slightly lower in the deeps, and the boundaries of the layers are lowered. S-wave velocities obtained from Love wave data are found to be larger than those from Rayleigh wave data, the difference being most pronounced in the basaltic layer. If this difference is attributed to anisotropy, the anisotropy coefficient = (SH - SV)/Smean is reasonable (2–3%) in the upper layers, and exceeds 9% in the basaltic layer.     

Upper mantle structure of the Apulian plate from Rayleigh waves

Phase velocities of Rayleigh waves for the Adriatic Sea area are obtained in the period range 25–190 sec along the path (l'Aquila-Trieste) AQU-TRI and 20–167 sec along the path (Trieste-Bari) TRI-BAI.The phase velocities are systematically higher than the known values for the surrounding regions. The data inversion indicates the presence of a lithosphere typical of stable continental areas with clear high-velocity lid (V _s 4.6 km/sec) overlying a well developed low velocity zone (V _s 4.2 km/sec).P. F. Geodinamica C.N.R., Roma Pubbl. N. 189.     

Sequential sporadic-E layers at low latitudes in the Indian sector

A study of the formation and movement of sequential Sporadic-E layers observed during the night-time hours at two Indian low-latitude stations, SHAR(dip 10°N) and Waltair (dip 20°N) shows that the layer are formed around 19:00 h. IST at altitudes of ≈180 km. They descend to the normal E-region altitude of about 100 km in three to four hours and becomes blanketing type of Es before they disappear. However, the absence of these descending layers at an equatorial station, Trivandrum (dip 2°N) gives the experimental evidence for wind shear theory. The meridional neutral wind derived from the height variation of the F-layer showed significant poleward wind during the descent of these layers. Hence it is inferred that these layers are formed as a consequence of the convergence of plasma by the poleward wind and the equatorward propagating gravity waves (inferred from the height fluctuations of F-layer).     

Normal moveout velocity for pure‐mode and converted waves in layered orthorhombic medium

We study the azimuthally dependent hyperbolic moveout approximation for small angles (or offsets) for quasi‐compressional, quasi‐shear, and converted waves in one‐dimensional multi‐layer orthorhombic media. The vertical orthorhombic axis is the same for all layers, but the azimuthal orientation of the horizontal orthorhombic axes at each layer may be different. By starting with the known equation for normal moveout velocity with respect to the surface‐offset azimuth and applying our derived relationship between the surface‐offset azimuth and phase‐velocity azimuth, we obtain the normal moveout velocity versus the phase‐velocity azimuth. As the surface offset/azimuth moveout dependence is required for analysing azimuthally dependent moveout parameters directly from time‐domain rich azimuth gathers, our phase angle/azimuth formulas are required for analysing azimuthally dependent residual moveout along the migrated local‐angle‐domain common image gathers. The angle and azimuth parameters of the local‐angle‐domain gathers represent the opening angle between the incidence and reflection slowness vectors and the azimuth of the phase velocity ψ_phs at the image points in the specular direction. Our derivation of the effective velocity parameters for a multi‐layer structure is based on the fact that, for a one‐dimensional model assumption, the horizontal slowness and the azimuth of the phase velocity ψ_phs remain constant along the entire ray (wave) path. We introduce a special set of auxiliary parameters that allow us to establish equivalent effective model parameters in a simple summation manner. We then transform this set of parameters into three widely used effective parameters: fast and slow normal moveout velocities and azimuth of the slow one. For completeness, we show that these three effective normal moveout velocity parameters can be equivalently obtained in both surface‐offset azimuth and phase‐velocity azimuth domains.     

2-D Crustal Velocity Structure Along Hirapur-Mandla Profile in Central India: An Update

The 2-D crustal velocity model along the Hirapur-Mandla DSS profile across the Narmada-Son lineament in central India (Murty et al., 1998) has been updated based on the analysis of some short and discontinuous seismic wide-angle reflection phases. Three layers, with seismic velocities of 6.5–6.7, 6.35–6.40 and 6.8 km s^–1, and upper boundaries located approximately at 8, 17 and 22 km depth respectively, have been identified between the basement (velocity 5.9 km s^–1) and the uppermost mantle (velocity 7.8 km s^–1). The layer with 6.5–6.7 km s^–1 velocity is thin and is less than 2-km deep between the Narmada north (at Katangi) and south (at Jabalpur) faults. The upper crust shows a horst feature between these faults, which indicates that the Narmada zone acts as a ridge between two pockets of mafic intrusion in the upper crust. The Moho boundary, at 40–44 km depth and the intra-crustal layers exhibit an upwarp suggesting that the Narmada faults have deep origins, involving deep-seated tectonics. A smaller intrusive thickness between the Narmada faults, as compared to those beyond these faults, suggests that the intrusive activities on the two sides are independent. This further suggests that the two Narmada faults may have been active at different geological times. The seismic model is constrained by 2-D gravity modeling. The gravity highs on either side of the Narmada zone are due to the effect of the high velocity/high density mafic intrusion at upper crustal level.     

Rupture process of the 19 August 1992 Susamyr, Kyrgyzstan, earthquake

The Susamyr earthquake of August 19, 1992 in Kyrgyzstan is one of the largest events (M_s = 7.4, M_b = 6.8) of this century in this region of Central Asia. We used broadband and long period digital data from IRIS and GEOSCOPE networks to investigate the source parameters, and their space-time distribution by modeling both body and surface waves. The seismic moment (M₀ = 6.8 × 10¹⁹ N m) and the focal mechanism were determined from frequency-time analysis (FTAN) of the fundamental mode of long period surface waves (100–250 s). Then, the second order integral moments of the moment-rate release were estimated from the amplitude spectra of intermediate period surface waves(40–70 s). From these moments we determined a source duration of 11–13 s, major and minor axes of the source of 30 km and 10–22 km, respectively; and an instant centroid velocity of 1.2 km/s. Finally, we performed a waveform inversion of P and SH waves at periods from 5–60 s. We found a source duration of 18–20 s, longer than the integral estimate from surface wave amplitudes. All the other focal parameters inverted from body waves are similar to those obtained by surface waves ( = 87° ± 6°, = 49° ± 6°, = 105° ± 3°, h = 14 ± 2 km, and M₀ = 5.8 ± 0.7 × 10¹⁹ N m). The initial rupture of this shallow earthquake was located at the south-west border of Susamyr depression in the western part of northern Tien Shan. A finite source analysis along the strike suggests a westward propagation of the rupture. The main shock of this event was preceded 2 s earlier by small foreshock. The main event was almost immediately followed by a very strong series of aftershocks. Our surface and body wave inversion results agree with the general seismotectonic features of the region.     

Two-dimensional kinematic and rheological modeling of the 1912 pyroclastic flow,Katmai, Alaska

Pyroclastic flow emplacement is strongly influenced by eruption column height. A surface along which kinetic energy is zero theoretically connects the loci of eruption column collapse with all coeval ignimbrite termini. This surface is reconstructed as a two-dimensional energy line for the 1912 Katmai pyroclastic flow in the Valley of Ten Thousand Smokes from mapped flow termini and the runup of the ignimbrite onto obstructions and through passes. Extrapolation of the energy line to the vicinity of the source vent at Novarupta suggests the eruption column which generated the ignimbrite eruption was approximately 425 m high. The 1912 pyroclastic flow travelled about 25 km downvalley. Empirical velocity data calculated from runup elevations and surveyed centrifugal superelevations indicate initial velocities near Novarupta were greater than 79–88 m s^–1. The flow progressively decelerated and was travelling only 2–8 m s^–1 when it crossed a moraine 16 km downvalley. The constant slope of the energy line away from Novarupta suggests the flow was systematically slowed by internal and basal friction. Using a simple physical model to calculate flow velocities and a constant kinetic friction coefficient (Heim coefficient) of 0.04 derived from the reconstructed energy line, the flow is estimated to have decelerated at an average rate of –0.16 m s^–2 and to have taken approximately 9.5 minutes to travel 25 km down the Valley of Ten Thousand Smokes. The shear strength of the flowing ignimbrite at the moraine was approximately 0.5 kPa, and its Bingham viscosity when it crossed the moraine was 3.5 × 10³ P. If the flow was Newtonian, its viscosity was 4.2 × 10³ P. Reynolds and Froude numbers at the moraine were only 41–62 and 0.84–1.04, respectively, indicating laminar, subcritical flow.     

Letter to the editor Complete maps of the aspect sensitivity of VHF atmospheric radar echoes

Using the MU radar at Shigaraki, Japan (34.85°N, 136.10°E), we measure the power distribution pattern of VHF radar echoes from the mid-troposphere. The large number of radar beam-pointing directions (320) allows the mapping of echo power from 0° to 40° from zenith, and also the dependence on azimuth, which has not been achieved before at VHF wavelengths. The results show how vertical shear of the horizontal wind is associated with a definite skewing of the VHF echo power distribution, for beam angles as far as 30° or more from zenith, so that aspect sensitivity cannot be assumed negligible at any beam-pointing angle that most existing VHF radars are able to use. Consequently, the use of VHF echo power to calculate intensity of atmospheric turbulence, which assumes only isotropic backscatter at large beam zenith angles, will sometimes not be valid.     

Substorm associated radar auroral surges: a statistical study and possible generation model

Substorm-associated radar auroral surges (SARAS) are a short lived (15–90 minutes) and spatially localised (5° of latitude) perturbation of the plasma convection pattern observed within the auroral E-region. The understanding of such phenomena has important ramifications for the investigation of the larger scale plasma convection and ultimately the coupling of the solar wind, magnetosphere and ionosphere system. A statistical investigation is undertaken of SARAS, observed by the Sweden And Britain Radar Experiment (SABRE), in order to provide a more extensive examination of the local time occurrence and propagation characteristics of the events. The statistical analysis has determined a local time occurrence of observations between 1420 MLT and 2200 MLT with a maximum occurrence centred around 1700 MLT. The propagation velocity of the SARAS feature through the SABRE field of view was found to be predominately L-shell aligned with a velocity centred around 1750 ms^–1 and within the range 500 m s^–1 and 3500 m s^–1. This comprehensive examination of the SARAS provides the opportunity to discuss, qualitatively, a possible generation mechanism for SARAS based on a proposed model for the production of a similar phenomenon referred to as sub-auroral ion drifts (SAIDs). The results of the comparison suggests that SARAS may result from a similar geophysical mechanism to that which produces SAID events, but probably occurs at a different time in the evolution of the event.     

Effect of basal guided waves on landslides

A landslide model riding on basal guided waves is investigated to explain lower net frictions at high slide velocities from the wave-theoretical point of view. It is shown that there is a wave propagated along the basal layer at the phase velocity equal to the slide velocity, as well as a guided wave with considerably higher phase velocities propagated likewise along the basal layer as a leaking mode at low slide velocities. With increasing slide velocity the phase velocity of the guided wave decreases until it is equal to that of the slide mass. Over this threshold slide velocity, a sonic boom is generated around the basal layer, and the shock contributes to a loosening of the slide mass into a fluidized state. Landslides on long slide-ways are more liable to exceed this threshold velocity since their slide velocities tend to be higher than those on short slide-ways of a similar shape. Hence, the reduction of net friction of landslides can possibly be better correlated with the lengths of slide-ways than with the volumes of landslides as is widely maintained.