Radio plasma fringes as guide stars: tracking the global tilt

Erez N Ribak
Department of Physics, Technion, Haifa 32000, Israel

Roberto Ragazzoni
Observatory of Padua, viccolo dell'Osservatorio 5, 35122 Padua, Italy

Vadim A Parfenov
S.I.Vavilov State Optical Institute, 12 Birzhevaya liniya, St Petersburg, 199034, Russia


We present a novel technique to alleviate the problem of the global tilt in artificial guide stars for adaptive optics. This technique is based on the registration of trails of radio-excited plasma spots caused by the atmospheric tilt. Following the time trace of the trails one can find and measure the tilt produced by atmospheric turbulent layers. Different methods were applied to estimate the extent of the trails. We describe results of computer simulations, showing the performance of the proposed approach.
In the past several years a number of 8-10 m ground-based telescopes have been built and more are now are under construction. All of them are still limited by the distorting effects of the Earth's atmosphere, which limits their effective coherent aperture, and are thus provided with adaptive optics. Adaptive optics allows to measure and undo the effects of turbulence in real time, delivering near-diffraction-limited performance at the infra-red and visible wavelengths. Many have detailed plans for adding artificial guide stars, which serve as reference beacons for precise measurement of the turbulence even when nearby stars are too weak for this purpose.
Unfortunately, artificial guide stars suffer from the problem of determination of the ascending beam global tilt, which is difficult to measure and remove from the descending beam. Because of reciprocity of propagation paths, a conventional synthetic beacon is unable to sense a full aperture tilt and can only be used to measure the higher-order distortions. For this reason, the problem of indeterminable tilt of artificial guide stars is recognised as a fundamental one to be overcome before applying a full adaptive optics correction for any type of seeing condition. The current solution is to use a further natural star just for measurement of the tilt, complicating the optics and processing of the adaptive optics.

Since laser guide stars (LGSs) play a key role in the achievement of full sky diffraction capabilities at visible wavelengths, the efforts of the world astronomical community in the past years have concentrated on the development of the LGS technology. However, this solution has many problems, both practical and conceptual, slowing its application. For this reason, one of authors has proposed an alternative approach to the creation of artificial reference source. This method is based on the use of visible plasma spots excited by fringes between intense radio waves (see below and in [Ribak 1997, 1998]). Three distant radio dishes are phased together to have 1 m fringes tens of kilometers above the telescope, the fringes spanning up to tens of meters. The air breaks down where the radio intensity is the highest, and creates plasma, visible in oxygen and nitrogen lines. Despite the attraction of the mentioned approach, the problem of lack of a useful tip-tilt signal is not solved in the case of radio guide stars (RGSs), just as it is not solved in the case of the LGSs. Therefore, the solution of the global tilt problem is very important for RGSs too.

The global tip-tilt problem has been attacked in several ways and now one can count more than ten techniques to sense it [Esposito 1998]. A couple of these techniques were tested experimentally on the sky. Namely, the statistical technique devised by Belen'kii [1995] was tested on Polaris using a Rayleigh laser beacon [Belenkii et al. 1999]. Shortly afterwards, a perspective based technique [Ragazzoni 1997] was tested by comparison on a couple of stars using a sodium laser guide star [Esposito et al. 2000]. Both experiments are open loop ones, i.e. no attempt has been made for real-time correction of tip-tilt, but data have been collected and later reduced in order to establish the merits of the technique. Although both experiments show a remarkable good agreement with the expected results, it is worth pointing out that tilt correlation does not exceed correction of more than seventy to eighty percent. This means that more work needs to be done in order to provide a reliable system able to produce high Strehl ratios. The Kolmogorov turbulence power spectrum is peaked at the very low frequencies, making the tip-tilt component very important. Hence only very accurate tip-tilt removal can prevent significant deterioration of the achievable Strehl ratio [Rigaut and Gendron 1992]. However, there are theoretical suggestions that in the case of favourably smaller outer scales the tip-tilt contribution can be so low that, under special but not impossible circumstances, there would be no relevant tip-tilt at all to be corrected.

Even RGSs cannot solve the problem in full, since the pointing accuracy of the radio beams is limited due to their large wave length. In addition, atmospheric effects are estimated to swing the interfering beams at about 0.1 m at the elevation of 50 km. In addition, local winds up to 40 m/s might carry the plasma spot away, although its elevation can be chosen to be minimally affected by the wind. To make things easier, since the pulse repetition rate of the radio beams is chosen to leave residual plasma for the next pulse (to save on power), this repetition rate can be tuned to have a known spot drift. See the book by Gurevich et al. [1997] for details on the above mentioned considerations.

