A comparison between mesospheric wind measurements made near Christchurch (44°S, 173°E) using the high resolution doppler imager (HRDI) and a medium frequency (MF) radar

We report on the comparison of winds measured by a medium frequency (MF) radar near Christchurch, New Zealand, and by the high resolution doppler imager (HRDI). Previous comparisons have demonstrated that there can be significant differences in the winds obtained by the two techniques, and our results are no different. However, these data show relatively good agreement in the meridional direction, but large differences in the zonal direction, where the radar is regularly measuring the zonal wind as too easterly. To do the comparison, overpasses from the satellite must be obtained when it is close to the radar site. The radar data are averaged in time around the overpass because we know the radars sample phenomena which have spatial and temporal scales which make them invisible to HRDI. There are a limited number of overpass comparisons which limit our confidence in these results, but a detailed analysis of these data show that the proximity of the overpass is often an important factor in the differences obtained. Other factors examined include the influence of the local time of the overpass, and the amount of radar data averaged around the overpass time.     

Letter to the editor A comparison of F-region ion velocity observations from the EISCAT Svalbard and VHF radars with irregularity drift velocity measurements from the CUTLASS Finland HF radar

During August 1998, the UK EISCAT special programme SP-UK-CSUB, which combines operation of both the mainland VHF and Svalbard UHF incoherent scatter radars, was run for several hours around magnetic midnight on four consecutive days. The CUTLASS Finland HF coherent scatter radar was, at these times, operating in a discretionary mode, sounding on all 16 beams, one at high-time resolution. This study presents a comparison of the velocities measured by coherent and incoherent techniques during the SP-UK-CSUB experiments. Agreement, particularly between the ion velocities measured by the EISCAT Svalbard radar and irregularity drift measurements by the Finland radar, is remarkable, thereby validating the scientific integrity of both data sets. This work highlights the substantive contribution to our understanding of the solar-terrestrial environment which can be made by use in concert of incoherent and HF coherent scatter radars.     

Space Plasma Exploration by Active Radar (SPEAR): an overview of a future radar facility

SPEAR is a new polar cap HF radar facility which is to be deployed on Svalbard. The principal capabilities of SPEAR will include the generation of artificial plasma irregularities, operation as an all-sky HF radar, the excitation of ULF waves, and remote sounding of the magnetosphere. Operation of SPEAR in conjunction with the multitude of other instruments on Svalbard, including the EISCAT Svalbard radar, and the overlap of its extensive field-of-view with that of several of the HF radars in the SuperDARN network, will enable in-depth diagnosis of many geophysical and plasma phenomena associated with the cusp region and the substorm expansion phase. Moreover, its ability to produce artificial radar aurora will provide a means for the other instruments to undertake polar cap plasma physics experiments in a controlled manner. Another potential use of the facility is in field-line tagging experiments, for coordinated ground-satellite experiments. Here the scientific objectives of SPEAR are detailed, along with the proposed technical specifications of the system.     

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.     

Extrapolating EISCAT Pedersen conductances to other parts of the sky using ground-based TV auroral images

Ionospheric conductivity is not very easily measured directly. Incoherent scatter radars perhaps offer the best method but can only measure at one point in the sky at any one time and are limited in their time resolution. Statistical models of average conductivity are available but these may not be applied to individual case studies such as substorms. There are many instances where a real-time estimate of ionospheric conductivity over a large field-of-view is highly desirable at a high temporal and spatial resolution. We show that it is possible to make a reasonable estimate of the noctural height-integrated Pedersen conductivity, or conductance, with a single all-sky TV camera operating at 557.7 nm. This is not so in the case of the Hall conductance where at least two auroral wavelengths should be imaged in order to estimate additionally the energy of the precipitating particles.     

Ground clutter cancellation in incoherent radars: solutions for EISCAT Svalbard radar

