Ion and Electron Characteristics in Auroral Density Cavities Associated with Ion Beams: No Evidence for Cold Ionospheric Plasma
Submitted to JGR, May 15, 1998

J. P. McFadden, C. W. Carlson, and R. E. Ergun

Space Sciences Laboratory, University of California, Berkeley, CA.
 
D. M. Klumpar
Lockheed Martin Palo Alto Research Laboratory, Palo Alto, CA
 
E. Moebius
University of New Hampshire, Durham, NH
 

Abstract

Low density cavities associated with upgoing ion beams were identified in the auroral regions over two decades ago. In order to understand the waves, double layers, and solitary structures observed within these cavities, accurate measurements of the plasma distribution function are required. Although measurements by DE-1 indicated that these cavities were composed primarily of hot plasma in the form of ion beams, plasma sheet ions, and inverted-V electrons, later reports from Viking showed these cavities contained a cold plasma component whose density was an order of magnitude larger than the hot component. Recent measurements by the FAST satellite contrast sharply with the Viking results and support the earlier DE-1 observations. Regions of upgoing ion beams observed by FAST are shown to contain little or no cold plasma. The hot electron densities (>100 eV) and the combined plasma sheet ion and upgoing ion beam densities (>30 eV) agree remarkably well. Furthermore, no cold (0-30 eV) ions are measured at low energies by the mass spectrometer which precludes the presence of significant (>20%) cold electrons to preserve charge neutrality. Characteristics of the plasma for 11 ion beam events are tabulated.
 
Introduction
The auroral upward current regions contain low density cavities associated with upgoing ion beams (Persoon et al., 1988; Strangeway et al., 1998). These cavities are thought to be caused by the evacuation of plasma by field aligned potential drops at lower altitudes (Mozer et al., 1977). Accurate measurement of the particle distributions are required in order to characterize the waves, double layers and solitary structures observed in these density cavities (Kintner et al., 1978; Temerin et al., 1982; Koskinen et al., 1990). Accurate particle measurements are also necessary to understand how the plasma is able to support the parallel electric fields. This paper will identify the primary particle distributions observed in these ion beam regions or density cavities. Throughout this paper the term density cavity and ion beam region can be considered synonymous, however we primarily use the latter since historically these regions have been identified by the ion sensors.
In an extensive study using both wave and particle data from DE-1, Persoon et al. (1988) concluded that the thermal (<18 eV) electrons and ions were minor components of the plasma in density cavities associated with upgoing ion beams. Later observations from the Viking spacecraft (Koskinen et al., 1990) indicated that ion beam regions contained a cold (few eV) ion population that was approximately an order of magnitude denser than the energetic ion beam. Observations near their instrument sensitivity threshold also indicated that the cold ions were drifting up the field line. No mass separation was performed and the cold ions were assumed to be hydrogen. A cold (few eV) electron component with roughly the same density was inferred from Langmuir probe measurements (Bostrom et al., 1988; Koskinen et al., 1989) that could provide charge neutrality. These Viking observations were made at 8000-11000 km altitudes, in the afternoon-evening sector, and in aurora with a total potential drop of about 1 kilovolt.
During the first year of the FAST mission, numerous ion beams were identified at altitudes >2500 km, typically in the nightside winter hemisphere. The FAST observations described below were selected from a set of orbits near midnight, in the northern winter hemisphere, at altitudes near 4000 km. Below we present detailed observations from three ion beam regions and general properties of the 11 events studied. Selection criteria for the events included ion beams lasting >10 seconds, valid mass spectrometer data, and low amplitude (<100 mV/m) electric fields. We will show that the hot ion and electron densities agree within measurement uncertainties (~20%), and that no evidence for cold plasma is observed. These results support the observations reported by Persoon et al. (1988) and contrast sharply with Viking results.
 