We suggest a new concept for absolute tilt retrieval based on RGS technology. Because of the large scale of the radio fringes, this solution has direct implications also for the fields of radio astronomy, radio interferometry and optical interferometry, laser communications and power beaming (see below). In all these fields, there is a need for a local reference frame to help phase the elements of the radio dish, or co-phase the different dishes to serve as one coherent unit. The added capability to track tilts as proposed here will enhance these devices.


The need for artificial guide stars has led to the suggestion to use inexpensive sodium lamps, and to image them on the sodium layer. Since this image is too large to center on, interfering radio waves could quench the radiation along narrow fringes. This idea led on to proposing that the interference itself could produce light by creating a very high local electric field in the crests of the fringes (Figure 1). Such a field creates a break down of the thin air at high elevations, and the radiative recombination will be the artificial beacon [Ribak 1997, 1998]. This is a kitchen tabletop experiment – every gas lighter produces a strong electric field, with a blue spark that flows through air as a result. These photons are the result of radiative relaxation of the excited plasma. The electric field is much lower at lower air pressures, but the lines are the same – nitrogen and oxygen lines. It turns out this idea was tested before for other purposes.

Baily [1937], Gurevich [1972], Utlaut [1975] and others proposed to create an artificial ionosphere for reflection of radio waves above and behind the horizon. In the simplest scheme, a radio beam heats up and ionizes or creates a break-down in the atmosphere. At an elevation of 30-60 km a patch of plasma is thus created. Efficient reflection of the distant radio waves occurs only at a very high plasma density, so the idea was improved to have two radio beams interfere to create fringes. Radiation density culminates in the fringes where ionization starts, and the few electrons multiply very fast. The ionized fringes serve as a radio Bragg reflector, a more efficient scheme than mere return from a crude plasma volume [Gurevich 1980]. In a laboratory experiment, clear visible fringes were apparent and served to trace the plasma distribution [Kuo and Zhang 1990]. Among other applications, the use of artificial ionization was suggested as a means to create ahigh nitrogen laser (100 kW inside 1 m2) in the atmosphere, shining down on the ground [Gurevich et al. 1997]. Fringes are created between moving microwave beams, with a crossing area moving down at the speed of light. The advancing breakdown regime is a source of first spontaneous and then stimulated emission. In laboratory tests 3 kW radiation at 337.1 nm were produced [Vikharev et al. 1991]. Full description of the physical processes involved can be found in Gurevich et al. [1997], and physical and engineering aspects are described in Gaponov-Grekhov and Granatshtein [1994].

The ionosphere around 250 km [Bernhardt et al. 1988] was studied by similar methods: radio heating of a large volume, by a single radio dish, creates plasma and cavitation of the heated volume. The result is enhanced air glow, as well as a typical radio signature, over tens of kilometers. Electrons ionized and/or accelerated in the heated volume spiral out of it, along the earth’s magnetic field, and create an illuminated regime when they hit denser neutral species. However, at these elevations (above 150 km) the diffusion of the plasma and the fast winds disperse and shift swiftly any created patterns. In these experiments, relatively low power was required (10 kW - 2 MW), in a continuous mode. At lower elevations, the wavelength has to be shorter (microwave), and because its absorption decreases as more plasma is created, it is necessary to utilize pulses rather than continuous radiation.

The ionospheric air glow is very patchy. This is due to the density variations of the excited volume. In addition, filamentation was observed in the heated area, as well as in laser sodium guide stars. It is not clear what is the typical spatial frequency of the patches and what is their contrast. They have been resolved down to the ten meter scale, consistent with theoretical predictions. The visible natural intensity (night glow) is usually few tens of Rayleighs, and in the experiment values above 250 Rayleighs were also produced.