Incoherent scatter radars measure ionosphere parameters using modified Thomson scatter from free electrons in the target (see e.g. Hagfors, 1997). The integrated cross section of the ionospheric scatterers is extremely small and the measurements can easily be disturbed by signals returned by unwanted targets. Ground clutter signals, entering via the antenna side lobes, can render measurements at the nearest target ranges totally impossible. The EISCAT Svalbard Radar (ESR), which started measurements in 1996, suffers from severe ground clutter and the ionosphere cannot be measured in any simple manner at ranges less than about 120–150 km, depending on the modulation employed. If the target and clutter signals have different, and clearly identifiable, properties then, in principle, there are always ways to eliminate the clutter. In incoherent scatter measurements, differences in the coherence times of the wanted and unwanted signals can be used for clutter cancellation. The clutter cancellation must be applied to all modulations, usually alternating codes in modern experiments, used for shorter ranges. Excellent results have been obtained at the ESR using a simple pulse-to-pulse clutter subtraction method, but there are also other possibilities.     

Imaging of structures in the high-latitude ionosphere: model comparisons

The tomographic reconstruction technique generates a two-dimensional latitude versus height electron density distribution from sets of slant total electron content measurements (TEC) along ray paths between beacon satellites and ground-based radio receivers. In this note, the technique is applied to TEC values obtained from data simulated by the Sheffield/UCL/SEL Coupled Thermosphere/Ionosphere/Model (CTIM). A comparison of the resulting reconstructed image with the input modelled data allows for verification of the reconstruction technique. All the features of the high-latitude ionosphere in the model data are reproduced well in the tomographic image. Reconstructed vertical TEC values follow closely the modelled values, with the F-layer maximum density (NmF2) agreeing generally within about 10%. The method has also been able successfully to reproduce underlying auroral-E ionisation over a restricted latitudinal range in part of the image. The height of the F2 peak is generally in agreement to within about the vertical image resolution (25 km).     

Acoustic tomography in the atmospheric surface layer

Acoustic tomography is presented as a technique for remote monitoring of meteorological quantities. This method and a special algorithm of analysis can directly produce area-averaged values of meteorological parameters. As a result consistent data will be obtained for validation of numerical atmospheric micro-scale models. Such a measuring system can complement conventional point measurements over different surfaces. The procedure of acoustic tomography uses the horizontal propagation of sound waves in the atmospheric surface layer. Therefore, to provide a general overview of sound propagation under various atmospheric conditions a two-dimensional ray-tracing model according to a modified version of Snell’s law is used. The state of the crossed atmosphere can be estimated from measurements of acoustic travel time between sources and receivers at different points. Derivation of area-averaged values of the sound speed and furthermore of air temperature results from the inversion of travel time values for all acoustic paths. Thereby, the applied straight ray two-dimensional tomographic model using SIRT (simultaneous iterative reconstruction technique) is characterised as a method with small computational requirements, satisfactory convergence and stability properties as well as simple handling, especially, during online evaluation.     

A new technique for evaluating mesospheric momentum balance utilizing radars and satellite data

A new method for evaluating momentum balance in the mesosphere using radar and satellite data is presented. This method is applied to radar wind data from two medium frequency installations (near Adelaide, Australia and Christchurch, New Zealand) and satellite temperature data from the Improved Stratospheric and Mesospheric Sounder (ISAMS). Because of limitations in data availability and vertical extent, the technique can only be applied to evaluate the momentum balance at 80 km above the radar sites for May 1992. The technique allows the calculation of the residual terms in the momentum balance which are usually attributed to the effects of breaking gravity waves. Although the results are inconclusive above Adelaide, this method produces values of zonal and meridional residual accelerations above Christchurch which are consistent with expectation. In both locations it is apparent that geostrophic balance is a poor approximation of reality. (This result is not dependent on a mismatch between the radar and satellite derived winds, but rather is inherent in the satellite data alone.) Despite significant caveats about data quality the technique appears robust and could be of use with data from future instruments.     

The magnetic field experiment onboard Equator-S and its scientific possibilities

