Engineering PapersSearch

Engineering topics

Hudson, M. K.

Publications and source records attributed to Hudson, M. K..

At least 19 records

Gradual Diffusion and Punctuated Phase Space Density Enhancements of Highly Relativistic Electrons: Van Allen Probes Observations

The dual-spacecraft Van Allen Probes mission has provided a new window into mega electron volt (MeV) particle dynamics in the Earth's radiation belts. Observations (up to E (is) approximately 10MeV) show clearly the behavior of the outer electron radiation belt at different timescales: months-long periods of gradual inward radial diffusive transport and weak loss being punctuated by dramatic flux changes driven by strong solar wind transient events. We present analysis of multi-MeV electron flux and phase space density (PSD) changes during March 2013 in the context of the first year of Van Allen Probes operation. This March period demonstrates the classic signatures both of inward radial diffusive energization and abrupt localized acceleration deep within the outer Van Allen zone (L (is) approximately 4.0 +/- 0.5). This reveals graphically that both 'competing' mechanisms of multi-MeV electron energization are at play in the radiation belts, often acting almost concurrently or at least in rapid succession.

Van Allen Probes

Latitudinal Density Dependence of Magnetic Field Lines Inferred from Polar Plasma Wave Data

Using observations of the electron density, n(sub e), based on measurement of the upper hybrid resonance frequency by the Polar spacecraft Plasma Wave Instrument, we have examined the radial density dependence along field lines in the outer plasmasphere and the near plasmatrough. Sampled L values range from 2.5 to 6.6. Our technique depends on the fact that Polar crosses particular L values at two different points with different radial distance R. In our plasmaspheric data set (n(sub e) > 100/cm3), we find that on average n(sub e) is flat along field lines from the equator up to the latitudes sampled by Polar (R approximately equal to or > 2.0). In the plasmatrough data set (n(sub e) < 100/cm-3), there is on average a mild radial dependence n(sub e) varies as R(exp -1.7).

Goldstein, J.

Ponderomotive effects on distributions of O(+) ions in the auroral zone

Test particle calculations are used to compute the effects of gravity and ponderomotive acceleration by shear Alfven wave oscillations on the distribution function of O(+) ions along auroral field lines, assuming an ionospheric Maxwellian source of the ions at 2000 km altitude with approximately 0.5 eV of thermal energy in the parallel component of velocity. The electric field model corresponds to a standing wave oscillation with a frequency approximately 1 Hz in the azimuthal direction superimposed on the background dipole field, in which the wave amplitude is either increasing or decreasing in time. The electric field is taken to be primarily in the perpendicular direction. The time varying wave produces broad distributions with widths of 2 to 10 times the initial 0.5-eV thermal energy of the Maxwellian source, and the density and flux of upward going O(+) ions at one Earth radius are both enhanced in this model. The oxygen ion distribution functions at 1 R(sub E) altitude resulting from interaction with waves whose amplitudes are increasing in time have a more gradual lower energy cutoff than do the distribution functions resulting from decaying waves. The high-energy part of the distribution functions in growing waves reflects the temperature of the Maxwellian source, while the high-energy part of the distributions resulting from decaying waves steepens with time, independent of the source temperature.

Witt, E.

Simulation of proton radiation belt formation during the March 24, 1991 SSC

The rapid formation of a new proton radiation belt at L approximately = 2.5 following the March 24, 1991 Storm Sudden Commencement (SSC) observed at the Combined Release and Radiation Effects Satellite (CRRES) satellite is modeled using a relativistic guiding center test particle code. The SSC is modeled by a bipolar electric field and associated compression and relaxation in the magnetic field, superimposed on a dipole magnetic field. The source population consists of both solar and trapped inner zone protons. The simulations show that while both populations contribute to drift echoes in the 20-80 MeV range, primary conditions is from the solar protons. Proton acceleration by the SSC differs from relativistic electron acceleration in that different source populations contribute and nonrelativistic conservation of the first adiabatic invariation leads to greater energization of protons for a given decrease in L. Model drift echoes and flux distribution in L at the time of injection compare well with CRRES observations.

Hudson, M. K.

An MHD simulation study of the poloidal mode field line resonance in the Earth's dipole magnetosphere