The choice of beacon elevations is important from the point of view of adaptive optics: the higher the better, provided the pattern is stable enough and of high enough contrast on the few arcsecond scale. Above 200 km the random pattern can still be used as a reference source, much in the same way that the solar surface serves for adaptive optics. Below 100 km the density of the atmosphere is high enough that the plasma mean free path (a few cm) is much shorter than the required fringe period (less than 2m [Baharav et al. 1994,6]). Thus the luminescent fringes follow the radio-heated ones. Experiments show that the most power efficient air breakdown occurs around 0.5-5 Torr [Kuo and Zhang 1990], corresponding to an elevation of 30-80 km. The sodium layer at 87-95 km is also a respectable compromise. In addition, the created plasma field can break up sodium dust particles into atoms and increase the returned signal [Wuerker 1997].

The sodium layer was powered by a 1.35 MHz transmitter with a resulting enhancement of the sodium lines by 100 Rayleighs [Lyakhov et al. 1984], corresponding to mv=17/arcsec2. No details were given for the actual power used or the spot size created, except that it was measured from a distance of 90 km, at an elevation of 80-100 km. Natural and laser guide stars need to be brighter than mv=12-14. To achieve this, single lasers need be brighter than 8-12 W. Laser fringe systems, which measure more than one atmospheric layer, require 200-500 W [Baharav et al. 1994,6]. However, these estimates can be reduced drastically by employing an optical fringe demodulation scheme [Ribak 2000] to mv=11-13 (11-16 W). This is close to the theoretical limit, set by the number of degrees of freedom, or r0 patches, to be measured. In two conjugate layers there are only 20-30% more r0 patches than one integrated volume. Thus by reducing the plasma spot size by a factor of one hundred with respect to Lyakhov et al. [1984], the integrated surface brightness could reach the required level for a full sky coverage. This reduction in size can be achieved by employing phased arrays and interfering beams, and wave lengths shorter by orders of magnitude.


Tilt, produced by atmospheric turbulent layers, can result in a trail of the spot of the artificial guide star visible from the telescope. In considering this effect, we pose a question: is it a positive or a negative effect and can it be used to break down the reciprocity and to sense a full aperture tilt with an artificial guide star?

To reply this question, let us suppose that the upward interfering beams exciting the beacon are characterised by an upward wandering b(t). In what follows, we only consider the central fringe; inclusion of more fringes (which move in unison) can only improve the method discussed below. Let us also suppose that the excited spot has a certain decay time t0. In other words, if one excites the beacon with a single short pulse of intensity I0 , one will observe a beacon intensity that changes with time accordingly to:

I = I0 exp(-t/t0)

For the mesospheric sodium layer, excited by laser light at the same D1 and D2 lines [Happer el al. 1994], t0 is of the order of pSec. The effect of the trails caused by turbulent atmosphere tilt is negligible and hence useless. In exciting the oxygen or nitrogen in a plasma cloud using radio waves, the decay time t0 can take much longer values. Because the recombination of plasma varies according to the density of the atmosphere, the decay time will depend upon the height of the excitation. At the altitude of 60 km the decay time of plasma is of the order of 0.01 sec, comparable with the time required for the excited beam to be registered. The repetition rate is chosen so as to leave exactly this residual plasma for the next pulse. Hence, we will see simultaneously the recently excited spot and some residue of previous spots next to it. As a time-series, this looks like a bright spot with a trail behind it. This is similar in appearance to an over-exposed image intensifier showing trails behind bright moving spots. Therefore, in the case of RGS the trails of the excited plasma spot, produced by tilt, will be visible from the ground-based telescope.

Since the elevation of the plasma can tuned, t0 will vary with it, for instance, to be of the same order of the Greenwood frequency. However, there is no need to have t0 locked to some specific frequency related to the atmospheric parameters.

Except for an ionospheric beacon [Ribak 1997, 1998], the spot and its trail are always contained within an isoplanatic patch. In this situation, we will see the evolution of the beacon. Because of the non-zero decay time t0, we will see a sort of evanescent trace that follows the beacon movements. The direction and extension of portion of the trail closer to the beacon is an indication of the variation of the upward tilt. A slowly moving upward tilt will produce a short trail, a fast moving upward tilt will produce a longer trail. For instance, looking at the brightest point in the trail (the current beacon), of brightness I, and at an adjacent point of brightness I/e, one will see exactly an angular separation of

s(t) = b(t) - b(t-t0)

This is, roughly speaking, the derivative of b(t). More precisely, this is a finite difference of b(t) with respect to t0 (Figure 2). In any case, by integration or numerical summation one should be able to retrieve b(t). Each measurement of the finite difference will be affected by some error and the integration error will grow with time leading, at a certain time, to some unacceptable value. If this happens after long time (compared to t0), one can lock onto a much fainter natural guide star and close the low frequency loop there. Such an approach was discussed in detail previously [Ragazzoni 1996].