Instrumentation
The FAST satellite contains top hat design, (Carlson et al., 1983; Carlson and McFadden., 1998a) electrostatic analyzers to measure the ion (3 eV - 25 keV) and electron (4 eV - 30 keV) distributions. The analyzers employ microchannel plate detectors in chevron configuration and biased to produce narrow pulse height distributions well above the preamplifier thresholds. The spacecraft is in a reverse cartwheel orbit and oriented relative to the magnetic field such that plasma analyzers measure full pitch angle distributions continuously. Analyzers include deflectors at their entrance aperture that steer their field of view by +/-10o to always include both the parallel and anti-parallel magnetic field line. Both analyzers form a 48 energy x 32 angle measurement each 78 ms. This study uses "survey data" formed by onboard averaging software which provides continuous coverage across the auroral regions. The electrostatic analyzers provide the bulk of the measurements used in this study, so we include a discussion of instrument calibration and measurement uncertainties below. A discussion of false counts in ion sensors (background, UV contamination, scattered electrons) is included in the "Observations" section where these sources of error are obvious from the data.
Large electric fields along the spacecraft velocity vector (approximately north-south) can produce ExB drifts that move field aligned ion beams out of the analyzer's 360ox6o planar (or conical when deflectors are used) field of view. This deflection introduces an error in the density calculation of ion beams. The deflection is most pronounced for oxygen beams, where a 100 mV/m field can deflect a 1 keV oxygen beam by about 5o for magnetic fields at 4000 km. Although the deflection of hydrogen beams is 1/4 that of oxygen beams at the same energy, hydrogen beams tend to be narrower (~10o FWHM) in pitch angle than oxygen beams (~30o FWHM) so the effect on the pitch angle distribution measurement is similar. A 5o deflection of a 1 keV proton beam requires a 400mV/m field at 4000 km altitudes, and reduces the calculated ion density to about 0.7 of the actual density for a 10o FWHM beam.
East-west electric fields also affect the ion density calculation by shifting the ion beam off the magnetic field line. This error increases the calculated beam density since a larger phase space volume is assumed in the density calculation. Although east-west electric fields tend to be smaller, errors caused by these fields can be larger. For the same 5o deflection of a 1 keV proton beam described above, the calculated ion beam density is about 1.6 times the actual density.
To reduce errors in the ion density calculation, slowly varying ion beam events with small electric fields (<100 mV/m during most of the ion beam) were selected for the study. Errors introduced by these fields are believed to be the major source of the small ~20% variations in the ratio of ion to electron densities seen below. It may be possible to calculate the ion density by correcting for large (>>100 mV/m) electric fields seen in the more turbulent ion beams, however this type of modeling is beyond the scope of this paper.
The ion and electron analyzer absolute geometric factors were estimated pre-flight from computer simulation, known grid transparencies, and published values of microchannel plate efficiency. The sensors were checked for uniformity using a Ni63 beta source and found to vary <20% around the analyzer field of view. The electron sensor was also tested in-flight for uniformity and small corrections in sensitivity with look direction are included in the calibration. The absolute geometric factor was tested in-flight by comparing the measured magnetic deflection across inverted-V arcs with the expected deflection determined from the measured electron flux and found to agree nearly exactly. (For examples see McFadden et al., 1998; Carlson et al., 1998b.)
The ion sensor calibration was tested by comparing the hot (>100 eV) electron and hot (>30 eV) ion densities in a selected ion beam (see Figure 1 and orbit 1616 in Table 1). A 35% adjustment to the pre-flight estimated sensitivity was required to increase the calculated ion density to that of the measured hot electron density. If cold electrons were present, a larger correction factor would be required. We attribute this correction to uncertainties in our knowledge of the absolute MCP efficiency for ions. Data from orbit 1616 were chosen because the plasma was primarily hydrogen (mass effects were nearly negligible on the density calculation) and because low level background counts due to scattered electrons in the ion sensor (see Observations section) could be ignored. In addition, no cold (<30 eV) ions were observed in the mass spectrometer, so the primary error in this normalization comes from a lack of knowledge of the cold (<100 eV) electron density (see Appendix: Cold Electron Measurements). The use of orbit 1616 data to normalize the ion sensor geometric factor would appear to invalidate this orbit as test for agreement between ion and electron densities. However, all 10 additional orbits showed good agreement between the hot ion and electron densities and a similar correction to the ion spectrometer's geometric factor would have been calculated for any of the 11 orbits. Since it is unlikely that the ratio of cold to hot electrons was the same for all 11 ion beam events, we have confidence that our normalization method is sound.
Finally, an additional source of error can arise from microchannel plate efficiency dependence upon energy and mass (Burrous et al., 1967; Gao et al., 1984). The FAST sensors include post-analyzer accelerations of ~450 V for electrons and ~1800 V for ions which should reduce these efficiency variations to the <10% level.

The FAST satellite also contains a Time-of-flight Energy-Angle imaging Mass Spectrometer (TEAMS) which measures the 3-D ion distribution each spin (Moebius et al., 1998a). These data, typically available with 5 second resolution, are used to obtain the relative H+, He+, and O+ densities across the ion beams. Pre-flight TEAMS calibrations are used for these data with estimated errors in relative sensitivity versus mass of <20% and with absolute sensitivity knowledge of ~50%. Since our calculations only use relative density versus mass to correct the upgoing ion beam density for mass, and since these corrections are typically a 25% density correction, TEAMS calibration errors should introduce <5% errors.

The relative ion densities vary slowly across arcs (Moebius et al., 1998b) and average ion density ratios were used to correct the higher time resolution ion spectrometer data for mass in the density calculation. The density ratio variations were small for most of the selected events, and we estimate <10% error introduced by changes in composition across these ion beams. However, two of the ion beams (orbits 1761 and 1793, Table 1) showed a larger variation in composition, which resulted in ~20% errors in the calculation during portions of these beams. A more sophisticated composition correction would have reduced this error, but was deemed unnecessary to improve the conclusions. The mass spectrometer data are also used to show that low level background counts in the ion spectrometer are due to internally scattered energetic electrons.

Observations