The poloidal mode field line resonance in the Earth's dipole magnetic field is investigated using cold plasma ideal MHD simulations in dipole geometry. In order to excite the poloidal mode resonance, we use either an initial or a continuous velocity perturbation to drive the system. The perturbation is localized at magnetic shell L = 7 with plasma flow in the radial direction (electric field component in the azimuthal direction). It is found that with the initial perturbation alone, no polodial mode resonance can be obtained and the initially localized perturbation spreads out across all magnetic L shells. With the continuous perturbation, oscillating near the poloidal resonance frequency, a global-scale poloidal cavity mode can be obtained. For the first time, a localized guided poloidal mode resonance is obtained when a radial component of electric field is added to the initial perturbation such that the curl of the electric field is everywhere perpendicular to the background dipole magnetic field. During the localized poloidal resonance, plasma vortices parallel/antiparallel to the background dipole magnetic field B(sub 0). This circular flow, elongated radially, results in twisting of magnetic field flux tubes, which, in turn, leads to the slowdown of the circular plasma flow and reversal of the plasma vortices. The energy associated with the localized poloidal resonance is conserved as it shifts back and forth between the oscillating plasma vortices and the alternately twisted magnetic flux tubes. In the simulations the eigenfunctions associated with the localized poloidal resonance are grid-scale singular functions. This result indicates that ideal MHD is inadequate to describe the underlying problem and nonideal MHD effects are needed for mode broadening.

Ding, D. Q.

Simulation of the prompt energization and transport of radiation belt particles during the March 24, 1991 SSC

We model the rapid (about 1 min) formation of a new electron radiation belt at L about or = 2.5 that resulted from the Storm Sudden Commencement (SSC) of March 24, 1991 as observed by the Combined Release and Radiation Effects Satellite (CRRES) satellite. Guided by the observed electric and magnetic fields, we represent the time-dependent magnetospheric electric field during the SSC by an asymmetric bipolar pulse that is associated with the compression and relaxation of the Earth's magnetic field. We follow the electrons using a relativistic guiding center code. The test-particle simulations show that electrons with energies of a few MeV at L greater than 6 were energized up to 40 MeV and transported to L about or = 2.5 during a fraction of their drift period. The energization process conserves the first adiabatic invariant and is enhanced due to resonance of the electron drift motion with the time-varying electric field. Our simulation results, with an initial W(exp -8) energy flux spectra, reproduce the observed electron drift echoes and show that the interplanetary shock impacted the magnetosphere between 1500 and 1800 MLT.

Li, Xinlin

Analysis of cyclotron harmonic emissions at the outer planets

The flyby missions of Voyagers 1 and 2 at Jupiter, Saturn and Uranus revealed intense waves above the electron gyrofrequency. Observation of waves at the upper hybrid frequency is often accompanied by power at adjacent electron Bernstein harmonics, and the relative power in these modes depends both on the density and temperature ratios of the cold background electron population and the hot magnetospheric electrons which drive the instability. A model of electron distributions which is consistent with observations is used for analysis of the excited waves, their dependence upon plasma parameters, and the time scales of the saturation processes. It is shown that in the presence of two-temperature electron distributions the linear excitation is due to a fluidlike coupling of two eigenmodes for perpendicular propagation and to kinetic destalization of oblique modes. The dependence of linear growth rates on propagation angle is presented, along with results from particle simulations. A quasi-linear diffusion time for relaxing the hot electron loss cone is calculated and compared with simulation results. This time scale is faster than for local saturation by heating the cold population, and also the convective amplification time scale, suggesting that the waves saturate at quasi-linear levels, while being convectively localized to the equatorial regions of the outer planetary magnetospheres.

Roth, I.

Lower hybrid waves excited through linear mode coupling and the heating of ions in the auroral and subauroral magnetosphere

Two related aspects of the linear mode conversion mechanism are considered: (1) the controlled heating of suprathermal ions in the ionosphere and magnetosphere over a powerful VLF/ELF transmitter using lower hybrid (LH) waves excited by the electromagnetic (EM) transmitter signals through linear mode conversion; (2) the excitation of intense LH waves in the auroral regions through the linear conversion of VLF/ELF EM auroral hiss, and the subsequent heating of ions by the excited LH waves. Experimental results are used to demonstrate that strong LH waves can be excited by EM waves down to frequencies as low as 2 kHz in the subauroral low-altitude magnetosphere. It is also shown that an approximately 250 kW VLF/ELF transmitter operating in the subauroral region at 2 kHz could excite sufficiently intense LH waves to heat suprathermal 1 eV H(+) ions and 16 eV O(+) ions to roughly 50 eV in a region extending in altitude from 1000 to 5000 km above the transmitter with a horizontal scale of about 500 km. If coherent wave stochastic heating occurs, the O(+) energy gain could be as large as 200 eV.

