Search

It is very well known that the effects of electromagnetic ion cyclotron (EMIC) waves on ring current (RC) ion and radiation belt (RB) electron dynamics strongly depend on such particle/wave characteristics as the phase-space distribution function, frequency, wave-normal angle, wave energy, and the form of wave spectral energy density. The consequence is that accurate modeling of EMIC waves and RC particles requires robust inclusion of the interdependent dynamics of wave growth/damping, wave propagation, and particles. Such a self-consistent model is being progressively developed by Khazanov et al. This model is based on a system of coupled kinetic equations for the RC and EMIC wave power spectral density along with the ray tracing equations. We will discuss the recent progress in understanding EMIC waves formation mechanisms in the inner magnetosphere. This problem remains unsettled in spite of many years of experimental and theoretical studies. Modern satellite observations by CRRES, Polar and Cluster still do not reveal the whole picture experimentally since they do not stay long enough in the generation region to give a full account of all the spatio-temporal structure of EMIC waves. The complete self-consistent theory taking into account all factors significant for EMIC waves generation remains to be developed. Several mechanisms are discussed with respect to formation of EMIC waves, among them are nonlinear modification of the ionospheric reflection by precipitating energetic protons, modulation of ion-cyclotron instability by long-period (Pc3/4) pulsations, reflection of waves from layers of heavy-ion gyroresonances, and nonlinearities of wave generation process. We show that each of these mechanisms have their attractive features and explains certain part experimental data but any of them, if taken alone, meets some difficulties when compared to observations. We conclude that development of a refined nonlinear theory and further correlated analysis of modern satellite and ground-based data is needed to solve this very intriguing problem.

It is well known that the effects of electromagnetic ion cyclotron (EMIC) waves on ring current (RC) ion and radiation belt (RB) electron dynamics strongly depend on such particle/wave characteristics as the phase-space distribution function, frequency, wave-normal angle, wave energy, and the form of wave spectral energy density. The consequence is that accurate modeling of EMIC waves and RC particles requires robust inclusion of the interdependent dynamics of wave growth/damping, wave propagation, and particles. Such a self-consistent model is being progressively developed by Khazanov et al. [2002 - 2007]. This model is based on a system of coupled kinetic equations for the RC and EMIC wave power spectral density along with the ray tracing equations. We will discuss the recent progress in understanding EMIC waves formation mechanisms in the inner magnetosphere. This problem remains unsettled in spite of many years of experimental and theoretical studies. Modern satellite observations by CRRES, Polar and Cluster still do not reveal the whole picture experimentally since they do not stay long enough in the generation region to give a full account of all the spatio-temporal structure of EMIC waves. The complete self-consistent theory taking into account all factors significant for EMIC waves generation remains to be developed. Several mechanisms are discussed with respect to formation of EMIC waves, among them are nonlinear modification of the ionospheric reflection by precipitating energetic protons, modulation of ion-cyclotron instability by long-period (Pc3/4) pulsations, reflection of waves from layers of heavy-ion gyroresonances, and nonlinearities of wave generation process. We show that each of these mechanisms have their attractive features and explains certain part experimental data but any of them, if taken alone, meets some difficulties when compared to observations. We conclude that development of a refined nonlinear theory and further correlated analysis of modern satellite and ground-based data is needed to solve this very intriguing problem.