Using the assumption that b(t) is always within an isoplanatic patch, the measurement of the trail extent and position angle is totally unaffected by the downward tilt, because both b(t) and b(t-t0) are observed at the same time and are in the same isoplanatic patch. Hence, this downwthas no effect on the differential evaluation of the derivative of the upward tilt.

Therefore, the full aperture tilt of the atmosphere can be found and measured following the time trace of the trails of the radio beacon. It means that the proposed technique permits the sensing of the full-aperture tilt of the atmospheric wave front distortions using a RGS, eliminating one of the fundamental limitations of adaptive optics.


We wish to find out the required intensities for our scheme to be applicable. The (linear) size of each fringe is , where is the radio wave length and b the base line. The centroiding error is , where N is the number of photons. The number of measurements is the ratio of the total measurement time to the single frame integration time t0 (which is the smaller of the Greenwood rate and the plasma decay rate). As the number of measurements grows as , this centroiding error grows as . The number of photons in a single integration, N, is related to the number of photons per second, N0, through N=N0 t0. The total centroiding error is


and it has to be smaller than ([Sivaramakrishnan et al. 1995]). Solving for the number of photons per second, N0, we get


To get the required brightness in Rayleighs we divide this number by 1010 photons for MKSA units, by the aperture area and by the fringe angular area:

For the sodium line, an integration time t0=0.01s and total integration time of =0.1s, we get 200,000 Rayleighs from a single fringe crest, much in line with the estimates in Section II. This integration time should allow for full sky coverage provided that a faint NGS in the isoplanatic patch is used to corroborate the measurement [Ragazzoni 1996]. Inclusion of all the crests will reduce this requirement even more. This is because the number of fringes in each dimension is proportional to the ratio of the base line b to the single aperture diameter d. For (b/d)2 fringe crests the error will drop by b/d and the total plasma pattern area will grow by (b/d)2, with an attendant drop in the necessary flux by a factor of (b/d)3.


To check out the proposed approach we ran many computer simulations. Specifically, we generated a series of RGS spots using the Kolmogorov spectrum for the upward tilt, common to the three radio beams involved. The drift in the tilt (and hence the trail) of the radio beams can be as small as 0.1 m, in which case it might be so stable that it does not require any improvement [Gurevich et al. 1997]. Even that small a drift can be detected by extending it intentionally, by adding an oscillatory angular term to the upward beams. Since the radio beams are interfered to create bright fringes, all one has to do is change the phase of the beams relative to each other. Notice, however, that the shift of the spot cannot be too large, since each radio pulse re-vitalises the plasma from the previous pulse, this being more power efficient than a continuous excitation [Gurevich et al. 1997]. The size of the plasma spots depends on the chosen wave length, elevation, and distance between transmitters, and can range between a few decimeters and a few decameters [Ribak 1997, 1998, Gurevich et al. 1997]. We used a gaussian spot profile, of standard deviation varying between 75 and 200 cm.

In the simulation, we created a single plasma spot; for three transmitters one would get at least three, around the telescope, but we only wish to prove the method here. The next shot will create a new nearby spot, while the older one is reduced by a decay factor, chosen between 0.1 and 0.7 (i.e. with decay times different from t0). The next shot will add yet another new spot, with the two older ones decayed again by the same factor. Some snapshots of these simulations are shown in Figure 3. Between shots, we moved the beam using a Kolmogorov spectrum of tips and tilts, with a preset standard deviation varying between 1 and 3 m. In analysis, we added Poisson noise to the detected signal: we assumed between 20 and 2000 photons at the brightest pixel. All calculations were done in terms of pixels only, where each pixel is equivalent to about 25 cm at the plasma elevation, or 0.5-1 arcseconds. Other noise sources, such as detector additive noise, or non-homogenous plasma, were neglected. Consequently we now have a sum of two gaussians, one strong and one weak, and even weaker residues of older gaussians (Figures 2, 3). We have to identify the gaussians, so as to detect the most recent one and the distance from it to the previous one. This is a difficult task, even without including noise in the signals. We tried a number of ways to analyse the data.