Figure 1 shows a section of an inverted-V arc (orbit 1616) containing ion beams that stretch over ~0.5o of invariant latitude. The top four panels show electron energy, electron pitch angle, ion energy and ion pitch angle spectrograms, respectively. The electrons have a spectral peak varying slowly between 3 and 1 keV and are isotropic in pitch angle except for the loss cone. Dropouts in the atmospheric secondary electrons (<100 eV, panel a) can be seen coincident with the appearance of the ion beam (panels c and d). The ion beam varies slowly in energy between 30 eV and 1 keV, and is collimated to within about 10o of the magnetic field. Panel e shows the hot (>100 eV) electron density calculated from the measured distributions. The low energy (<100 eV) electron measurement is dominated by spacecraft photoelectrons and scattered secondary electrons precluding any reliable measurement of this density (see Appendix: Cold Electron Measurements). Panels f and g show the hot ion beam (>30 eV, 150o-180o pitch angle) and plasma sheet ion (>30 eV, 0o-150o pitch angle) densities calculated from the measured distributions. As discussed below, the mass spectrometer shows that no cold (<30 eV) ions are present in the ion beam region. The density calculation requires a mass assumption (n ~ (count rate) x  ) since the ion spectrometer is not mass resolving. The plasma sheet ion density is calculated assuming hydrogen, which is consistent with the mass spectrometer data. The ion beam density has been calculated using a mass,  , given by:
where np, nh, no are the relative densities of hydrogen, helium, and oxygen, determined by the mass spectrometer and averaged over the ion beam. For this event, the above correction is  . Panel h shows the ratio of hot ion to hot electron densities which shows remarkably good agreement between the densities inside the ion beam even though the beam density (panel f) varies by a factor of ~5 and the electron density (panel f) varies by a factor of ~2.
Now consider the measured ion distributions and why we introduced a cutoff of 30 eV in the ion density calculation. Figure 2a shows a contour plot of the ion differential energy flux (proportional to sensor count rate) for a 10 second average during the ion beam in Figure 1. The beam at 180o pitch angle, and the hot, nearly isotropic, plasma sheet ions are easily identified. The average counts per energy-angle bin below 100 eV is about 1 in the contour plot. The slight enhancement in low energy counts at pitch angles between 180o and 270o is due to scattered ultraviolet into the sensor. Figure 2b and 2c show the differential flux and count rate for this event, summed over pitch angle. Note that the count rate (2c) is roughly constant, or independent of energy, below 50 eV and forms a differential flux spectra (2b) proportional to 1/energy. Figure 2d shows the partial ion density at each energy step assuming H+ (blue line, , where Ei is the center energy of each energy step and  is the energy width), and the partial integral of the density starting at the highest step (green line,  , where Emax is the highest energy step). One can see that most of the density resides in the hot component, however as much as 20% appears to exist in the cold ions (4-30 eV). These low energy counts are about an order of magnitude above the microchannel plate noise levels. However, as shown below, these counts are due to a small number of scattered energetic electrons and a small amount of solar ultraviolet and do not represent actual cold plasma.
In addition to the ion spectrometer, the FAST satellite has a mass spectrometer which covers roughly the same energy range. The overall geometric factor of the mass spectrometer is about one half that of the ion spectrometer, producing count rates that differ by a factor of ~2 between the instruments. Figure 3a shows the count rate, summed over angle, from the mass spectrometer for the same interval as Figure 2. Total counts from the ion beam are smaller by a factor of ~50 for the mass spectrometer since it samples the ion beam only twice per spin whereas the ion spectrometer samples the beam continuously. A better comparison uses the nearly isotropic plasma sheet ion count rates (>5 keV) which differ by ~2 as expected for the two instruments.

In contrast to the ion spectrometer, the mass spectrometer measures no counts below about 40 eV. This is because the mass spectrometer is not susceptible to scattered electrons (or sunlight) due to coincidence detection in the time-of-flight electronics and due to a large (-20 kV) post acceleration between the analyzer and time-of-flight detector (Moebius et al., 1998a). Examination of ion beam events has shown that the ion spectrometer's low energy count rate increases with electron flux and electron energy, whereas the mass spectrometer observes no counts at these energies. For the event in Figure 1, false counts due to scattered electrons are rather small because most of the inverted-V electrons are repelled by the -1800 V bias on the front of the ion spectrometer's microchannel plate. In energetic arcs (~10 keV), the scattered electron count rate can be an order of magnitude larger. For these reasons, density calculations using the ion spectrometer do not include the low energy (<30 eV) counts due to scattered electrons, which is consistent with no counts in the mass spectrometer data at these energies.

Figure 3b uses the same format as Figure 2d to display the partial ion density at each energy step assuming H+ (blue line), and the partial integral of the density starting at the highest step (green line) for the TEAMS measurement. Here we have summed the counts at all masses and assumed H+ to mimic the ion spectrometer measurement, but have used the full 3-D TEAMS measurement. A comparison of the two density calculations show quite good agreement (only 20% difference when the <30 eV counts are ignored) between instruments with different sensitivities, field-of-views, and detection systems.

Before proceeding with additional examples, a final word should be said about the possible presence of cold ions whose energy is below the mass spectrometer's lower energy threshold of 2-5 eV (which depends upon spacecraft potential), and the possibility of cold electrons below our 100 eV cutoff. First, it would be quite remarkable if this plasma existed, since it would require both the cold electron and cold ion populations to be present with the same density, independent of hot plasma densities. However, the proof comes from the presence of large (~1 V/m) perpendicular electric fields in the stronger ion beams not used in this study. The ExB drifts from these fields would cause cold hydrogen (oxygen) to have energies of tens (hundreds) of eV in the spacecraft frame. Low energy cold ions moving perpendicular to B would be easily detected by FAST and have not been observed during these strong electric fields. The lack of cold ions precludes the presence of significant cold electrons to preserve charge neutrality. In addition, Strangeway et al. (1998) has shown that the hot electron density is consistent with measured plasma wave cutoffs, and Ergun et al. (1998) has shown that the AKR spectrum is consistent with a hot, relativistic electron population, providing further proof that the cold plasma is a small fraction of the total density.

