the effects of the ionosphere have always been an obstacle for ground - based radio frequency observations of astronomical sources .
this is especially true for interferometric observations used to make images of objects with relatively small angular sizes .
vhf and uhf interferometers such as the very large array ( vla ) in new mexico , the westerbork synthesis radio telescope ( wsrt ) in the netherlands , the giant metrewave radio telescope ( gmrt ) in india , and the australia telescope compact array ( atca ) , among others , are all affected by the ionosphere in the same way .
the fundamental principles for the operation of an interferometer are well described in the literature ( e.g. , * ? ? ?
briefly , an interferometer measures the time - averaged correlation of the complex electric fields measured at pairs of antennas pointed at a particular object .
these correlations , or `` visibilities '' provide a measure of the spectrum of the sky brightness distribution at different spatial frequencies .
these frequencies are essentially the difference between the position vectors of the two antennas normalized by the wavelength of the observation in a coordinate system based on the position of the object in the sky .
these spatial frequencies are commonly referred to as @xmath2 , @xmath3 , and @xmath4 , and the coordinate system used is defined such that @xmath2 is the spatial frequency in the east - west direction , @xmath3 corresponds to the north - south direction , and @xmath4 gives the spatial frequency along the line of sight to the object .
thus , as the object moves through the sky , the @xmath2,@xmath3,@xmath4 coordinates of a pair of antennas , or `` baseline , '' changes .
the visibility measured for a specific baseline at a particular time is given by @xmath5 where @xmath6 denotes solid angle , @xmath7 the frequency of the observed signal , and @xmath8 is the intensity on the sky at a position given by the direction cosines @xmath9 , @xmath10 , and @xmath11 which are measured relative to the position of the observed object on the sky .
generally , interferometers are `` fringe - stopped , '' that is , the measured visibilities are multiplied by a factor of @xmath12 so that visibilities from a source at the center of the field of view will have a phase of zero ( i.e. , the source will not produce fringes ) .
this is mainly done because with fringe - stopping and for small fields of view ( i.e. , @xmath13 ) , equation ( 1 ) becomes a simple two - dimensional fourier transform such that the observed visibilities , measured as a function of @xmath2 and @xmath3 , may be converted to maps of intensity on the sky using standard numerical methods .
as signals from astronomical sources pass through the ionosphere , a phase term is added given by @xmath14 where @xmath15 is the electron density and @xmath16 denotes the path - length through the ionosphere .
thus , the phase of the observed visibility for a baseline is altered by the difference between the ionospheric phase terms observed by the two antennas , which is proportional to the difference in the total electron content ( tec ) observed along the lines of sight from the two antennas to the observed object .
to first order , if these phase terms are not removed , then during the fourier inversion process involved in making an image of the sky , the ionospheric phase terms have the effect of changing the apparent positions of objects in the image plane .
when higher order ionospheric effects begin to dominate , objects begin to appear distorted in the image plane , and in extreme cases , almost disappear .    therefore , to make an image , one must remove these phase terms through a calibration process which estimates the required phase corrections .
typically , after initial calibrations for instrumental effects are performed , one usually uses some form of a procedure referred to as `` self - calibration '' @xcite .
this involves dividing the visibilities by an assumed sky model for the observed field of view which removes any contribution by the sky brightness distribution to the observed visibility phases . following this ,
a linear fit is used to determine the complex gain for each antenna .
since the interferometer only provides phase differences between antenna pairs , the absolute phase of the complex gain for each antenna can not be determined by this fitting process .
the general practice is to choose one reference antenna for the interferometer and to set the phase for its complex gain to zero .    since for @xmath17 antennas
, there are @xmath18 baselines , this is generally an over - determined problem and can be done over relatively short time periods , depending on the brightness of the source .
one may also use this calibration to make an image , deconvolve the image to produce a better sky model , and then repeat the process until it converges .
this has been shown to be a rather robust procedure for determining the complex antenna gains ( see * ? ? ? * and references therein ) .    according to equations ( 2 ) and ( 3 ) , the phase corrections obtained from the determined antenna gains are essentially measurements of the difference between the tec along the lines of sight of a particular antenna and that of the reference antenna .
these equations also demonstrate that the effect of the ionosphere will be more substantial for vhf observations . given the size of available astronomical vhf interferometers ( baselines ranging from @xmath19 km to as large as @xmath20 km ) , the robustness of the self - calibration procedure on relatively small time scales ( typically @xmath21 minute , but as small as a few seconds for extremely bright sources ) , and the sensitivity of such interferometers to relatively small tec variations ( fluctuations in differential tec of as small as 0.001 tecu ) , these instruments are capable of studying tec fluctuations on substantially finer scales than many others available .
subsequently , previous work has been done using astronomically motivated observations ( e.g. , * ? ? ?
* ; * ? ? ?
* ) as well as observations geared toward studying the ionosphere ( e.g. , * ? ? ?
* ; * ? ? ?
* ) to explore the phenomenology of the ionosphere on these fine scales .
much of this work has been performed with the vla ( latitude@xmath22 n and longitude@xmath23 w ) .
this is in part due to the fact that the vla is relatively well suited to the study of the ionosphere because its 27 antennas are distributed in a `` y''-shape which allows it to probe structures along three different directions .
arrays such as the wsrt and atca are so - called `` east - west '' arrays because the antennas are aligned along a single east - west axis and can therefore only observe ionospheric fluctuations along one dimension .
the vla is also unique in that the antennas are moved from time to time among four different configurations , referred to a , b , c , and d. the configurations go from larger and more spread out to smaller and more compact . at its largest in the a configuration , each of the array s three arms is about 20 km long with the shortest antenna separations being about 1.5 km . in the most compact configuration ,
d , the arms are at most 0.6 km long and the shortest spacings are about 0.04 km .
this allows the vla the ability to study the ionosphere over a wider range of physical scales than other similar interferometers .
finally , the vla had a somewhat unique vhf system in place that allowed observations to be made simultaneously at 74 mhz and 327 mhz using a pair of dipole antennas ( one for each band ) mounted near the prime focus of each antenna .
recently , the vla electronics and receivers were upgraded to establish the new expanded vla ( evla ) which does not include the old vla vhf system .
however , a new and improved vhf system is now being developed and will be available in the near future .
past work using the vla and other similar instruments has led to interesting results .
these include discoveries such as the new class of magnetic - eastward - directed ( med ) waves , predominantly found at night , discovered by @xcite as well as larger statistical studies such as the measurements by @xcite of the dependence of differential ionospheric refraction on relatively large angular scales ( @xmath24 ) using data from a 74 mhz all - sky survey which showed large dependences on time of day .
however , there is much more to be learned from these types of data , especially at smaller time and amplitude scales .
therefore , we are embarking on a program utilizing the vla data archives ( https://archive.nrao.edu ) that seeks to push this type of analysis to even finer scales .
we will use previously unexplored data sets of the brightest vhf objects .
we will apply new techniques for calibration and the mitigation of radio frequency interference ( rfi ) to other data sets to significantly improve their sensitivity to small amplitude tec fluctuations as well as fluctuations occurring on smaller time and spatial scales than could be explored previously .    here
, we describe the first step in this program , a thorough evaluation of the ionospheric information contained within a single , relatively long vhf observation of one of the brightest radio sources in the sky with the vla .
this data set provides the opportunity to develop and establish techniques for processing , analyzing , and interpreting similar data in future segments of our program .
this paper focusses on the data selection and calibration as well as the post - processing done on the phase information extracted from this exemplar data set to measure tec gradients with the array . in a companion paper @xcite ,
we detail new techniques for spectral analysis of these data .
in the study presented in this paper , we have sought to explore the kinds of ionospheric phenomena capable of being observed with the vla at the lowest fluctuation levels and the smallest times scales possible . to ensure that we were able to determine ionospheric phase terms to the the level of accuracy needed at the smallest available time intervals , including the resolution of any @xmath25 ambiguities
, we chose to use a data set focused on cygnus a ( or `` cyg a '' ; also known as 3c405 ) . with a total intensity of more than 17000 jy ( @xmath26 w m@xmath27 hz@xmath28 ) at 74 mhz ,
cyg a is one of the two brightest sources in the sky at frequencies below 100 mhz .
its brightness distribution on the sky is also known relatively well , implying that it can be used within the self - calibration procedure described briefly in  1 to solve for the phase corrections ( both instrumental and ionospheric ) for each antenna over relatively short time intervals .
the data set we have selected consists of simultaneous 74 mhz and 327 mhz dual polarization observations of cyg a with the vla over a period of more than 12 hours on 12 - 13 august , 2003 ( vla program number ak570 ) . during this period , there was a moderate amount of geomagnetic activity ( @xmath29 index @xmath304 ) and low solar activity ( @xmath31 sfu ; 1 sfu@xmath32 w m@xmath33 hz@xmath28 ) .
for this data , the vla was in the a configuration with the vhf dipole system available on all but three antennas ( antenna numbers 11 , 13 , and 15 ; see  1 and fig .
[ layout ] ) .
the exact layout of the antennas is shown in fig .  [ layout ] . the data set
contains a 35 minute scan ( i.e. , a contiguous block of observing time ) , and a longer , nearly 13 hour scan which , taken together , cover a range of local time from roughly @xmath34 on 12 august to @xmath35 the next morning . for the 74 mhz band ,
the total bandwidth was 1.5 mhz ; it was 3 mhz for the 327 mhz band .
for both bands , the time - averaged visibilities were measured over intervals of 6.67 seconds simultaneously in both rr and ll polarizations .
the processing of vhf interferometric observations from the vla or other similar interferometers is , in the average case , a lengthy and difficult process .
one usually must deal with significant sources of rfi , which typically appear much stronger on shorter baselines . for the vla ,
the field of view at 74 mhz is more than ten degrees in diameter , and the ionospheric conditions may vary significantly from one part of the field of view to another , which requires special calibration techniques .
one must also typically produce an image for the entire field of view , regardless of which object(s ) one is interested in because the removal of sidelobes ( i.e. , secondary peaks of the impulse response in the image plane ) produced by other objects is crucial to minimizing noise in the final image .
fortunately , the unusually high intensity of cyg a makes it substantially brighter than any source of rfi or any other object in the field of view .
this makes the calibration process much simpler and more straight forward .
the calibration of the data was performed using standard tasks within the astronomical image processing system ( aips ; * ? ? ?
* http://www.aips.nrao.edu ) .
the first step in the calibration was to determine the bandpass response of each band ( aips task ` bpass ` ) .
this was done within one minute time intervals by dividing the data by the visibilities measure within a few central channels where the response typically peaks .
the relative amplitude and phase responses across each bandpass were then measured and interpolated onto the full set of time steps . after trimming the first and last few channels from each band where the response drops substantially , the data were corrected for the bandpass responses .
following this , at 6.67 second intervals ( i.e. , the shortest possible for these data ) , the data for each band and polarization were used with model visibilities computed using images of cyg a presented by @xcite ( currently publicly available at http://lwa.nrl.navy.mil/tutorial/ ) to compute the phase corrections for each antenna at each time step .
this was done with the aips task ` calib ` which does a series of consistency and sanity checks as it determines the solutions needed to fit the measured visibilities to the model and flags antennas and intervals that appear spurious or of poor quality .
we relaxed some of the criteria for these checks , namely the minimum number of antennas ( we used four ) required and the minimum signal to noise ratio ( we used three ) because of the relatively short time intervals used . as noted above , an interferometer like the vla can only measure relative phases , implying that for this calibration to work , a reference antenna must be chosen whose phase is arbitrary and is subsequently set to zero .
typically , the antenna closest to the center of the vla is not used because the effects of rfi tend to be the worst for this antenna . however , as noted above , cyg a is bright enough that this is not a consideration for our data set , and we therefore used this antenna as our reference ( antenna 9 ; see fig .
[ layout ] ) .      following the determination of phase corrections using standard aips routines ,
several steps were taken to extract ionospheric information from the phase data .
this was done using ad hoc , python - based software .
the phase corrections measured by the calibration process contain contributions from several effects , the ionosphere being the largest , especially at 74 mhz .
in particular , the phase difference between two antennas within a single baseline is given by @xmath36 where @xmath37 denotes the difference in the ionospheric phases along the lines - of - sight of the two antennas given by equation ( 3 ) , @xmath38 represents difference in the instrumental effects of the two antennas , @xmath39 is the contribution to the phase difference from the structure of the observed source , and @xmath40 is the contribution from @xmath25 ambiguities
. we can remove @xmath39 by dividing the observed visibilities by a model of cyg a and @xmath40 by having short enough time sampling to `` unwrap '' the phases ( see below ) .
the dual 74 and 327 mhz observing mode is very useful for removing instrumental effects since the ionospheric phase simply scales with wavelength [ see equations ( 2 ) and ( 3 ) ] , whereas instrumental effects do not necessarily .
the instrumental components of the phase corrections include errors in the delays added to individual antenna signals used to fringe - stop the visibilities ( see  1 ) and offsets between the antenna pointing and the actual source position .
these two effects and @xmath25 ambiguities were dealt with in three separate steps .
first , the aips task ` calib ` flags antennas that are spurious at each time step so that the phase corrections are `` missing '' for some antennas at a relatively small fraction ( @xmath41 ) of time steps . generally , a small amount of missing data is not problematic .
however , to facilitate the use of fast fourier transforms ( ffts ) within our spectral analysis of the data presented in a companion paper which work best with evenly sampled data , these missing time steps were filled in .
we have done this by linearly interpolating the real and imaginary parts of the complex antenna gains , @xmath42 , computed by ` calib ` for all antennas onto a common time - step grid consisting of 7299 steps spaced by 6.67 seconds , and then recomputing the phase corrections \{@xmath43$]}.    next , the time sampling of the data ( 6.67 seconds ) was sufficiently short that the phase , @xmath44 , as a function of time for each antenna could be `` unwrapped '' in the conventional way , i.e. , by correcting phase jumps of more than @xmath45 radians by adding or subtracting @xmath46 .
there was one caveat to this process , however . for each antenna
, there were a few ( ranging from zero to five ) times steps where , for one reason or another , the phase correction was either spurious or represented a real and very short jump in the instrumental phase , appearing as sharp spikes in the unwrapped phases .
since there were a total of 7299 time steps , a few short jumps in phase would not be an issue except for their effect on the unwrapping process .
any of these spikes can cause an artificial large phase jump if it is included in the unwrapping process , which we have illustrated in the upper panel of fig .
[ despike ] .    to combat this
, we wrote a simple algorithm that computes the difference between @xmath47 at each time step and the value for the next time step , where @xmath44 is the wrapped phase .
any time step where the absolute value of this difference was more than ten times the standard deviation among all time steps for a give antenna was flagged and not included in the unwrapping process .
these empty time steps were then filled by linearly interpolating the unwrapped phase data for the un - flagged time steps . it should be noted that these spikes only occupy 12 time steps ( 6.6713.34 seconds ) and that the surrounding data are otherwise well - behaved , making interpolation a reasonable and straightforward solution to eliminating instances of such spikes .
an example of how the data were flagged is illustrated in the middle panel of fig .
[ despike ] while the resulting unwrapped phases are shown in the lower panel . from this result
, one can see that the spikes are not always completely removed from the data , but the goal of eliminating their effect on the unwrapping process has been achieved .
a number of instrumental effects can contribute to the phase corrections , including errors in the fringe - stopping process ( see  1 ) and offsets between the position of the source and that of the observed field center ( i.e. , pointing errors ) .
these effects are generally stable in time , changing insignificantly over periods of days ( see , e.g. , * ? ? ?
. however , for the vla , the instrumental phase is known to occasionally have short jumps that one must be wary of .
we have found such a jump in our data occurring at a local time of about @xmath48 on 12 august .
the jump is most obvious when one plots the difference between the 74 mhz phase and the 327 mhz phase scaled by a factor of ( 327/74 ) so that the ionospheric phases cancel out ( note , this only gives us the scaled difference between the 327 and 74 mhz instrumental phases , not the instrumental phases themselves ) .
we have plotted this difference as a function of time for antenna 14 for both polarizations in fig .
[ pjump ] to show the location of what we will refer to as the `` phase jump region . ''
while being fairly subtle in the ll polarization , it is quite obvious in the rr polarization .
the jump lasted about 6 minutes and can be seen in the data for several antennas . to deal with this
, we have treated the phase jump region , as well as the time periods before and after it , as separate scans , assuming that each scan has its own instrumental phase .
this basically resulted in us treating the data as if it contained four scans instead of two .    to remove instrumental effects ,
we have used a kind of continuum subtraction process . within this process , we have treated ionospheric fluctuations as features superimposed on a smooth continuum consisting of the instrumental phases , which vary relatively slowly with time as well as any slowly varying component of the ionospheric phase . with the current data
, we unfortunately do not have the means to separate the slowly varying component of the ionospheric phase from the instrumental effects and can therefore only measure fluctuations in tec on relatively short ( @xmath19 hour ; see below ) time scales . in the future , with instruments with larger bandwidths , it will be possible to use the wavelength dependence of the phase corrections to separate these effects since the instrumental phases are @xmath49 and the ionospheric phase is @xmath50 .
we have chosen to perform our continuum subtraction process for each antenna , band , and polarization by smoothing the unwrapped phases with a one - hour - wide boxcar which appeared to preserve any apparent fluctuations while giving a good representation of the continuum .
for the first scan and the phase jump region , we simply subtracted a single mean value from all the phases since they are each shorter than an hour .
we did the same for the first and last hour of each of the remaining two scans so that the same filter width would be used for all times steps .    following this
, we found that the position offset component of the instrumental phases presented a problem for this process near the edges of each scan . this is because , depending on the antenna , these phases can vary significantly over one hour , especially at 327 mhz . from the fringe - stopped ( see  1 ) version of equation ( 1 ) , one can see that a position offset in the direction cosines @xmath9 and @xmath10 of @xmath51 and @xmath52 will produce an additional phase of @xmath53 .
since @xmath2 and @xmath3 are normalized by the observed wavelength , any such phase will be 4.4 times larger for the 327 mhz band .
in addition , the offsets can be different for each band and this difference can vary with baseline .
this is due to a number of factors including the fact that different model images were used for the bands which may not be exactly aligned and that cyg a has a significant amount of resolved structure which larger baselines are more sensitive to , especially at 327 mhz where the angular resolution in the image plane is 4.4 times better than that at 74 mhz .
we have demonstrated this in the upper panel of fig .
[ edge ] where we have plotted @xmath54 for antenna 3 , ll polarization as a function of time . with @xmath55 scaled by ( 327/74 ) , the ionospheric phases are removed and all that is left is the difference between the instrumental phases for the two bands .
the data follow a smooth curve which is inconsistent with the known behavior of vla instrumental phases , especially at 74 mhz .
furthermore , the curve that the data follows is easily fit by a linear combination of the un - normalized versions of the @xmath2 and @xmath3 coordinates ( see the red curve in the upper panel of fig .   [ edge ] ) .
a single baseline ( antenna 3 with the reference antenna ) will sweep out an ellipse in the @xmath2,@xmath3-plane because of the rotation of the earth ( see , e.g. , * ? ? ?
therefore , this is exactly what one would expect for a scenario where there is a single position offset for each of the two bands during each scan .    to show the effect of the time dependence of the position offset phase on our continuum determination process , we have plotted the continuum - subtracted version of @xmath54 versus time in the middle panel of fig .
one can see from this plot that within the first scan and within the last hour of the second and fourth scans where a single mean continuum value was subtracted from each , the gradient of the position offset phases has introduced an artificial difference which increases / decreases with time . since the position offset phase is much larger for the 327 mhz band , we have introduced the following additional step for the continuum subtraction of the 327 mhz data . within the first scan , the phase jump region , and the first and last hours of the other two scans
, we have fit a linear combination of the un - normalized versions of @xmath2 and @xmath3 to @xmath54 separately for each time range . within each of these time periods
, we used the mean values for @xmath56 , @xmath57 , and the @xmath2,@xmath3 fit to construct a time - variable continuum for the 327 mhz data .
the benefits of this are illustrated in the bottom panel of fig .
[ edge ] where we have plotted the continuum - subtracted version of @xmath54 , this time , including the additional computation for the time - variable 327 mhz continuum with the first / last hour of each scan .
one can see that the roughly linear features seen in the middle panel of fig .
[ edge ] have been removed and that the remaining difference between the 74 mhz and 327 mhz continuum - subtracted phases is essentially noise .    following the application of the corrections detailed above , we used equation ( 3 ) to convert the continuum subtracted phases for each antenna , band , and polarization to values of differential tec , or @xmath0 .
then , at each time step and antenna , we computed the median @xmath0 among the four values ( i.e. , two bands and two polarizations ) as well as the median absolute deviation ( mad ) as an estimate of the uncertainty in the median .
to increase the reliability of the mad computations , we included with each time step the four nearest time steps ( using their own individual median values ) for a total of 20 data points per time step . in both computations ( median and mad ) , the median was used to minimize the effects of any spurious data which remained .
the resulting @xmath0 values are plotted as functions of time for each antenna in the northern arm in fig .
[ north ] , the southeastern arm in fig .
[ east ] , and the southwestern arm in fig .
[ west ] along with the mad values to illustrate the relatively high precision to which @xmath0 has been measured .
the typical @xmath0 uncertainty , represented by the mad computations , is about @xmath1 tecu , demonstrating the remarkable ability of the vla to detect extremely small tec fluctuations , even on time scales @xmath58 seconds when an object as bright as cyg a is used .
with the fully reduced @xmath0 data , including robust estimates of the uncertainties , we are in a position to explore a wide range of ionospheric phenomena .
first , we note that since the vla measures differential tec values between antenna pairs , it is essentially only sensitive to changes in the tec gradient . given the geometry of the array ( see fig .  [ layout ] ) , we can not numerically compute the tec gradient at each antenna location from our data , and measuring the full tec gradient requires a somewhat ad hoc approach .
such measurements are crucial for any analysis of observed tec fluctuations because without modification , the set of @xmath0 time series can only be spectrally analyzed for specific assumed pattern models ( e.g. , a single plane wave , see * ? ? ? * ) .
this is different from the normal mode of operation for radio interferometers in which standard techniques are used to invert and de - convolve sparsely sampled visibility data to make an image . as equation ( 1 ) demonstrates , the observed visibilities are functions of the _ differential _ antenna positions , @xmath2 , @xmath3 , and @xmath4 . for small fields of view
, the contribution of the @xmath4 term is negligible if fringe - stopping is applied ( see  1 ) .
thus , even for a y - shaped array , reasonably good @xmath2,@xmath3-coverage can be obtained .
this is improved further by the rotation of the earth which causes each baseline to sweep out an ellipse in the @xmath2,@xmath3-plane ( see  2.2 and * ? ? ?
using this fact to obtain better @xmath2,@xmath3 coverage is sometimes referred to as `` earth rotation synthesis . ''
in contrast , the gradient of an arbitrary set of tec fluctuations varies over the array as a function of the _ actual _ antenna positions projected onto the ionosphere .
improved spatial coverage can be obtained by exploiting the change in the apparent position of the observed source ( i.e. , rotation of the earth ) and the movement of the fluctuations themselves , effectively converting temporal baselines into spatial ones . however , since the tec fluctuations presumably have a distribution of speeds and directions , this is not as straightforward as in earth rotation synthesis .
one must decompose the time series into temporal spectral modes and then analyze how the properties of each mode vary across the array to extract the size , speed , and direction of the dominant pattern(s ) for that mode ( such spectral techniques are detailed in a subsequent paper ) .
therefore , we must still contend with data that has been sampled in a y - pattern which can not be inverted in a straightforward manner .
we have consequently developed two ad hoc techniques designed to provide measurements of the tec gradient time series over the full array and along each of the vla arms . before implementing either technique , we first had to perform two basic geometric corrections to the data so that the measured tec gradients would correspond to vertical tec gradients as closely as possible .
first , to ensure that our characterization of the shape of the observed tec surface is physically meaningful , we needed to project the antenna pattern displayed in fig .
[ layout ] onto the locations where the lines of sight of the antennas pass through the ionosphere , or `` pierce - points . ''
second , we needed to compute the slant - to - vertical tec corrections for the line of sight to cyg a as its apparent position on the sky changed throughout the observation . for a plane parallel approximation , both of these corrections are relatively straightforward .
however , since cyg a was as low as @xmath59 above the horizon during the observing run , a plane parallel approximation was far from valid at all time steps .
we have therefore computed the two required geometric corrections using a spherical model detailed in appendix a.    within this model , the ionosphere was approximated with a thin shell located at the height of the maximum electron density , or `` peak height '' ( see * ? ? ?
* ; * ? ? ?
we obtained estimates of the peak height as a function of time by using the international reference ionosphere ( iri ; * ? ? ? * ) software , inputting the date and time of our observations and the latitude and longitude of the vla .
we then re - determined the peak heights using the latitudes and longitudes of the pierce - points .
we found that additional iterations of this process only marginally changed the results and chose to use one iteration only .
the final peak heights used are plotted in the upper panel of fig .
[ zion ] along with the projected separations from the array center for the farthest antennas of each arm ( antennas 1 , 7 , and 22 ; see fig .  [ layout ] ) and the corresponding slant - tec corrections .      after applying the geometric corrections to the antenna positions and the @xmath0 measurements
, we sought to characterize the full two - dimensional tec gradient observed by each antenna at each time step .
rather than assume a particular dominant structure ( e.g. , a plane wave ) , we simply assumed that since the array is smaller than many transient ionospheric phenomena , the observed tec surface at any time step could be approximated with a low - order taylor series .
we examined many time steps and found that a second order , two - dimensional taylor series adequately approximated the amount of curvature in the tec surface detected by the vla .
this taylor series has the following form @xmath60 where @xmath16 and @xmath61 are the north - south and east - west antenna positions , respectively , projected onto the surface of the ionosphere at the estimated peak height . to maximize the amount of data used to constrain the parameters of each fit , we used the difference between @xmath0 for each of the 300 unique antenna pairs at each time step
. thus , the form of equation ( 5 ) actually fit to the data was @xmath62 where the @xmath63 and @xmath64 subscripts denote the values for the @xmath65 and @xmath66 antennas , respectively .
we also utilized some standard sigma - clipping during the fitting process for each time step by computing the rms of the fit residuals , rejecting all antenna pairs with absolute residuals @xmath67 , and repeating 50 times .
as many as about 10 and as few as zero were rejected for any given time step .
we note that each time step was fit independently to preserve the presence of any small - scale spatial / temporal fluctuations .
the fitted coefficients as a function of time are plotted in fig .
[ pfit ] . from these
, one can see the same large amplitude and period fluctuations at the beginning of the observing run that are visible in the individual antenna data plotted in fig .
[ north][west ] .
note that they are not quite as large here because of the applied slant - tec correction discussed above . here
, we can see that they are most visible in the @xmath68 coefficient which is the partial derivative of the tec surface at the center of the array along the north - south direction . with the plots in fig .
[ pfit ] , one can also see the same thing beginning to happen near the end of the run toward dawn .
this is qualitatively consistent with the known behavior of medium - scale traveling ionospheric disturbances ( mstids ) which are prevalent near sunrise and sunset ( e.g. , * ? ? ?
during the middle of the night , the second - order terms become more significant .
while the polynomial - based measurements provide useful information about the variation of the full two - dimensional tec gradient , they neglect the ability of the vla in its a configuration to detect fluctuations on scales as small as a few kilometers . in principle , one could do this by simply increasing the order of the polynomials used .
however , it is likely that the small - scale structures observed do not span the array .
this implies that such fits would not yield accurate representations of the full tec gradient at each antenna , especially those near the ends of the arms ( see fig .
[ layout ] ) .
therefore , instead of using higher order polynomial fits , we have opted for an alternative approach to make full use of the data .
this complementary method computes the projection of the tec gradient ( or , the spatial derivative of @xmath0 ) along each vla arm .
the projected gradient was computed at each time step separately for the antennas of each arm using simple three - point lagrangian interpolation .
given the typical @xmath0 precision of @xmath1 tecu and the mean separation between antennas of 2.5 km , the precision of these projected tec gradient measurements is typically about @xmath69 tecu km@xmath28 . the time series of the projected gradient at each antenna is plotted in fig .
[ tecgrad ] .
we have also plotted in red the projected gradient computed using the polynomial coefficients plotted in fig .
one can see that for the larger amplitude , longer period disturbances , the polynomial fits largely recover the structure observed using the data for individual antennas .
however , during the middle of the night , there appears to be a significant amount of smaller - scale structure missed by the polynomial fits that can only be observed using the individual antenna gradients , especially for the shortest baselines near the center of the array .
our exploration of a long , vhf observation of cyg a with the vla has successfully demonstrated the power of this instrument to characterize a variety of transient ionospheric phenomena .
for this observation , the typical @xmath70 uncertainty in the @xmath0 measurements was @xmath1 tecu , yielding more than an order of magnitude better sensitivity to tec fluctuations than can be achieved with gps - based relative tec measurements ( see , e.g. , * ? ? ? * ) .
large amplitude , long period waves reminiscent of mstids are visible within the @xmath0 data near dusk and dawn as well as other times intermittently throughout the night .
the polynomial - based approach we have detailed in  3.2 appears to be able to recover the properties of the full two - dimensional tec gradients associated with these relatively large disturbances as they passed over the array .
this information can be used to estimate the size , speed and directions of such patterns down to scales of roughly half the size of the array ( @xmath71 km ) .
this is demonstrated in more detail in a subsequent paper describing the spectral analysis of these data .
in addition , our approach of measuring the projected gradient at each antenna along each arm has shown that there are smaller - scale tec fluctuations observed throughout the night , most prominently after midnight local time .
thus , the vla can be used to simultaneously study fine - scale ionospheric dynamics .
this may include a host of phenomena such as the small - scale distortions / structure within mstid wavefronts , turbulent fluctuations from ion - neutral coupling within the lower ionosphere / thermosphere , and sporadic - e ( @xmath72 ) layers . in the case of @xmath72
, @xcite demonstrated with a combination of vla data and optical observations that many of the small - scale fluctuations seen by the vla during summer nighttime are likely associated with @xmath72 layers .
@xcite showed that the tec gradients caused by these layers are typically @xmath73 tecu km@xmath28 which is easily detectable using the arm - based gradient method .
a specialized spectral analysis technique has also been developed for these measurements and will likewise be detailed in the companion manuscript to this paper .
two basic geometric corrections must be applied to the antenna positions and @xmath0 measurements so that they more accurately represent the actual conditions within the ionosphere .
since we have used observations of cyg a that include times when it is relatively close to the horizon , we can not use a plane - parallel approximation .
instead , we have used a thin shell approximation for the ionosphere where the shell is located at the height of maximum electron density , @xmath74 , as computed by the iri software for the dates and times of the observations ( see  3.1 and fig .
[ zion ] ) .
the full spherical corrections used are detailed below .
first , the positions of the antennas on the ground must be converted to projected positions within the ionosphere which , for a non - plane - parallel atmosphere , change with the elevation of the observed source . for a spherical shell , we may define a `` pierce - point '' for each antenna where its line of sight to the source passes through the ionosphere .
the positions of these pierce - points relative to that for the center of the array can then be used as their projected ionosphere positions .
[ geo ] provides a schematic illustration ( not to scale ) of how these positions are determined .
we first define a set three position vectors , @xmath75 , @xmath76 , and @xmath77 , which define the positions of the array center / antenna , the ionosphere pierce - point , and the observed source , respectively , with the center of the earth as the origin of the coordinate system .
next , we note that the vast majority of astronomical sources , including cyg a , are essentially infinitely far away , which implies that the line of sight from the array center / antenna location to the source is essentially parallel to that from the center of the earth to the sources , or @xmath78 . if we define a `` left - handed '' coordinate system such that @xmath75 points along the @xmath79-axis , then the source position is given by @xmath80 where @xmath81 is the angular elevation of the source , @xmath82 is the azimuthal angle measured from north though east , and the @xmath83 and @xmath84 unit vectors point toward the north and east , respectively , as viewed from the array / antenna . combining this with the assumption of parallel lines of sight to the source and the fact that the length of @xmath76 is set to @xmath85 yields the following expression @xmath86^{2 } \nonumber   + \left [ r \ ; \mbox{cos}(h)\mbox{sin}(a ) + y_{a } \right ] ^{2 } \nonumber   + \left [ r \ ; \mbox{sin}(h ) + z_{a } + \mbox{r}_{\mbox \scriptsize earth } \right ] ^{2}\ ] ] where @xmath87 , @xmath88 and @xmath89 are the coordinates of the antenna relative to the array center and @xmath90 . since the antenna positions are known , @xmath91 is the only unknown variable .
equation ( a2 ) can then be rewritten as a quadratic equation and solved for @xmath91 keeping in mind that @xmath92 which allows one to compute the @xmath16 , @xmath61 , and @xmath79 coordinates of @xmath76 for the array center and each antenna in the current coordinate system . following this ,
a coordinate rotation was performed such that @xmath76 for the array center pointed along the @xmath79-axis and the @xmath16 and @xmath61 axes pointed toward north and east , respectively , as viewed from the location on the earth directly below the array center pierce - point .
these rotated coordinates were then taken to be the @xmath16 and @xmath61 antenna positions projected onto the ionosphere thin shell for each time step .
[ geo ] shows a graphical representation of these computations for the array center ( in black ) and for an exemplar antenna ( in gray ) .    the second correction deals with the fact that the path length through the ionosphere is increased when the observed source is closer to the horizon .
for a thin spherical shell , it is increased by a factor of @xmath93 where @xmath94 is the angle between the line of sight from the vla to the source and a line from the ionosphere pierce - point to the location on the earth directly below it . in the schematic in fig .
[ geo ] , @xmath94 is the angle between the position vectors @xmath76 and @xmath95 .
therefore , to compute the factor needed to correct our @xmath0 measurements , @xmath96 , we simply computed the dot product between these two vectors and divided by the product of their lengths , @xmath97 .    finally , while computing the above geometric corrections , we also computed estimates of the apparent motion of cyg a within the coordinate system of each time step .
this was done to estimate the degree of doppler shifting of the temporal / spatial frequencies of any detected wave phenomena .
we did this for each time step by recomputing the position of the array center pierce - point for the two nearest time steps within the coordinate system of the current time step .
these positions were then used to numerically compute the time derivatives of the @xmath16 and @xmath61 coordinates of the array center pierce - point to obtain the north - south and east - west components of the sidereal velocity .
these are plotted in fig .
[ svel ] as functions of time along with a histogram for the azimuth angle ( measured north through east ) of the sidereal velocity vector for the entire observing run .
one can see from this figure that the velocities were sometimes significant , especially when the source was at lower elevations .
in addition , while the motion is generally from east to west , as one would naively assume , there is a significant spread in position angle of more than @xmath98 .
the authors would like to thank the referees for useful comments and suggestions .
basic research in astronomy at the naval research laboratory is supported by 6.1 base funding .
the vla was operated by the national radio astronomy observatory which is a facility of the national science foundation operated under cooperative agreement by associated universities , inc .
part of this research was carried out at the jet propulsion laboratory , california institute of technology , under a contract with the national aeronautics and space administration .
coker , c. , s.  e.  thonnard , k.  f.  dymond , t.  j.  w.  lazio , j.  j.  makela , and p. j. loughmiller ( 2009 ) , simultaneous radio interferometer and optical observations of ionospheric structure at the very large array _ rad .
, rs0a11    cornwell , t.  and e.  b.  fomalont ( 1999 ) , self - calibration , in _ synthesis imaging in radio astronomy ii , asp conference ser .
180 , edited by g.  b.  taylor , c.  l.  carilli , and r.  a.  perley , pp . 187199 , asp , san francisco , calif .
helmboldt , j.  f. , t.  j.  w.  lazio , h.  intema , and k.  dymond ( 2011 ) , a new technique for spectral analysis of ionospheric tec fluctuations observed with the very large array vhf system : from qp echoes to mstids _ rad .
, _ in press                    thompson , a.  r.  ( 1999 ) , fundamentals of radio interferometry , in _ synthesis imaging in radio astronomy ii , asp conference ser .
_ , vol . 180 , edited by g.  b.  taylor , c.  l.  carilli , and r.  a.  perley , pp . 187199 , asp , san francisco , calif .         for antenna 4 , rr polarization , as a function of time ( relative to midnight , 13 august ) when no attempt was made to `` de - spike '' the phase data ( see  2.2 )
middle panel : the difference between the cosine of the wrapped phase at a time step i , @xmath99 , and the next time step as a function of time used to find spikes in the phase data .
flagged spikes are highlighted with circles .
lower panel : the phase as a function of time , unwrapped with the flagged spikes excluded.,width=576 ]    , and the 327 mhz phase , @xmath55 , with the 327 mhz phase scaled such that the ionospheric components of the two phases canceled out for the rr ( upper ) and ll ( lower ) polarizations for antenna 14 .
the region where the instrumental phases have changed relatively abruptly ( i.e , the `` phase jump '' region first referred to in  2.2 ) is shaded in gray.,width=576 ]    , and the 327 mhz phase , @xmath55 , with the 327 mhz phase scaled such that the ionospheric components of the two phases canceled out for antenna 3 , ll polarization .
a fit to this data of a linear combination of the spatial frequencies @xmath2 and @xmath3 [ see equation ( 1 ) ] which assumes a single position offset ( see  2.2 ) is also plotted in red .
middle : the same as in the upper panel , but for the phases with the smoothed ( with a one - hour wide box ) phases subtracted from each band .
lower : the same as the middle panel , but with a fit of a linear combination of @xmath2 and @xmath3 used in the continuum subtraction of the 327 mhz phases at the edges of each scan ( see the discussion in
2.2).,width=576 ]    ) of the vla , the the difference between the tec fluctuation ( i.e. , above or below the mean tec ) along its line of sight and that measured along the reference antenna s line of sight , @xmath100 . the estimated uncertainty ( see  2.2 )
is plotted in red in each panel.,width=576 ]    ) of the vla , the the difference between the tec fluctuation ( i.e. , above or below the mean tec ) along its line of sight and that measured along the reference antenna s line of sight , @xmath100 . the estimated uncertainty ( see  2.2 ) is plotted in red in each panel.,width=576 ]    ) of the vla , the the difference between the tec fluctuation ( i.e. , above or below the mean tec ) along its line of sight and that measured along the reference antenna s line of sight , @xmath100 . the estimated uncertainty ( see  2.2 ) is plotted in red in each panel.,width=576 ]          ) at each antenna as a function of local time ( black curves ) .
plotted in red are the projected tec gradients computed using the fitted polynomial coefficients plotted in fig .
[ pfit ] . in each panel ,
the antenna number and the mean separation among the nearest antennas ( @xmath101 ) are printed.,width=576 ]