(a) We first identified the brightest point in the pattern. This is the approximate location of the current spot. We subtract a gaussian from the intensity pattern, and locate the previous spot as being the next brightest point again. This is similar to the CLEAN algorithm used to deconvolve radio maps [Thompson el al. 1986, Perley el al. 1985], except that we stop after only two iterations. The resolution is limited to integer pixels.

(b) Next, we tried to use optimisation: we fitted one gaussian near the brightest pixel, and subtracted it as before. The older spot was located by finding the center of gravity of the residual.

(c) The first two methods suffer from bias in the choice of the first spot, a bias introduced by the nearby second spot. To avoid this bias, we fitted two gaussians simultaneously. The optimisation process involved in the fitting is rather slow.

(d) Another way to tell the new and old spots apart is by comparing the first moment of the intensity pattern and the first moment the square of the intensity pattern. Because the older spot is weaker, it will show less on the squared image. Let d be the decay ratio, and the first and second intensity moments be i1 and i2 . Call for short the estimates for the previous, current, and separation , a, b, and s = b - a, respectively. Now we get

s = (i1 - i2)/[d2/(1 + d2) – d1/2/(1 + d1/2)]

b = i1 - s d2/(1 + d2)

(e) In an alternative analysis, we found the first and second coordinate moments of the elongated spot. Assuming a gaussian profile for each spot, we found the analytical values of these moments along the two Cartesian axes. Thus we could find the locations of the two recent spots and hence their separation. Including the location of the brightest spot in the pattern broke the central symmetry of these moments. Let v be the variance of the gaussian spot, x1 and x2 the first and second moments, and x3 the location of the brightest spot. The old and new spots are at

a = x1 - sign (x3- x1) [(-x12 - v + x2)/d]1/2

b = x1 (1 + d) - a d

(f) Since the spots tend to merge, it is possible to sharpen the images digitally. We used a Wiener deconvolution scheme, in which the theoretical shape of the spot was used as the object, and the noise has Poisson statistics [Christou el al. 1986, Ribak 1986]. In the Fourier domain the deconvolved image is

Ideconv = Isource S* / (SS* + N),

where Ideconv, Isource and S are the Fourier transforms of the deconvolved image, the original image, and the spot profile. N is the mean power spectrum of the noise, and * symbols conjugation. After deconvolution, the brightest and next brightest spots are identified, as in (a). The accuracy is limited again to integer pixels.

After identification of the two last spots in the pattern, we calculated their distance. We integrated this quantity for a large number of pulses. Although we have quite good knowledge of the statistics of the beam wander, we did not include them in the calculation as a prior. This allows to include other sources of tilt, such as mechanical vibr, as well as intentional scanning of the beams. In the simulation, we used the known gaussian spot shape for fitting it in the different meth. In the experiment, this profile will have to be measured.

In Figure 3 we see the true location of the trail as created by the radio beacon in the sky, and the calculated location from the data reduction using the different schemes. As shown next, the method works best when the trail extent is the longest and the spots are nearly separated, when the decay in time is the least, and when the noise is negligible. Identification of the brightest location is always quite accurate, but that of the previous one is lacking. This is the main source of error. We should remark that discretisation errors, in some of the methods, grow with the square root of the integration time. Still, we could see that the numerical integration did follow the actual beam with time even for quite unfavourable conditions.

We compared the six methods for different decay values. In this case, the gaussian spot had a standard deviation of 3 pixels, and the jitter a standard deviation of 4 pixels. The mean signal at the brightest pixel was twenty photons, and the simulation included 200 steps. If the true and estimated translation vectors (between old and new spot, Figure 4) at step i were ti and ei, what is shown is [| ti - ei |2 /?i]1/2. The results are shown in Figure 5: the two methods with the least error (b and d) are also the methods that tend to slightly underestimate the translation vector - the other methods slightly overestimate it. Corresponding results were obtained for different parameters of the beam.VI PARAMETERS OF THE METHOD

The process of the creation of trails of RGSs may be more complicated than simulated here, since atmospheric wind dispersion and general diffusion can also result in this effect. These processes tend to spread even the single beacon [Gurevich et al. 1997]. In particular, the maximum speed of winds at altitude of about 60 km (where inter-pulse time is about 10 ms) and latitudes between 30 degrees north and 30 degrees south (where most telescopes are located) can reach a reasonable value of about 20 m/s. Tuning of the plasma to an altitude of slower wind is possible. However, in this study we have not introduced any influence of wind on the beacon layer for simplicity. We believe that wind, if constant in direction and intensity, will introduce an additional apparent constant value in the derivative estimate. This can be removed easily in the real-time reduction process, because it is a constant offset in the derivative, estimated at each update loop on the faint natural guide star.