The special feature of the ringcore fluxgate magnetometer on Equator-S is the high time and field resolution. The scientific aim of the experiment is the investigation of waves in the 10–100 picotesla range with a time resolution up to 64 Hz. The instrument characteristics and the influence of the spacecraft on the magnetic field measurement will be discussed. The work shows that the applied pre- and inflight calibration techniques are sufficient to suppress spacecraft interferences. The offset in spin axis direction was determined for the first time with an independent field measurement by the Equator-S Electron Drift Instrument. The data presented gives an impression of the accuracy of the measurement.     

MIDAS-W: a workstation-based incoherent scatter radar data acquisition system

The Millstone Hill Incoherent Scatter Data Acquisition System (MIDAS) is based on an abstract model of an incoherent scatter radar. This model is implemented in a hierarchical software system, which serves to isolate hardware and low-level software implementation details from higher levels of the system. Inherent in this is the idea that implementation details can easily be changed in response to technological advances. MIDAS is an evolutionary system, and the MIDAS hardware has, in fact, evolved while the basic software model has remained unchanged. From the earliest days of MIDAS, it was realized that some functions implemented in specialized hardware might eventually be implemented by software in a general-purpose computer. MIDAS-W is the realization of this concept. The core component of MIDAS-W is a Sun Microsystems UltraSparc 10 workstation equipped with an Ultrarad 1280 PCI bus analog to digital (A/D) converter board. In the current implementation, a 2.25 MHz intermediate frequency (IF) is bandpass sampled at 1 s intervals and these samples are multicast over a high-speed Ethernet which serves as a raw data bus. A second workstation receives the samples, converts them to filtered, decimated, complex baseband samples and computes the lag-profile matrix of the decimated samples. Overall performance is approximately ten times better than the previous MIDAS system, which utilizes a custom digital filtering module and array processor based correlator. A major advantage of MIDAS-W is its flexibility. A portable, single-workstation data acquisition system can be implemented by moving the software receiver and correlator programs to the workstation with the A/D converter. When the data samples are multicast, additional data processing systems, for example for raw data recording, can be implemented simply by adding another workstation with suitable software to the high-speed network. Testing of new data processing software is also greatly simplified, because a workstation with the new software can be added to the network without impacting the production system. MIDAS-W has been operated in parallel with the existing MIDAS-1 system to verify that incoherent scatter measurements by the two systems agree. MIDAS-W has also been used in a high-bandwidth mode to collect data on the November, 1999, Leonid meteor shower.     

EDI electron time-of-flight measurements on Equator-S

We present the first electron time-of-flight measurements obtained with the Electron Drift Instrument (EDI) on Equator-S. These measurements are made possible by amplitude-modulation and coding of the emitted electron beams and correlation with the signal from the returning electrons. The purpose of the time-of-flight measurements is twofold. First, they provide the drift velocity, and thus the electric field, when the distance the electrons drift in a gyro period becomes sufficiently large. Second, they provide the gyro time of the electrons emitted by the instrument, and thus the magnitude of the ambient magnetic field, allowing in-flight calibration of the flux-gate magnetometer with high precision. Results of both applications are discussed.     

ST radar observations of atmospheric waves over mountainous areas: a review

Lee and mountain waves are dominant dynamic processes in the atmosphere above mountain areas. ST VHF radars had been intensively used to investigate these wave processes. These studies are summarized in this work. After discussing features of long-period quasi-stationary lee waves, attention is drawn to the frequent occurrence of freely propagating waves of shorter periods, which seem to be more common and characteristic for wave processes generated over mountainous areas. Characteristics of these waves such as their relation to the topography and background winds, the possibility of trapping by and breaking in the tropopause region and their propagation into the stratosphere is investigated. These orographically produced waves transport energy and momentum into the troposphere and stratosphere, which is considered an important contribution to the kinetic energy of the lower atmosphere. The occurrence of inertia-gravity waves in the stratosphere had been confused with lee waves, which is discussed in conclusion. Finally further questions on mountain and lee waves are drawn up, which remain to be solved and where investigations with ST radars could play a fundamental role.     