This paper continues presentation and discussion of the results from our new global self-consistent theoretical model of interacting ring current ions and propagating electromagnetic ion cyclotron waves [Khazanov et al., 2006]. To study the effects of electromagnetic ion cyclotron wave propagation and refraction on the wave induced ring current precipitation and heating of the thermal plasmaspheric electrons, we simulate the May 1998 storm. The main findings after a simulation can be summarized as follows. Firstly, the wave induced ring current precipitation exhibits quite a lot of fine structure, and is highly organized by location of the plasmapause gradient. The strongest fluxes of about 4 x 10(exp 6) (cm(raised dot) s(raised dot) sr(raised dot) (sup -1)) are observed during the maill and early recovery phases of the storm. The very interesting and probably more important finding is that in a number of cases the most intense precipitating fluxes are not connected to the most intense waves in simple manner. The characteristics of the wave power spectral density distribution over the wave normal angle are extremely crucial for the effectiveness of the ring current ion scattering. Secondly, comparison of the global proton precipitating patterns with the results from RAM [Kozyra et al., 1997a] reveals that although we observe a qualitative agreement between the localizations of the wave induced precipitations in the models, there is no quantitative agreement between the magnitudes of the fluxes. The quantitative differences are mainly due to a qualitative difference between the characteristics of the wave power spectral density distributions over the wave normal angle in RAM and in our model. Thirdly, the heat fluxes to plasmaspheric electrons caused by Landau resonate energy absorption from electromagnetic ion cyclotron waves are observed in the postnoon-premidnight MLT sector, and can reach the magnitude of 10(exp 11) eV/(cm(sup 2)(raised dot)s). The Coulomb energy degradation of the RC H(+) and O(+) ions maximizes at about 10(exp 11) (eV/(cm(sup 2) (raised dot) s), and typically leads to electron energy deposition rates of about 2(raised dot) 10(exp 10) (eV/(cm(sup 2)(raised dot)s) which are observed during two periods; 32-48 hours, and 76-86 hours after 1 May, 0000 UT. The theoretically derived spatial structure of the thermal electron heating caused by interaction of the ring current with the plasmasphere is strongly supported by concurrent and conjugate plasma measurements from the plasmasphere, ring current, and topside ionosphere [Gurgiolo et al., 2005]. Finally, the wave induced intense electron heating has a structure of the spot-like patches along the most enhanced density gradients in the plasmasphere boundary layer and can be a possible driver to the observed but still not explained small-scale structures of enhanced emissions in the stable auroral red arcs.

The further development of a self-consistent theoretical model of interacting ring current ions and electromagnetic ion cyclotron waves (Khazanov et al., 2003) is presented In order to adequately take into account wave propagation and refraction in a multi-ion magnetosphere, we explicitly include the ray tracing equations in our previous self-consistent model and use the general form of the wave kinetic equation. This is a major new feature of the present model and, to the best of our knowledge, the ray tracing equations for the first time are explicitly employed on a global magnetospheric scale in order to self-consistently simulate the spatial, temporal, and spectral evolution of the ring current and of electromagnetic ion cyclotron waves To demonstrate the effects of EMIC wave propagation and refraction on the wave energy distribution and evolution, we simulate the May 1998 storm. The main findings of our simulation can be summarized as follows. First, owing to the density gradient at the plasmapause, the net wave refraction is suppressed, and He+-mode grows preferably at the plasmapause. This result is in total agreement with previous ray tracing studies and is very clearly found in presented B field spectrograms. Second, comparison of global wave distributions with the results from another ring current model (Kozyra et al., 1997) reveals that this new model provides more intense and more highly plasmapause-organized wave distributions during the May 1998 storm period Finally, it is found that He(+)-mode energy distributions are not Gaussian distributions and most important that wave energy can occupy not only the region of generation, i.e., the region of small wave normal angles, but all wave normal angles, including those to near 90 . The latter is extremely crucial for energy transfer to thermal plasmaspheric electrons by resonant Landau damping and subsequent downward heat transport and excitation of stable auroral red arcs.

The further development of a self-consistent theoretical model of interacting ring current ions and electromagnetic ion cyclotron waves [Khazanov et al., 2003] is presented. In order to adequately take into account the wave propagation and refraction in a multi-ion plasmasphere, we explicitly include the ray tracing equations in our previous self-consistent model and use the general form of the wave kinetic equation. This is a major new feature of the present model and, to the best of our knowledge, the ray tracing equations for the first time are explicitly employed on a global magnetospheric scale in order to self-consistently simulate spatial, temporal, and spectral evolutions of the ring current and electromagnetic ion cyclotron waves. To demonstrate the effects of EMIC wave propagation and refraction on the EMIC wave energy distributions and evolution we simulate the May 1998 storm. The main findings of our simulation can be summarized as follows. First, due to the density gradient at the plasmapause, the net wave refraction is suppressed, and He(+)-mode grows preferably at plasmapause. This result is in a total agreement with the previous ray tracing studies, and very clear observed in presented B-field spectrograms. Second, comparison the global wave distributions with the results from other ring current model [Kozyra et al., 1997] reveals that our model provides more intense and higher plasmapause organized distributions during the May, 1998 storm period. Finally, the found He(+)-mode energy distributions are not Gaussian distributions, and most important that wave energy can occupy not only the region of generation, i. e. the region of small wave normal angles, but the entire wave normal angle region and even only the region near 90 degrees. The latter is extremely crucial for energy transfer to thermal plasmaspheric electrons by resonant Landau damping, and subsequent downward heat transport and excitation of stable auroral red arcs.

The self-consistent treatment of the RC ion dynamics and EMlC waves, which are thought to exert important influences on the ion dynamical evolution, is an important missing element in our understanding of the storm-and recovery-time ring current evolution. Under certain conditions, relativistic electrons, with energies greater than or equal to 1 MeV, can be removed from the outer radiation belt by EMlC wave scattering during a magnetic storm (Summers and Thorne, 2003; Albert, 2003). That is why the modeling of EMlC waves is critical and timely issue in magnetospheric physics. This study will generalize the self-consistent theoretical description of RC ions and EMlC waves that has been developed by Khazanov et al. [2002, 2003] and include the heavy ions and propagation effects of EMlC waves in the global dynamic of self-consistent RC - EMlC waves coupling. The results of our newly developed model that will be presented at Huntsville 2006 meeting, focusing mainly on the dynamic of EMlC waves and comparison of these results with the previous global RC modeling studies devoted to EMlC waves formation. We also discuss RC ion precipitations and wave induced thermal electron fluxes into the ionosphere.

We demonstrate extreme ionospheric response to the large interplanetary electric fields during the "Halloween" storms that occurred on October 29 and 30, 2003. Within a few (2 - 5) hours of the time when the enhanced interplanetary electric field impinged on the magnetopause, dayside total electron content increases of approx.40% and approx.250% are observed for the October 29 and 30 events, respectively. During the Oct 30 event, approx.900% increases in electron content above the CHAMP satellite (approx.400 km altitude) were observed at mid-latitudes (+/-30 degrees geomagnetic). The geomagnetic storm-time phenomenon of prompt penetration electric fields is a possible contributing cause of these electron content increases, producing dayside ionospheric uplift combined with equatorial plasma diffusion along magnetic field lines to higher latitudes, creating a "daytime super-fountain" effect.

Modeling results of the inner magnetosphere showing the influence of the ionospheric conductance on the inner magnetospheric electric fields during the April 17, 2002 magnetic storm are presented. Kinetic plasma transport code results are analyzed in combination with observations of the inner magnetospheric plasma populations, in particular those from the IMAGE satellite. Qualitative and quantitative comparisons are made with the observations from EW, MENA, and HENA, covering the entire energy range simulated by the model (0 to 300 keV). The electric field description, and in particular the ionospheric conductance, is the only variable between the simulations. Results from the data-model comparisons are discussed, detailing the strengths and weaknesses of each conductance choice for each energy channel.

The goal of that proposal was to examine the relationship between solar wind drivers and ring current dynamics through data analysis and numerical simulations. The data analysis study was a statistical examination (via superposed epoch analyses) of a solar cycle's worth of storm data. Solar wind data, geophysical indices, and geosynchronous plasma data were collected for every time period with Dst< -50 nT from 1989 through 2002, and the storm list now exceeds 400 entries. This work was first conducted by a summer undergraduate student, Mr. John Vann (University of Kansas), with funding from the NSF Research Experience for Undergraduates program. It was then continued by a University of Michigan graduate student, Mr. Jichun Zhang. Mr. Zhang is now in his fourth year at U-M and is progressing very well toward a PhD in space science. His dissertation will be based on his data analysis and modeling efforts using this geomagnetic storm database. The results of the data analysis study have been the focus of several conference presentations, and the first manuscript has just been published. Two additional papers are presently being prepared, one on average (superposed) solar wind features for various storm subsets (e.g., intense storms at solar maximum), and another on geosynchronous plasma features for these same storm subsets. The latter result was highlighted by the TR&T program director in his presentation at the COSPAR meeting this summer.

Under the proper conditions, guided plasmaspheric hiss is shown to be more efficient than Coulomb collisions at scattering electrons in the superthermal energy range of 50 to 500 eV. Broadband, whistler mode hiss becomes guided by plasma density gradients, intensifying the wave energy densities and focusing the wave normal angles. These waves are shown to interact through Cherenkov (Landau) resonance with electrons below 500 eV, and the presented equatorial plane timescales for pitch angle, energy, and mixed diffusion are shown to be faster than Coulomb collision timescales for typical values at the inner edge of the plasmapause and in detached plasma regions. In the latter case, energy diffusion timescales of less than 100 s for small pitch angle electrons between 250 and 500 eV indicate that these waves have the potential to dramatically change the distribution function.

The unusually high solar wind pressure and strongly southward IMF on May 4, 1998, pushed the magnetopause well into the geosynchronous orbit which exposed the POLAR satellite to the magnetosheath and solar wind. We use a gasdynamic convected field model to predict the magnetosheath quantities and then compare them with the in situ observations. The model prediction helps to reduce the uncertainty in the timing of the solar wind arrival time and provides a reference value for each physical parameter. It also helps to resolve the location of the satellite during strong magnetic fluctuations near the magnetopause. The plasma measurements from the TIDE instrument, in conjunction with the magnetometer measurements, indicate that there is a magnetospheric boundary layer during the event. There are also transient signatures near the magnetopause which may be caused by magnetospheric flux transfer events.

Two time-dependent kinetic models of superthermal electron transport are combined to conduct global calculations of the nonthermal electron distribution function throughout the inner magnetosphere. It is shown that the energy range of validity for this combined model extends down to the superthermal-thermal intersection at a few eV, allowing for the calculation of the en- tire distribution function and thus an accurate heating rate to the thermal plasma. Because of the linearity of the formulas, the source terms are separated to calculate the distributions from the various populations, namely photoelectrons (PEs) and plasma sheet electrons (PSEs). These distributions are discussed in detail, examining the processes responsible for their formation in the various regions of the inner magnetosphere. It is shown that convection, corotation, and Coulomb collisions are the dominant processes in the formation of the PE distribution function and that PSEs are dominated by the interplay between the drift terms. Of note is that the PEs propagate around the nightside in a narrow channel at the edge of the plasmasphere as Coulomb collisions reduce the fluxes inside of this and convection compresses the flux tubes inward. These distributions are then recombined to show the development of the total superthermal electron distribution function in the inner magnetosphere and their influence on the thermal plasma. PEs usually dominate the dayside heating, with integral energy fluxes to the ionosphere reaching 10(exp 10) eV/sq cm/s in the plasmasphere, while heating from the PSEs typically does not exceed 10(exp 8) eV/sq cm/s. On the nightside, the inner plasmasphere is usually unheated by superthermal electrons. A feature of these combined spectra is that the distribution often has upward slopes with energy, particularly at the crossover from PE to PSE dominance, indicating that instabilities are possible.

We have a new model of ion transport that we have applied to the problem of plasmaspheric flux tube refilling after a geomagnetic disturbance. This model solves the Fokker-Planck kinetic equation by applying discrete difference numerical schemes to the various operators. Features of the model include a time-varying ionospheric source, self-consistent Coulomb collisions, field-aligned electric field, hot plasma interactions, and ion cyclotron wave heating. We see refilling rates similar to those of earlier observations and models, except when the electric field is included. In this case, the refilling rates can be quite different that previously predicted. Depending on the populations included and the values of relevant parameters, trap zone densities can increase or decrease. In particular, the inclusion of hot populations near the equatorial region (specifically warm pancake distributions and ring current ions) can dramatically alter the refilling rate. Results are compared with observations as well as previous hydrodynamic and kinetic particle model simulations.

Under the proper conditions, guided plasmaspheric hiss is shown to be more efficient than Coulomb collisions at scattering electrons in the superthermal energy range of 50 to 500 eV. Broadband, whistler mode hiss becomes guided by plasma density gradients, intensifying the wave energy densities and focusing the wave normal angles. These waves are shown to interact through Cherenkov (Landau) resonance with electrons below 500 eV, and the presented equatorial plane timescales for pitch angle, energy, and mixed diffusion are shown to be faster than Coulomb collision timescales for typical values at the inner edge of the plasmapause and in detached plasma regions. In the latter case, energy diffusion timescales of less than 100 s for small pitch angle electrons between 250 and 500 eV indicate that these waves have the potential to dramatically change the distribution function.

A bounced-averaged ring current kinetic model for arbitrary pitch angle, including losses due to charge exchange and Coulomb collisions along ion drift paths, is developed and solved numerically. Results from simplifield model runs, intended to illustrate the effects of adiabatic drifts and collisional losses on the proton population, are presented. The processes of: (1) particle acceleration under the conditions of time-independent magnetospheric electric fields; (2) a predominant loss of particles with small pitch angles due to charge exchange; and (3) a buildup of a low-energy population caused by the Coulomb drag energy degradation, are discussed.

Protons that are convected into the inner magnetosphere in response to enhanced magnetic activity can resonate with ducted plasmaspheric hiss in the outer plasmasphere via an anomalous Doppler-shifted cyclotron resonance. Plasmaspheric hiss is a right-hand-polarized electromagnetic emission that is observed to fill the plasmasphere on a routine basis. When plasmaspheric hiss is confined within field-aligned ducts or guided along density gradients, wave normal angles remain largely below 45 deg. This allows resonant interactions with ions at typical ring current and radiation belt energies to take place. Such field-aligned ducts have been observed both within the plasmasphere and in regions outside of the plasmasphere. Wave intensities are estimated using statistical information from studies of detached plasma regions. Diffusion coefficients are presented for a range of L shells and proton energies for a fixed wave distribution. Harmonic resonances in the range N = +/-100 are considered in order to include interactions between hiss at 100 Hz to 2 kHz frequencies, and protons in the energy range between approximately 10 keV and 1000 keV. Diffusion timescales are estimated to be of the order of tens of days and comparable to or shorter than lifetimes for Coulomb decay and charge exchange losses over most of the energy and spatial ranges of interest.

The decay of the major ion species which constitute the ring current is studied by solving the time evolution of their distribution functions during the recovery phase of a moderate geomagnetic storm. In this work, only equatorially mirroring particles are considered. Particles are assumed to move subject to E x B and gradient drifts. They also experience loses along their drift paths. Two loss mechanisms are considered: charge exchange with neutral hydrogen atoms and Coulomb collisions with thermal plasma in the plasmasphere. Thermal plasma densities are calculated with a plasmaspheric model employing a time-dependent convection electric field model. The drift-loss model successfully reproduces a number of important and observable features in the distribution function. Charge exchange is found to be the major loss mechanism for the ring current ions; however the important effects of Coulomb collisions on both the ring current and thermal populations are also presented. The model predicts the formation of a low-energy (less than 500 eV) ion population as a result of energy degradation caused by Coulomb collision of the ring current ions with the plasmaspheric electrons; this population may be one source of the low-energy ions observed during active and quiet periods in the inner magnetosphere. The energy transferred to plasmaspheric electrons through Coulomb collisions with ring current ions is believed to be the energy source for the electron temperature enhancement and the associated 6300 A (stable auroral red (SAR) arc) emission in the subauroral region. The calculated energy deposition rate is sufficient to produce a subauroral electron temperature enhancement and SAR arc emissions that are consistent with observations of these quantities during moderate magnetic activity levels.

A set of observations describing ionospheric conditions, magnetospheric populations, and 6300-A emission intensities on stable auroral red (SAR) arc field lines during the solar minimum 19-24 Sept. 1984 magnetic storm period prompted a study of solar cycle and magnetic storm phase variations in SAR arc emissions and their magnetospheric energy source. It was found that medium-energy H(+) was significantly enhanced during the main phase compared to the late recovery phase of the 19-20 Sept. 1984 storm. Enhanced heating of the thermal electron plasma caused by this population resulted in more than an order of magnitude greater SAR arc emissions in the main phase compared to the recovery phase. O(+) was found to be the dominant energy source for SAR arcs in the late recovery phases of storms in the 19-24 Sept. period.

A statistical study of the seasonal variations of the subauroral electron temperature enhancement was undertaken using data from the Langmuir probe experiment on the DE 2 satellite throughout most of the mission (1981-1982). In the winter hemisphere the nighttime background electron temperature is the highest and the magnitude of the peak Te responds most weakly to the geomagnetic activity. This behavior can be explained by seasonal trends in the nighttime downward heat flux due to conjugate photoelectrons. Moreover, model results indicate that a factor of about three increase in heat inflow during equinox relative to solstice is required to raise the electron temperature to a given level. This is a consequence of the higher electron densities at the Te peak near equinox. The Te peak occurs on field lines which thread the outer plasmasphere in the vicinity if the plasmapause and thus can be used as a tracer of the plasmapause position.

The elevation of the subauroral electron temperature is one of the phenomena showing the energy transfer from the magnetosphere and the response of the ionosphere. This study addresses solar cycle variations in the subauroral Te peak by comparing observations of the subauroral peak by the Atmosphere Explorer C (AE-C) satellite near solar minimum (1974 and 1977) with similar observations by the Dynamics Explorer 2 (DE 2) satellite near solar maximum (1981-1982). Te peaks with magnitudes sufficient to produce observable stable auroral red arc emissions occurred more frequently during solar maximum than in solar minimum, but the variation in the magnitudes and positions of these peaks with magnetic activity did not change significantly with solar cycle. These results are discussed in terms of the solar cycle changes in the ionosphere and the magnetospheric energy source.

Data from low altitude plasma instrument (LAPI) on Dynamics Explorer 2 document a population of high-energy (up to 800 eV) atmospheric photoelectrons that has not been reported in the published literature. The source of these photoelectrons is postulated to be the soft X-ray portion of the whole sun spectrum. This conclusion is supported by sunrise-sunset characteristics that track those of the classical (below 60 eV) EUV-produced photoelectrons, and theoretical results from two models that incorporate the soft X-ray portion of the solar spectrum. The models include K-shell ionization effects and predict peaks in the photoelectron spectrum due to Auger electrons emitted from oxygen and nitrogen. The peak for nitrogen is observed as predicted, but the peak for oxygen is barely observable. Excellent quantitative agreement is achieved between theory and experiment by using reasonable adjustments to the few published soft X-ray spectra based on solar activity. The upflowing energetic photoelectrons provide a heretofore unknown source of electrons to the magnetosphere. They occur whenever and wherever the sun is up, that is, at all invariant latitudes. Their density is low, but they are steady and ubiquitous. If scattering and trapping occur on closed field lines, then photoelectrons could contribute as a significant particle source and thus represent a new facet of magnetosphere-ionosphere coupling.

A case study has been carried out using data from the Dynamics Explorer 1 and 2 spacecraft to study the effect of Coulomb interactions between ring current and suprathermal O(+) and thermal protons on the plasmasphere. Results from a one-dimensional plasmaspheric model suggest that heating due to Coulomb collisions may be sufficient to raise the ion and electron temperatures to observed values. The resultant high temperature produced enhancements in the model O(+) and O(++) densities in agreement with observations.

The solar wind interacts very strongly with the extensive cometary coma, and the various interaction processes are initiated by the ionization of cometary neutrals. The main ionization mechanism far outside the cometary bow shock is photoionization by solar extreme ultraviolet radiation.Electron distributions measured in the vicinity of comets Halley and Giacobini-Zinner by instruments on the VEGA and ICE spacecraft, respectively, are used to calculate electron impact ionization frequencies. Ionization by electrons is of comparable importance to photoionization in the magnetosheaths of Comets Halley and Giacobini-Zinner. The ionization frequency in the inner part of the cometary plasma region of comet Halley is several times greater than the photoionization value. Tables of ionization frequencies as functions of electron temperature are presented for H2O, CO2, CO, O, N2, and H.

Observations of thermal and energetic ion populations by the Dynamics Explorer 1 satellite in the magnetospheric energy source region as well as nearly simultaneous Langmuir probe measurements of enhanced electron temperatures by Dynamics Explorer 2 within the stable auroral red (SAR) arc at F region heights are used to examine the role of heavy ions in the formation of SAR arcs. It is found that sufficient energy is transferred to the electron gas at high altitudes via Coulomb collisions between the observed ring current ions and thermal electrons to support the enhanced F region electron temperatures measured on these field lines. The latitudinal variation in the electron heating rates calculated using observed ion populations is consistent with the observed variation in electron temperature across the SAR arc. In all cases, ring current O(+) is the major source of energy for the SAR arcs.

Results of coincident measurements by ground-based photometers and the low-altitude plasma instrument on board the Dynamics Explorer 2 satellite are presented that demonstrate the association of precipitating low-energy electrons with stable auroral red (SAR) arcs. A search of available data has yielded 23 instances of DE 2 overflights during the presence of SAR arcs being monitored by the photometers. For each case, downward fluxes of electrons are found to be enhanced along field lines penetrating the arcs in relation to regions both north and south of the features. Modeling of the atmospheric response to these influxes indicates that these electrons can represent a major source of the energy required to establish temperature profiles within the ionospheric electron gas that are sufficient to produce the recorded 6300-A emission rates. The sensitivity of these results to uncertainties of the assumed spacecraft potential and thermospheric composition has been investigated and found to be important, but does not alter the conclusion that precipitating electrons are a fundamental link in the production of SAR arcs.