How does the distance between the old and new plasma spot affect our results? Of course, the further the points are the easier they are to separate numerically. If this separation is too small, we can enhance it by scanning the relative phases of the interfering beams, and thus moving the fringes across the interference pattern. Unfortunately artificial sweeping of the spots, or by atmospheric propagation phase and local wind, cannot exceed a fraction of the spot size: beyond this distance, the seed of ionised air from the last pulse will not be enough to initiate the new breakdown. This is especially important at lower elevations, where recombination is fast and complete. In the simulation we replaced intentional sweeping by random tilts, created by pure atmospheric spectrum of intentionally large extent.

The effect of noise was also investigated. Since there are only rough estimates for the number of photons available (Sections II, IV and Gurevich et al. [1997]), we included in the simulation from a few to an infinite number of photons. We found out that down to 20 photons (at the brightest pixel) we could still locate the current and previous beacon, but below that the signal deteriorated too fast. In addition, Poisson noise caused over- or under-estimation of the locations and translation (difference) vectors, because this noise was not included in the simple model of the system and the resultant estimator methods.

The effect of decay of the plasma between pulses is also important. Since there are again only rough estimates for this decay, we assumed the intensity of the previous pulse is between 0.1 and 0.7 of the recent one. The stronger the decay, the easier it is to find the location of the two spots. This is partly because in all six methods we assume only two such spots, and neglect the rest of the history. However, when the signal is too weak, a strong decay means low signal-to-noise ratio for the previous spot.

As a result of these simulations, we can say that for each elevation one should be able to find the optimal choice of parameters. At low elevations (30-50 km), the decay is rapid, so the pulses should be close. If the overlap is then too much, constant wind or artificial sweeping are required. Further up more power is needed, the plasma decays more slowly, and the pulses can be spaced apart, but the weaker signal might be a limitation. Above 100 km the plasma pattern is set by other parameters [Gurevich et al. 1997] such as magnetic field lines and is much less controllable.


We showed that the scheme of tracking the trails of plasma spots in order to reduce the tilt problem is indeed valid for the simple assumptions made here. We tested a number of methods to reduce the data and showed that it is possible to follow the tilting beam with time. The most accurate method was fitting of a gaussian to the plasma spot, and calculation of the center of gravity of the residual. Better results are possible by extension of the fit to even older spots, and by fine-tuning of the chosen parameters for the fit.

We have presented a novel technique that allows solving the global tilt problem in adaptive optics. This technique can be exploited not only in optical astronomy, but also in long base line optical interferometry [Gavel 1998], laser communications and laser power beaming to satellites. Another important application is in radio astronomy and in radio interferometry. In radio astronomy one needs to phase segmented mirrors and to sense the wave front [Serabyn et al. 1991]. One of the most severe problems in radio interferometry is the lack of good phase calibration between the dishes [Thompson el al. 1986, Perley el al. 1985]. The existence of reference hot spots is a significant step towards solving this problem. As in the visible regime, it should be set outside the central lobe of the radio dish or interferometer, but still within their field of view.


V R Baily (1937), Resonance in interaction of radio waves, Nature 139, 68, 838.

Y Baharav, E N Ribak, and J Shamir (1994,6), Atmospheric tomography using a fringe pattern in the sodium layer, Optics Letters 19, 242-4; Wide field analysis of turbulence layers using fringes in the mesosphere, Journal of the Optical Society of America A 13, 1083-1097.

M S Belen'kii (1995), Full-aperture tilt measurement technique with a laser guide star, SPIE 2471, 289-300.

M S Belen'kii, S J Karis, J M Brown, and R Q Fugate (1999), Experimental validation of a technique to measure tilt from a laser guide star, Optics Letters 24, 637-9.

P A Bernhardt, L M Duncan, and C A Tepley, (1988), Artificial airglow excited by high-power radio waves, Science 242, 1022-7.

J C Christou, E K Hege, J D Freeman, and E Ribak (1986), A self-calibrating shift-and-add technique for speckle imaging, Journal of the Optical Society of America A 3, 204-9.

S Esposito (1998), Techniques to solve the tilt indetermination problem: method, limitations, SPIE 3353, 468-76.

