All-Day Synthesis Issues

From EOVSA Wiki
Revision as of 15:53, 7 December 2019 by Dgary (Talk | contribs)

(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to: navigation, search

Correcting for Solar Rotation

Deconvolving solar images integrated over a substantial period of time presents a complex challenge due to at least two effects:

  • Solar rotation: Sources near disk center are carried from east to west at a rate of about 9 arcsec/hour. Those nearer the limb move very little in the plane of the sky. The rate of rotation is also latitude dependent due to differential solar rotation. This is made far more serious by the fact that the psf (point-spread-function, aka synthesized beam) also varies with time.
  • Intrinsic variations versus time: Even without actual flaring, evolution of active region sources in shape and or brightness can occur. In applications not involving aperture synthesis, such evolution is simply averaged out in the familiar way. However, when coupled with the changing psf due to solar rotation more serious problems can be presented.

One way to overcome these problems, when reasonably good uv-coverage is available, is simply to make images over a shorter time (e.g. 10 min, during which the rotation is less than 2 arcsec), then differentially rotate the images and average them in the image plane. There is some relatively minor problem in doing this, since the emission is distributed at varying heights and one must generally choose a reference height at which to do the rotation. Luckily, the emission near the limbs does not move in the plane of the sky, and so it is not a bad approximation to rotate emission on the disk and leave the emission above the limb unchanged. This approach is worthwhile with EOVSA when there is relatively bright active region emission that is the focus of the science study, but the fainter disk or off-limb emission is not well imaged in such short integrations.

In cases where all of the emission is weak, long integrations with EOVSA are required, and this page describes an alternative method for improving image deconvolution in that case. This approach involves creating a spatially variable psf, and hence deconvolving model sources in various parts of the solar disk with the appropriate psf for that location. Software does not yet exist for the entire process, but some exploration is possible with the current software.

Description of the Problem

If solar rotation is ignored in an 8-hour integration, say, then sources near disk center will shift 72 arcsec! If the standard psf is used to clean such sources, two errors result: (1) the source is smeared in the east-west direction, appearing larger and fainter than it should, and (2) the psf sidelobes used in the deconvolution are not those created by the moving source, so much larger residuals (and possibly false sources) are left in the image.

If the correct psf could be calculated and used, however, then applying it would reduce the residuals. Additionally the smearing of the source could be eliminated by restoring with the ideal psf core, which would also restore the true brightness. The problem, though, is that the correct psf varies with location on the solar disk.

Calculating the Correct PSF

No matter the shape of the source, assuming for now that it does not evolve with time, the smeared source is a convolution of the true source with an evolving psf that is created by changing uv sampling function of the array. In calculating the standard psf, the "point" in the point-spread-function is considered to be a unit amplitude point source at the phase center (amplitude 1, phase 0). If instead we consider a moving point source (amplitude 1, but shifting phase) when calculating the sythesized beam, then this would create a psf appropriate to that source. This amounts to creating a model visibility database with unit amplitudes, but properly calculated phases.

Consider the case of a time-dependent shift of the point source in RA and Dec, . The phase shift to be applied to a uv point with coordinate is . Thus, the visibility to be calculated is just .