Figure 4 shows an example of an ion beam event from orbit 1676 with the same format as Figure 1. The top 2 panels show an isotropic inverted-V electron distribution whose energy is just above 1 keV. The ion beam (panels c and d) varies between 100 eV and 2 keV. Low energy ion conics are also observed at 10:05:41-45 and after 10:05:49. Panel e shows the hot (>100 eV) electron density and panels f and g show the ion beam (>30 eV, 150o-180o pitch angle) and plasma sheet ion (>30 eV, 0o-150o pitch angle) densities. The ion beam density calculation used  as determined from the mass spectrometer data. The jump in plasma sheet ion density (panel g) at 10:05:41-:45 and after 10:05:49 is caused by ion conics included in the integral. Panel h shows that the ratio of hot ion to hot electron density is about 1 when the ion beam is present.

The smaller flux of inverted-V electrons in Figure 4a, as compared to Figure 1a, resulted in a factor of 3 reduction in scattered electrons into the ion sensor. Figure 5 shows a 10 second average of the ions in a format identical to Figure 2. Figure 5c shows a relatively constant count rate below 50 eV, similar to that in Figure 2c but with a lower rate. Mass spectrometer data for this beam shows no ions below 40 eV indicating these counts are instrumental. As can be seen from Figure 5d, the counts resulting from scattered electrons would have resulted in only a ~20% error for the ion density.

Figure 6 shows a third example of an ion beam event, this time from orbit 1771. The format is the same as Figures 1 and 4. The top 4 panels show an isotropic inverted-V electron population with an energy of ~1 keV, and an ion beam varying from about 100 eV to 1 keV. The mass spectrometer gave a beam mass of  for the beam density calculation (panel f). The plasma sheet ion density (panel g) is almost negligible (~0.05 cm-3) for this event. Although the ion beam density (panel f) varies by nearly a factor of 8, panel h shows that the ratio of hot ion to hot electron densities remains near one across the event.

A comparison of ion beams on the three orbits (Figures 1, 4, and 6) shows a number of differences and similarities between the ion beams. The density was highest for orbit 1616, varying from about 1 to 3 cm-3, whereas the beam density on orbit 1676 remained below 0.6 cm-3. Orbit 1771 showed the largest variations in density swing from 0.2 to 1.1 cm-3. Slightly more oxygen was present during orbits 1676 and 1771, as seen from the beam mass correction or Table 1. The additional oxygen also shows up in contour plot 5a as the wider pitch angle portion of the beam at higher energies. In all three events the plasma sheet ions have a nearly constant density across the event, whereas the ion beam shows large density variations. To first order the beam density variations reflect conservation of flux, where a relatively constant flux of upgoing ions into the acceleration region should show a density inversely proportional to the  .

Table 1 summarizes characteristics of the electrons and ions for 11 ion beam events investigated in detail. All these events give good agreement between the hot ion and hot electron densities. Orbit 1623 was the lowest energy event and the lower limit for the ion beam calculation was dropped to 15 eV for this event. Most of the ion beam events had densities less than 1 cm-3, and the lowest observed density was 0.13 cm-3. Variations in the O+/H+ and He+/H+ beam density ratios were significant (greater than a factor of 2) for many of the events, especially orbit 1761 where a factor of 5 variation in O+/H+ was observed across the ion beam. However, the mass correction to the beam density for this event,  , varied by only +21% to -22% across the event due to additional changes in He+/H+, and was relatively constant across most of the event giving good agreement between ion and electron densities. A more sophisticated program that uses a time varying,  , across the beams could have produced limited improvement in a couple of the density ratio calculations but was unnecessary for this study. Table 1 provides theorists a range of parameters that can be used to model the ion beam regions.

Discussion