Analysis of thick, non-planar boundaries using the discontinuity analyser

The advent of missions comprised of phased arrays of spacecraft, with separation distances ranging down to at least mesoscales, provides the scientific community with an opportunity to accurately analyse the spatial and temporal dependencies of structures in space plasmas. Exploitation of the multi-point data sets, giving vastly more information than in previous missions, thereby allows unique study of their small-scale physics. It remains an outstanding problem, however, to understand in what way comparative information across spacecraft is best built into any analysis of the combined data. Different investigations appear to demand different methods of data co-ordination. Of the various multi-spacecraft data analysis techniques developed to affect this exploitation, the discontinuity analyser has been designed to investigate the macroscopic properties (topology and motion) of boundaries, revealed by multi-spacecraft magnetometer data, where the possibility of at least mesoscale structure is considered. It has been found that the analysis of planar structures is more straightforward than the analysis of non-planar boundaries, where the effects of topology and motion become interwoven in the data, and we argue here that it becomes necessary to customise the analysis for non-planar events to the type of structure at hand. One issue central to the discontinuity analyser, for instance, is the calculation of normal vectors to the structure. In the case of planar and ‘thin’ non-planar structures, the method of normal determination is well-defined, although subject to uncertainties arising from unwanted signatures. In the case of ‘thick’, non-planar structures, however, the method of determination becomes particularly sensitive to the type of physical sampling that is present. It is the purpose of this article to firstly review the discontinuity analyser technique and secondly, to discuss the analysis of the normals to thick non-planar structures detected in magnetometer data.     

Validation of the CUTLASS HF radar gravity wave observing capability using EISCAT CP-1 data

Quasi-periodic fluctuations in the returned ground-scatter power from the SuperDARN HF radars have been linked to the passage of medium-scale gravity waves. We have applied a technique that extracts the first radar range returns from the F-region to study the spatial extent and characteristics of these waves in the CUTLASS field-of-view. Some ray tracing was carried out to test the applicability of this method. The EISCAT radar facility at Tromsø is well within the CUTLASS field-of-view for these waves and provides a unique opportunity to assess independently the ability of the HF radars to derive gravity wave information. Results from 1st March, 1995, where the EISCAT UHF radar was operating in its CP-1 mode, demonstrate that the radars were in good agreement, especially if one selects the electron density variations measured by EISCAT at around 235 km. CUTLASS and EISCAT gravity wave observations complement each other; the former extends the spatial field of view considerably, whilst the latter provides detailed vertical information about a range of ionospheric parameters.     

4D tropospheric tomography using GPS slant wet delays

Tomographic techniques are successfully applied to obtain 4D images of the tropospheric refractivity in a local dense network of global positioning system (GPS) receivers. We show here how GPS data are processed to obtain the tropospheric slant wet delays and discuss the validity of the processing. These slant wet delays are the observables in the tomographic processing. We then discuss the inverse problem in 4D tropospheric tomography making extensive use of simulations to test the system and define the resolution and the impact of noise. Finally, we use data from the Kilauea network in Hawaii for February 1, 1997, and a local 4 × 4 × 40 voxel grid on a region of 400 km² and 15 km in height to produce the corresponding 4D wet refractivity fields, which are then validated using forecast analysis from the European Center for Medium Range Weather Forecast (ECMWF). We conclude that tomographic techniques can be used to monitor the troposphere in time and space.     

Inversion, error analysis, and validation of GPS/MET occultation data