Bell, T. F.

Lower hybrid drift instability at the inner edge of the ring current

Excitation of electrostatic waves around and above the lower hybrid frequency at the inner edge of the ring current is investigated. The model consists of cold plasmaspheric ions; hot, dilute, inhomogeneous ring current ions; and cold charge neutralizing electrons. It is shown that in the presence of a sufficiently large density gradient the coupling of the lower hybrid with the drift cyclotron harmonic mode destabilizes the plasma. Unstable interactions also occur at higher cyclotron harmonics when the lower hybrid is replaced by an ion Bernstein mode. Additional power at frequencies above the lower hybrid frequency is attributed to the thermal enhancement by the hot ring current ions, which increases with their density gradient. Particle simulations are conducted, and comparison to AMPTE data is discussed.

Roth, I.

Particle simulation of ion heating in the ring current

Heating of heavy ions has been observed in the equatorial magnetosphere in GEOS 1 and 2 and ATS 6 data due to ion cyclotron waves generated by anisotropic hot ring current ions. A one-dimensional hybrid-Darwin code has been developed to study ion heating in the ring current. Here, a strong instability and heating of thermal ions is investigated in a plasma with a los cone distribution of hot ions. The linear growth rate calculation and particle simulations are conducted for cases with different loss cones and relative ion densities. The linear instability of the waves, the quasi-linear heating of cold ions and dependence on the thermal H(+)/He(+) density ratio are analyzed, as well as nonlinear parallel heating of thermal ions. Effects of thermal oxygen and hot oxygen are also studied.

Qian, S.

Generation models of electron conics

Electron distribution functions (EDFs) with a peak oblique to the magnetic field, adjacent to but distinct from loss-cone features, have been observed by the DE 1, Viking, and S3-3 spacecraft in passes through the nightside auroral zone, polar cap, dayside cusp, and extended dayside auroral oval. Using particle simulations, two types of wave excitation and particle acceleration mechanisms which may contribute to producing these electron conic distributions are investigated. The first involves excitation of upper-hybrid waves by the electron loss cone, and the subsequent perpendicular heating of the background and thermal electrons. The second involves excitation of downward propagating parallel modes by an auroral electron beam which frequently accompanies the upflowing electron conics. These modes provide parallel acceleration, which modifies the GDF. Those electrons which are not lost to the atmosphere and mirror back up the magnetic field line give rise to enhancements in the GDF at the edge of the loss cone.

Roth, I.

Ionospheric dynamo theory for production of far ultraviolet emissions on Uranus