S Esposito, R Ragazzoni, A Riccardi, C O'Sullivan, N Ageorges, M Redfern, and R Davies (2000), Absolute tilt from a Laser Guide Star: a first experiment, Experimental Astronomy (in press).

A V Gaponov-Grekhov, and VL Granatstein, eds (1994), Applications of high-power microwaves, Artech House, Boston.

D T Gavel (1998), Wide-baseline optical interferometry with laser guide stars, SPIE 3350, 793-9.

A V Gurevich (1972), Isothermal ionization of the lower ionosphere under the influeof radiowaves, Geomagnetizm I Aeronomiya 12, 631 (Geomagnetism and Aeronomy 12, 556).

A V Gurevich (1980), An ionized layer in the gas (ithe atmosphere), Soviet Physics Uspekhi, 23, 862-5.

A V Gurevich, N D Borisov, G M Milikh (1997), Physics of microwave discharges, Gordon and Breach Science Publishers.

W Happer, G J MacDonald, C E Max, and F J Dyson (1994), Atmospheric-turbulence compensation by resonant optical backscattering from the sodium layer in the upper atmosphere, Journal of the Optical Society of America A 11, 263-76.

S P Kuo and Y S Zhang (1990), Bragg scattering of electromagnetic waves by microwave-produced plasma layers, Phys Fluids B 2, 667-73.

S B Lyakhov, G G Managadze, M S Petrov, and I S Shlyuger (1984), Observation of pulsating photoluminescence of the lower ionosphere, synchronuous with pulsed action of a powerful radio wave, JETP Letters 38, 658-661.

R A Perley, F R Schwab, and A H Bridle (1985), Synthetic Imaging, NRAO Workshop 13.

R Ragazzoni (1996), Propagation delay of a laser beacon as a tool to retrieve absolute tilt measurements, Astrophysical Journal 465, L73-5.

R Ragazzoni (1997), Robust tilt determination from laser guide stars using a combination of different techniques, Astronomy and Astrophysics 319, L9-12.

E Ribak (1986), Astronomical imaging by filtered, weighted shift-and-add, Journal of the Optical Society of America A 3, 2069-76.

E Ribak (1997), Tomographic measurement of the atmosphere by artificial plasma fringes, European Southern Observatory Proceeding 55, 186-91.

E Ribak (1998), Alternative artificial guide stars for adaptive optics. SPIE 3353, 320-9.

E Ribak (2000), Measuring turbulence in an extended field of view. SPIE 4167-18.

F Rigaut, E Gendron (1992), Laser guide star in adaptive optics: the tilt determination problem, Astronomy and Astrophysics 261, 677-84.

E Serabyn, T G Phillips, and C R Masson (1991), Surface figure measurements of radio telescopes with a shearing interferometer, Applied Optics 30, 1227-41.

A Sivaramakrishnan, R J Weyman, and J W Beletic (1995), Measurement of the angular correlation of stellar centroid motion, The Astronomical Journal 110, 430-8

A R Thompson, J W Moran, and G W Swenson (1986), Interferometry and synthesis in radio astronomy, Wiley.

W F Utlaut (1975), Ionospheric modification by high-power HF transmitters. A potential extended range VHF-UHF communications and plasma physics research, IEEE 63, 1022.

R F Wuerker (1997), Sporadic sodium layer formation during an auroral storm, European Southern Observatory Proceeding 55, 16-21.

Figure 1: Radio guide stars are created as radio beams interfere and the air breaks down at the crests of the fringes.

Figure 2: Decay and motion of the plasma spot. On the left are shown the recent few spots, and on the right their (noiseless) appearance. The two bars and arrows are estimates from this observed pattern, where the global tilt is not known.

Figure 3: Sample simulation step: the location of the previous (cross) and current (box) beacons are the same but the decays are different (top to bottom). From left: noiseless signal, signal with noise, true locations, and results of the six estimators. The mean number of photons at the brightest pixel was twenty.

Figure 4: Calculation of the error. At each step the difference between the estimated trace and the true trace is taken. This is the current error, equivalent to subtraction of estimated and true displacement from the previous step.


Figure 5: Comparison of methods at different decay values (0.1 to 0.7 for each method), for signal up to 20 photons/pixel. The standard deviation of the beam jitter was 4 pixels, and that of the beacon size 3 pixels. Better models of the physical processes should improve all results.