The global positioning system meteorology (GPS/MET) experiment was the first practical demonstration of global navigation satellite system (GNSS)-based active limb sounding employing the radio occultation technique. This method measures, as principal observable and with millimetric accuracy, the excess phase path (relative to propagation in vacuum) of GNSS-transmitted radio waves caused by refraction during passage through the Earth’s neutral atmosphere and ionosphere in limb geometry. It shows great potential utility for weather and climate system studies in providing an unique combination of global coverage, high vertical resolution and accuracy, long-term stability, and all-weather capability. We first describe our GPS/MET data processing scheme from excess phases via bending angles to the neutral atmospheric parameters refractivity, density, pressure and temperature. Special emphasis is given to ionospheric correction methodology and the inversion of bending angles to refractivities, where we introduce a matrix inversion technique (instead of the usual integral inversion). The matrix technique is shown to lead to identical results as integral inversion but is more directly extendable to inversion by optimal estimation. The quality of GPS/MET-derived profiles is analyzed with an error estimation analysis employing a Monte Carlo technique. We consider statistical errors together with systematic errors due to upper-boundary initialization of the retrieval by a priori bending angles. Perfect initialization and properly smoothed statistical errors allow for better than 1 K temperature retrieval accuracy up to the stratopause. No initialization and statistical errors yield better than 1 K accuracy up to 30 km but less than 3 K accuracy above 40 km. Given imperfect initialization, biases ≫ 2 K propagate down to below 30 km height in unfavorable realistic cases. Furthermore, results of a statistical validation of GPS/MET profiles through comparison with atmospheric analyses of the European Centre for Medium-range Weather Forecasts (ECMWF) are presented. The comparisons indicate the high utility of the occultation data in that very good agreement of upper troposphere/lower stratosphere temperature (better than 1.5 K rms, ≪ 0.5 K bias) is found for a region (Europe+USA) where the ECMWF analyses are known to be good, but poorer agreement for a region (Southern Pacific) where the analyses are known to be degraded.     

21世纪重力与磁法勘探的展望

对21世纪重力与磁法勘探的仪器，数据处理技术，解释理论与方法，应用领域等方面的发展方向进行了分析与展望，发展航空标量，矢量，梯度重力测量和航空全梯度磁力测量，三分量磁力测量，提高综合信息采集能力；开展卫星重磁测量，综合卫星，航空，地面重磁测量资料研究地球结构与构造；发展高精度数据处理技术；重磁异常弱信号的提取，不同深度重磁异常的划分，低纬底变倾角化磁极以及位场面延拓；发展复杂条件下三维重磁场多参数综合反演可视化技术以及快速自动反演技术；探索磁性多参数的应用新领域，充分发挥磁法在环境污染调查中的作用并开拓应用新领域。     

Plasma density over Svalbard during the ISBJØRN campaign

In 1997, reliable operation of the EISCAT Svalbard Radar (ESR) was achieved and a rocket launching facility at Ny Ålesund on Svalbard (79°N, 12°E) (SVALRAK) was established. On 20 November, 1977, the first instrumented payload was launched from SVALRAK. Although the payload configuration had been flown many times previously from Andøya Rocket Range on the Norwegian mainland, this presented an unprecedented in situ determination of positive ion density over Svalbard. Simultaneously, ESR measured similar density profiles but in a higher altitude regime. We have combined the ESR measurements with ionosonde data to establish a calibration and subsequently combined the ground-based and in situ determined profiles to give a composite positive ion density profile from the mesosphere to the thermosphere.     

Investigation of the relationship between optical auroral forms and HF radar E region backscatter

The SuperDARN HF radars have been employed in the past to investigate the spectral characteristics of coherent backscatter from L-shell aligned features in the auroral E region. The present study employs all-sky camera observations of the aurora from Husafell, Iceland, and the two SuperDARN radars located on Iceland, Þykkvibær and Stokkseyri, to determine the optical signature of such backscatter features. It is shown that, especially during quiet geomagnetic conditions, the backscatter region is closely associated with east-west aligned diffuse auroral features, and that the two move in tandem with each other. This association between optical and radar aurora has repercussions for the instability mechanisms responsible for generating the E region irregularities from which radars scatter. This is discussed and compared with previous studies investigating the relationship between optical and VHF radar aurora. In addition, although it is known that E region backscatter is commonly observed by SuperDARN radars, the present study demonstrates for the first time that multiple radars can observe the same feature to extend over at least 3 h of magnetic local time, allowing precipitation features to be mapped over large portions of the auroral zone.