The results presented above show convincing evidence that the hot electron and hot ion densities in auroral ion beam regions provide the bulk of the total density. The measurements were made by high resolution electron and ion spectrometers, with mass corrections provided by a time-of-flight mass spectrometer. The primary errors in this measurement appear to be introduced by large, low frequency electric fields, or electrostatic shocks, that deflect the ion beam transverse to the field line due to ExB drifts. These deflections result in <20% variations in the calculated density ratio on short time scales, however the density ratio averaged over the events agree to better than 10%. The cold ion (<30 eV) measurement is difficult with a spectrometer because a small fraction of energetic electrons can scatter through the analyzer producing false counts that corrupt the low energy measurement. This contamination is not present in the mass spectrometer data, and these data were used to show that no measurable ion fluxes were present below 30 eV. The detection of small fluxes of cold (<100 eV) electrons is very difficult to make with a spectrometer due to spacecraft photoelectrons and secondary electrons produced in the analyzer aperture (see Appendix: Cold Electron Measurements), so only the hot electron density was calculated. However, to preserve charge neutrality, we infer that no significant cold electron population is present during ion beams. This result is consistent with ion beam regions being produced by quasi-static parallel electric fields, which prevent cold electrons from penetrating into the density cavities associated with ion beams. In addition, wave observations by Strangeway et al. (1998) and Ergun et al. (1998) have shown that plasma wave cutoffs are also consistent with a hot (>100 eV) electron population providing the bulk of the total density.
The evidence for a lack of cold plasma in the density cavity supports the DE-1 results presented by Persoon et al. (1988), but differs from previous reports of Viking data. Koskinen et al. (1990) report evidence of a low energy, cold ion component moving up the field line during ion beam events. The density of the cold ion component was estimated to be an order of magnitude larger than the hot ion beam. Figure 6 of the Koskinen et al. (1990) paper shows ion spectra that resulted in the identification of cold plasma. The spectra below 100 eV varies as 1/energy, which implies a constant count rate similar to Figures 2c and 5c of this paper. Koskinen et al. (1990) conclude that the measurement is valid since the low energy counts are significantly above their instrument sensitivity. However, no discussion of possible sources of false counts is included in the paper, specifically the instrument's sensitivity to scattered energetic electrons. The Viking measurement of the low energy ion flux is an order of magnitude larger than the false count rate in the FAST ion spectrometer, and 2 orders of magnitude above the sensitivity of the FAST mass spectrometer which saw no flux at these energies.
To provide charge neutrality, Koskinen et al. (1990) report that a cold electron component was also present during the ion beams with about the same density as the cold ion component. This result was based upon Langmuir probe sweeps near the ion beams which inferred a 1-10 eV temperature for the electrons (Bostrom et al., 1988). However, Hilgers et al. (1992) shows that the Viking Langmuir probe current is dominated by photoelectrons from the wire boom within ion beam regions. They conclude that only an upper limit to the thermal electron density can be determined and that cold electrons are at most of the order of the energetic electron density.
It may be possible that differences between the FAST and Viking observations are due to other effects. Both data sets were measured at similar local times, with Viking observations in the afternoon-evening sector and FAST observations at ~22 MLT. Viking measurements were made at higher altitudes (8000-11000 km), typically above most of the acceleration region, whereas FAST observations (~4000 km) are near the bottom of the acceleration region. However, there are no reported observations or theories that suggest the plasma density should increase with altitude inside the ion beam regions, with the lowest densities at the bottom of the density cavities. The Viking observations were made near solar minimum, as were the FAST observations. Koskinen et al. (1990) state that the events were from "relatively weak acceleration events associated with low geomagnetic activity". As seen from Table 1, most of the FAST events were also during low activity with relatively weak acceleration (~1-2 keV ion beams). In both cases the O+ beams tended to have a higher energy and broader pitch angle width than the H+ beams. Koskinen et al. (1990) report a relatively small (<30%) heavy ion density in the upflowing ion beams. As seen from Table 1, the FAST observations differ in that the sum of the He+ and O+ densities typically exceeded the H+ density for the 11 events. These are typical of ion beams reported by Moebius et al. (1998b) in a statistical study of FAST ion beam events. The Viking events were also selected for a study of solitary waves. The 11 FAST events were examined and solitary wave structures were found in all the beams. From the above analysis and comparison, which suggests measurements of similar phenomena under similar conditions, we conclude that the interpretation of Viking measurements as observations of cold plasma in ion beam regions are in error.
The new FAST observations, which find only hot plasma in the ion beam regions, have strong implications for the density profile and plasma waves in the upward current regions. Without cold plasma, certain wave modes will disappear and other modes grow without heavy damping by the cold electrons. Studies of wave modes and solitary structures within ion beam regions that assumed a cold background plasma (Marchenko and Hudson, 1995, and references therein) need to be revisited in light of these new measurements. The need to postulate anomalous resistivity to prevent current runaway may also be unnecessary since the parallel acceleration region appears to contain only hot plasma. The measurements are consistent with simple acceleration by parallel electric fields which prevent cold ionospheric and secondary electrons from penetrating into the acceleration region, and accelerate all the cold ions (or cool ion conics) into energetic ion beams. The density profile along the magnetic field line appears to be a simple function of the parallel electric field profile, the high altitude plasma distribution, and the ionospheric plasma distributions. It may be possible to determine the electric field profile directly from the low and high altitude plasma distributions, combined with charge neutrality and current continuity in a region carrying an upward current. It is hoped that these new measurements will allow a consistent theory to be developed to explain the parallel electric fields and waves in the upward current regions.
 
