A new numerical code for simulating current-driven instabilities on auroral field lines

A new code is described with open boundary conditions for simulating electrostatic instabilities driven by field-aligned currents in narrow channels whose widths are smaller than the plasma skin depth. The reflection of Alfven waves at regions containing plasma waves coupling the electron and ion motion is found to cause a reduction of the field-aligned current and limits it to a critical value which is marginal for the excitation of electrostatic ion-cyclotron waves at a height of 1-2 R(e). Electron heating due to current-driven plasma turbulence causes the overall plasma density in the system to decrease. The plasma density in the region of upward current is enhanced relative to that in the region of downward current, and the current density in upward-directed current channels can be larger than in downward-directed ones.