A model is presented to explain diffuse FUV emissions from the outer planets, specifically Uranus, in excess of those diffuse emissions that are currently explainable by scattering of sunlight and/or excitation by photoelectrons. These electroglow emissions in H Ly-alpha and H2 bands, which occur in the sunlit hemisphere slightly above the homopause, appear to require particle excitation in the 10- to 50-eV range. An in situ mechanism for accelerating photoelectrons (and ions is proposed, involving neutral wind dynamo generation of field-aligned currents analogous to what occurs in the earth's equatorial E and F regions. Sufficiently strong field-aligned currents are found in the model calculation for Uranus to produce a potential drop of about 100 eV or greater between the F peak and homopause, concentrated at lower altitudes, and capable in principle of accelerating photoelectrons (and ions) to the 10- to 50-eV energies required to explain the observed emissions. The fact that the excitation and ionization cross sections are larger than elastic scattering cross sections in an H2 atmosphere at these energies makes in situ acceleration feasible for the production of UV on the outer planets.

Hudson, M. K.

Oxygen acoustic solitary waves in a magnetized plasma

Ion-acoustic solitary waves in a magnetized plasma containing an arbitrary mixture of H(+) and O(+) ions are studied. A nonlinear wave equation has been derived from the Poisson-Vlasov equations, including a uniform magnetic field and dissipation due to reflected electrons. When dissipation is ignored, the equation has soliton solutions associated with both oxygen and hydrogen acoustic modes, which can be either rarefactive or compressive depending on the ion concentrations and the electron/ion temperature ratio and, more weakly, on the bulk drifts of the species. If electron reflection is included, the solitary wave can be intensified. Under somewhat restrictive conditions the oxygen solitary wave is rarefactive and propagates with a velocity comparable to that observed by the Viking satellite. The three-dimensional solitons obey a relation of scales parallel to the magnetic field and in the transverse direction. Computer simulations of one-dimensional versions of the nonlinear wave equation are presented.

Qian, S.

Effects of ion two-stream instability on auroral ion heating

The ion two-stream instability develops between upflowing H(+) and O(+) in the auroral acceleration region, if the two species are accelerated to different velocities in a semistationary electrostatic potential. Parallel propagating modes are unstable only for small relative drifts, while oblique modes remain unstable at larger drifts and, consequently, higher altitudes. The effects of oblique modes on ion heating are investigated, and it is found that O(+) heats primarily by trapping in the direction of wave propagation during the early, coherent phase of wave growth in the temporal evolution problem. After undergoing trapping along the magnetic field direction, H(+) heats primarily by quasi-linear diffusion when the spectrum broadens nonlinearly. The parallel ion distributions evolve by formation of a high-energy (low-energy) tail on oxygen (hydrogen) along with bulk slowing of hydrogen, which results in more energetic oxygen than hydrogen. Total heating is greatest in the direction of wave propagation, and this effect is greater for hydrogen than for oxygen.

Roth, I.

Simulations of active ion injection experiments on ARCS 3

The ARCS 3 sounding rocket experiment used two argon guns oriented parallel and perpendicular to the geomagnetic field. The parameters of this experiment have been used in a series of particle simulations in order to interpret ion beam interaction with the multispecies ionospheric plasma. Along with the lower hybrid mode between the electron and the hydrogen gyrofrequencies, ion-ion hybrid modes between the hydrogen and helium gyrofrequencies and between the helium and oxygen gyrofrequencies are identified in the simulations. Besides these cold plasma modes, warm Bernstein modes of the respective ion species are identified. The parallel argon beam injections excite lower hybrid waves perpendicular to the magnetic field because of the finite thermal spread of the parallel beam. Perpendicular ion heating of background hydrogen, helium, and oxygen is comparable in the simulations. Comparisons are made with the wave and particle data from the ARCS 3 rocket experiment.

Hudson, M. K.

Dynamics of localized ion-acoustic waves in a magnetized plasma

The evolution of negative potential pulses in a magnetized plasma is studied. A three-dimensional nonlinear ion-acoustic wave equation, including nonstationary effects of reflected electrons, has been derived from the Poisson-Vlasov equations with uniform magnetic field. The studies show that a negative potential pulse can be enhanced by drifting electrons. The growing pulse develops asymmetrically with an oscillatory precursor and a local potential jump resembling the early phase of weak double layer formation.

Qian, S.

Linear stability of the H(+) - O(+) two-stream interaction in a magnetized plasma

The Bergman and Lotko (1986) analysis of the linear stability of the H(+) - O(+) two-stream interactions for parallel propagating wave modes is extended to an analysis of the stability of wave modes that propagate at an angle to the terrestrial magnetic field. The analaysis provides an insight into the manner in which the oblique modes linearly couple and become unstable; the conditions under which the particular couplings are important; the qualitative behavior of the frequency, wave vector, and growth rate of the most unstable mode; and the changes to these introduced by finite ion temperatures. Consideration is also given to the relative importance of obliquely propagating and parallel propagating waves for the condition when the velocity difference between two ion beams is below the upper limit.

Bergmann, R.

Weak double layers in the auroral ionosphere

Previous work on the evolution of weak double layers in a hydrogen plasma was extended to include H(+) and O(+) with relative drift. The relative drift between hydrogen and oxygen ions due to a quasi-static parallel electric field gives rise to a strong linear fluid instability which dominates the ion-acoustic mode at the bottom of the auroral acceleration region. This ion-ion instability can modify ion distributions at lower altitudes and the subsequent nonlinear evolution of weak double layers at higher altitudes in the ion-acoustic regime. Ion hole formation can occur for smaller relative electron-ion drifts than seen in previous simulations, due to the hydrogen-oxygen two-stream instability. This results in local modification of the ion distributions in phase space, and a partial filling of the valley between the hydrogen and oxygen peaks, which would be expected at higher altitudes on auroral field lines. The observed velocity diffusion does not necessarily preclude ion hole and double layer formation in hydrogen in the ion-acoustic regime. These simulation results are consistent with the experimentally measured persistence of separate hydrogen and oxygen peaks, and the observation of weak double layers above an altitude of 3000 km on auroral field lines.

Hudson, M. K.