Appendix: Cold Electron Measurements
Low energy electron measurements in the auroral zone are extremely difficult due to spacecraft photoelectron and secondary electron contamination. This appendix is included to demonstrate the difficultly in this measurement, to identify the dominant sources of low energy electron counts, and to justify the 100 eV cutoff that was used. Figure 7 shows contour plots of the electron differential energy flux (proportional to count rate) both inside (7a) and outside (7b) the ion beam of Figure 1. There are 6 primary sources of counts below 100 eV: background microchannel plate noise, solar ultraviolet scattering into the sensor, spacecraft generated photoelectrons, internally scattered energetic electrons from the auroral electron beam, trapped plasma electrons outside the loss/antiloss cone, and atmospheric secondary electrons in the loss/antiloss cone. Only the last two sources of counts represent plasma that should be used in the density calculation. Microchannel plate background noise caused by radioactive decay and penetrating radiation are negligible. Orientation of the analyzer and blackened, scalloped analyzer hemispheres reduce solar ultraviolet to negligible levels also. Spacecraft photoelectrons can be seen in both distributions near 90o and 270o pitch angles. These electrons disappear when the spacecraft is in shadow, but dominate the low energy measurement in sunlight. Scattered electrons from the primary auroral electron beam dominate the count rate during ion beams (7a) at pitch angles near 0o and 180o. Atmospheric secondary and trapped electrons dominate the count rate when the spacecraft is below the ion beam regions (7b), except at the lowest energies and near 90o pitch angles where spacecraft photoelectrons always dominate.
From the above example, one can see that a quantitative low energy electron measurement is difficult at best. Since low energy counts can dominate the density calculation due to weighting by 1/ in the density calculation, a spectrometer determined density measurement that includes low energy electrons requires close scrutiny. Figure 7c, formed from the distribution in 7a within the ion beam, shows the partial electron density at each energy step (blue line, , where Ei is the center energy of each energy step and  is the energy width), and the partial integral of the density starting at the highest step (green line,  , where Emax is the highest energy step). The contribution of photoelectrons and internally scattered secondary electrons to the density clearly becomes significant below 50 eV. The plateau of the green curve between 50 and 300 eV suggests minimal contribution of real electrons to the total density at these low energies. This can be contrasted with Figure 7d, formed from the distribution in 7b, which shows that atmospheric secondary electrons between 50 and 300 eV contribute significantly to the total local density outside ion beams.
Although we neglect the low energy (<100 eV) electrons in the density calculation due to instrumental reasons, this assumption can be justified if the ion beams are formed by quasi-static parallel electric fields. In this case cold plasma should be excluded from the ion beam regions. One can also estimate the density of neglected atmospheric secondary and trapped electrons with energy <100 eV within an ion beam by folding an electron distribution measured outside the ion beam through a 1 kilovolt potential. The potential eliminates portions of the initial distribution below 1 keV, and adiabatically maps the measured distribution between 1 and 1.1 keV to the 0 to 100 eV energies. The calculated <100 eV density for this mapped distribution is only a few percent of the hot (>100 eV) electron density. Since the observed hot electron density agrees with the ion density as shown in Figures 1, 4, and 6, our assumption of negligible cold (<100 eV) plasma appears to be justified. We note that the assumption that the <100 eV electrons can be neglected is not valid when the ion beam has a low energy since significant secondary electrons can then penetrate the potential below the spacecraft. This effect can be seen in Figure 1 at 20:53:34, where the ion beam density (panel f) goes up, and is not balanced by the hot electrons (panel e).
The analysis of FAST data was supported by NASA grant NAG5-3596.
 
References
Bostrom, R., G. Gustafsson, B. Holback, G. Holmgren, H. Koskinen, and P. Kintner, Characteristics of solitary waves and weak double layers in the magnetospheric plasma, Phys. Rev. Lett., 61, 82, 1988.
 
Burrous, C. N., A. J. Lieber, and V. T. Zavansteff, Detection efficiency of a continuous channel electron multiplier for positive ions, Rev. Sci. Instrum., 38, 1477, 1967.
 
Carlson, C. W., D. W. Curtis, G. Paschmann, and W. Michael, An instrument for rapidly measuring plasma distribution functions with high resolution, Adv. Space Res., 2, 67, 1983.
 
Carlson, C. W., and J. P. McFadden, Design and applications of imaging plasma instruments, AGU Monog. Meas. Techn. Space Plasmas, ed. R. Pfaff, in press, 1998a.
 
Carlson, C. W., J. P. McFadden, R. E. Ergun, M. Temerin, W. Peria, F. S. Mozer, D. M. Klumpar, E. G. Shelley, W. K. Peterson, E. Moebius, R. Elphic, R. Strangeway, C. Cattell, R. Pfaff, FAST observations in the downward auroral current region: energetic upgoing electron beams, parallel potential drops, and ion heating, Geophys. Res. Lett., in press, 1998b.
 
Ergun, R. E., C. W. Carlson, J. P. McFadden, F. S. Mozer, G. T. Delory, W. Peria, C. C. Chaston, M. Temerin, R. Elphic, R. Strangeway, R. Pfaff, C. A. Cattell, D. Klumpar, E. Shelley, W. Peterson, E. Moebius, L. Kistler, FAST satellite wave observations in the AKR source region, Geophys. Res. Lett., in press, 1998.
Gao, R. S., P. S. Gibner, J. H. Newman, K. A. Smith, and R. F. Stebbings, Absolute and angular efficiencies of a microchannel-plate position-sensitive detector, Rev. Sci. Instrum., 55, 1756, 1984.

Gurnett, D. A. and L. A. Frank, VLF hiss and related plasma observations in the polar magnetosphere, J. Geophys. Res., 77, 172, 1972.

Hilgers, A., B. Holback, G. Holmgren, and R. Bostrom, Probe measurements of low plasma densities with applications to the auroral acceleration region and auroral kilometric radiation sources, J. Geophys. Res., 97, 8631, 1992.

Kintner, P. M., M. C. Kelley, and F. S. Mozer, Electrostatic hydrogen cyclotron waves near one earth radius altitude in the polar magnetosphere, Geophys. Res. Lett., 5, 139, 1978.

Koskinen, H., R. Bostrom, and B. Holback, Viking observations of solitary waves and weak double layers on auroral field lines, in Ionosphere-Magnetosphere-Solar Wind Coupling Processes, edited by T. Chang, G. B. Crew, and J. R. Jasperse, p.147, Scientific, Cambridge, Mass., 1989.

Koskinen, H., R. Lundin, B. Holback, On the plasma environment of solitary waves and weak double layers, J. Geophys. Res., 95, 5921, 1990.

Marchenko, V. A., M. K. Hudson, Beam driven acoustic solitary waves in the auroral acceleration region, J. Geophys. Res., 100, 19791, 1995.

McFadden, J. P., C. W. Carlson, R. E. Ergun, C. C. Chaston, F. S. Mozer, M. Temerin, D. M. Klumpar, E. G. Shelley, W. K. Peterson, E. Moebius, L. Kistler, R. Elphic, R. Strangeway, C. Cattell, R. Pfaff., Electron modulation and ion cyclotron waves observed by FAST, Geophys. Res. Lett., in press, 1998.

Moebius, E., L. M. Kistler, M. A. Popecki, K. N. Crocker, M. Granoff, Y. Jiang, E. Sartori, V. Ye, H. Reme, J. A. Sauvaud, A. Cros, C. Aoustin, T. Camus, J.-L. Medale, J. Rouzaud, C. W. Carlson, J. P. McFadden, D. Curtis, H. Heetderks, J. Croyle, C. Ingraham, B. Klecker, D. Hovestadt, M. Ertl, F. Eberl, H. Kastle, E. Kunneth, P. Laeverenz, E. Seidenschwang, E. G. Shelley, D. M. Klumpar, E. Hertzberg, G. K. Parks, M. McCarthy, A. Korth, H. Rosenbauer, B. Grave, L. Eliasson, S. Olsen, H. Balsiger, U. Schwab, M. Steinacher, The 3-D plasma distribution function analyzers with time-of-flight mass discrimination for Cluster, FAST, and Equator-S, AGU Monog. Meas. Techn. Space Plasmas, ed. R. Pfaff, in press, 1998a.

Moebius, E. L. Tang, L. M. Kistler, M. Popecki, E. J. Lund, D. Klumpar, W. Peterson, E. G. Shelley, B. Klecker, D. Hovestadt, C. W. Carlson, R. Ergun, J. P. McFadden, F. Mozer, M. Temerin, C. Cattell, R. Elphic, R. Strangeway, R. Pfaff, Species dependent energies in upward directed ion beams over auroral arcs as observed with FAST TEAMS, Geophys. Res. Lett., in press, 1998b.

Mozer, F. S., C. W. Carlson, M. K. Hudson, R. B. Torbert, B. Parady, J. Yatteau, and M. C. Kelley, Observations of paired electrostatic shocks in the polar magnetosphere, Phys. Rev. Lett., 38, 292, 1977.

Persoon A. M., D. A. Gurnett, W. K. Peterson, J. H. Waite, Jr., J. L. Burch, and J. L. Green, Electron density depletions in the nightside auroral zone, J. Geophys. Res., 93, 1871, 1988.

Strangeway, R. J., L. Kepko, R. C. Elphic, C. W. Carlson, R. E. Ergun, J. P. McFadden, W. J. Peria, G. T. Delory, C. C. Chaston, M. Temerin, C. A. Cattell, E. Moebius, L. M. Kistler, D. M. Klumpar, W. K. Peterson, E. G. Shelley, R. F. Pfaff, FAST observations of VLF waves in the auroral zone: evidence of very low plasma densities, Geophys. Res. Lett., submitted, 1998.

J. P. McFadden, C. W. Carlson, R. E. Ergun, Space Sciences Laboratory, University of California, Berkeley, CA 94720-7450. (e-mail: mcfadden@ssl.berkeley.edu)

D. M. Klumpar, Lockheed Martin Palo Alto Research Laboratory, Palo Alto, CA 94304

E. Moebius, L. Kistler, U. of New Hampshire, Durham, NH 93824

(Received May?, 1998; accepted?)

Copyright 1998 by the American Geophysical Union.

Figure Captions

Figure 1.
The figure shows a comparison of the hot (>100 eV) electron and hot (>30 eV) ion densities across an ion beam. Panels a and b show electron energy and pitch angle spectrograms. Panels c and d show ion energy and pitch angle spectrograms. Panels e, f and g are the integrated hot electron density, ion beam density corrected for beam composition, and hot plasma sheet ion density. Panel h is the ratio of ion to electron density, which is nearly one across the ion beam.
Figure 2.
Figure 2a is a contour plot of the ion differential energy flux (proportional to count rate) averaged over the beam in Figure 1. The 10 second average has broadened the ion beam's parallel energy. Isotropic plasma sheet ions, minus the loss cone, are seen at 10 keV. Counts below 30 eV are primarily due to internally scattered inverted-V electrons, along with some scattered solar ultraviolet (between 180o and 270o). These false counts have a constant count rate below 30 eV (see Figure 2c) and produce a spectra (Figure 2b) proportional to 1/energy. Figure 2d shows the partial density (blue line, , where Ei is the center energy of each energy step and  is the energy width), and the partial integral of the density starting at the highest step (green line,  , where Emax is the highest energy step). From 2d it can be seen that the low energy false counts would produce a ~20% error in the total density if included.
Figure 3.
Figure 3a is the count rate from the FAST mass spectrometer (TEAMS) averaged over the same 10 second interval as Figure 2. No counts below 30 eV are measured by this instrument, proving that the constant count rate seen at low energies in Figure 2c is not real. The count rates for the nearly isotropic plasma sheet ions (>5 keV) differ by a factor of two between Figures 3a and 2c, as expected for the two instruments' relative geometric factors. The ion beam is not prominent in the mass spectrometer data since it samples the beam only twice per spin. The partial density calculation in Figure 3b is nearly identical to Figure 2d except for the error at low energy (<30 eV) caused by false counts in the spectrometer data.
Figure 4.

The figure shows a comparison of hot (>100 eV) electron and hot (>30 eV) ion densities for orbit 1676 with the same format as Figure 1. Panels e, f, and g show the electron, ion beam and plasma sheet ion densities respectively. Panel h shows the ratio of ion to electron density is near one during the ion beam.

Figure 5.

The figure shows ion characteristics during the ion beam from Figure 4. The format is the same as Figure 2 and displays many of the same characteristics. The ion beam and plasma sheet ions are easily seen in Figure 5a. A weaker flux of inverted-V electrons (see Figure 4) produces a smaller false count rate in Figure 5c at energies <30 eV than is seen in Figure 2c.

Figure 6.

The figure shows a comparison of hot (>100 eV) electron and hot (>30 eV) ion densities for orbit 1771 with the same format as Figures 1 and 4. Panels e, f, and g show the electron, ion beam and plasma sheet ion densities respectively. Although the total density (panel e) varies by a factor of ~7 across this beam event, the ratio of ion to electron density remains near one (panel h).

Figure 7.

Contour plots of the electron differential energy flux (proportional to count rate) both inside (Figure 7a) and outside (Figure 7b) the ion beam region of Figure 1. Figures 7c and 7d show the partial density (blue line,  , where Ei is the center energy of each energy step and  is the energy width), and the partial integral of the density starting at the highest step (green line,  , where Emax is the highest energy step) for the contour plots 7a and 7b, respectively. Photoelectrons (<100 eV near -90o, 90o, 270o pitch angles) dominate the density calculation if included in the sum. Atmospheric secondary electrons are seen in Figure 7b at energies <200 eV and angles -50o to 50o and 130o to 230o, and secondary electrons produced in the analyzer are seen at these angles at <50 eV in Figure 7a.
 

 

   Table 1

Orbit 
   #
Ratio 
Ni/Ne
KP Electron 
density 
Range 
(cm-3)
Peak Energy 
(kev)
Ion 
Beam 
Density 
(cm-3)
O+/H+ 
Ratio*
He+/H+ 
Ratio*
Energy 
range 
(keV)
PS 
Density 
(cm-3)
PS 
Temp 
(keV)
1525 1.13+/-0.11 1 .55-.85 1 .3-.6 0.28 0.67 .3-1. .35-.45 10
1616 0.99+/-0.08 1 1.2-2.5 1-2 .4-2.2 0.36 0.65 .1-1. .7-.8 5
1623 0.98+/-0.13 1- .9-2.5 .8 .45-2.3 0.54 0.34 .1-.3 .3 2
1633 1.00+/-0.14 1 .4-.7 2-3 .15-.6 1.28 0.47 .1-3 .15 5
1676 0.93+/-0.09 1- .33-.6 1.5 .15-.45 0.87 0.20 .1-2 .15-.2 3
1692 0.83+/-0.12 2+ .4-1.2 1-2 .15-.75 0.87 0.54 .3-2 .2-.35 6
1740 0.98+/-0.16 2- .13-.40 1-2 0.08-.2 1.09 0.38 1-2 0.04-.1 5
1761 0.98+/-0.16 1- .13-.28 2-3 .05-.1 1.44 0.46 1-5 0.06-.1 7
1771 0.88+/-0.13 1- .2-1.2 1-2 .15-.9 0.54 0.57 .3-2 .05-.15 3
1793** 1.04+/-0.14 1 .25-1.0 3-10 .15-6. 1.42 0.23 .3-20. .1-.25 10
1804 0.93+/-0.13 1 .2-1.0 1-10 .1-.5 1.11 0.31 .2-8 .05-.2 4

 

All ion beam events are between 21 and 1 MLT, at altitudes ~4000 km.

* Averages over beam event. Factor of 2 to 5 variations may occur across a particular beam event. The largest variation occurred across orbit 1761 with O+/H+ changing from 0.5 to 2.5.

** Only 2/3 of orbit 1793 beam event is used due to large electric fields.