Search NASA⌕ Search

SEARCH · Search NASA

Results for “IONOSPHERIC HEATING”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRS2006

Self-Consistent Model of Magnetospheric Ring Current and Propagating Electromagnetic Ion Cyclotron Waves: Waves, Precipitating Ring Current Ions, and Thermal Electron Heating - 2

This paper is dedicated to further presentations and discussions of the results from our new global self-consistent theoretical model of interacting ring current ions and electromagnetic ion cyclotron waves [Khazanov et al., 2006; here referred to as Paper 1]. 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 [for details see Paper 1]. To demonstrate the effects of the EMIC wave propagation and refraction on the RC proton precipitations and heating of the thermal plasmaspheric electrons we simulate the May 1998 storm. The main findings of our simulation can be summarized as follows. Firstly, the wave induced precipitations have a quite fine structure, and are highly organized by location of the plasmapause gradient. The strongest fluxes of about 4 (raised dot) 10(exp 6) [(cm (raised dot) s (raised dot) sr)(sup -l)] are observed during the main 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 simply connected to the most intense EMIC waves. The character of the EMIC wave power spectral density distribution over the equatorial wave normal angle is an extremely crucial for the effectiveness of the RC ion scattering. Secondly, comparison of the global proton precipitating patterns with the results from other ring current model [Kozyra et al., 1997] reveals that although we observe a qualitative agreement between localizations of the wave induced fluxes in the models, there is no quantitative agreement between the magnitudes of these fluxes. These differences are mainly due to a qualitative difference between the characters of the EMIC wave power spectral density distributions over the equatorial wave normal angle. Finally, the two energy sources to the plasmaspheric electrons are considered; (i) the heat fluxes caused by the EMIC wave energy absorption due to Landau resonance, and (ii) the heat fluxes due to Coulomb energy degradation of the RC o(+) ions. The heat fluxes caused by the EMIC wave energy absorption due to Landau resonance are observed in the postnoon-premidnight MLT sector, and maximize at the magnitude of 10l1 (eV/(cm(sup 2)(raised dot) s) at L=3.25, MLT=22 at 3400 UT after 1 May, 0000 UT. The greatest Coulomb energy deposition rates are about 2 (raised dot) 10(sup 10)(eV/(cm(sup 2)(raised dot) s) and 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 RC with plasmasphere is strongly supported by concurrent and conjugate plasma measurements from the plasmasphere, the RC, and the topside ionosphere [Gurgiolo et al., 20051.

Khazanov, G. V.↗
NASA NTRS1981

Upper atmosphere of Jupiter - A post Voyager perspective

Some significant upper atmosphere features of Jupiter are discussed with reference to the Voyager observations of 1979. From the time of Pioneer observations in 1973-1974, the Jovian upper atmospheric temperature has increased by about 30%, the eddy diffusion coefficient at the homopause decreased by a factor of 100, the equatorial disk Lyman alpha intensity increased by a factor of 30, and the equatorial ionosphere became more extensive and showed diurnal variation. Important potential candidates for upper atmosphere energetics are thought to be the penetration of the magnetospheric soft electrons, Joule heating mechanism, and auroral electrons.

Atreya, S. K.↗
NASA NTRS2018

Collisionless Relaxation of the Ion Ring Distribution in Space Plasma

Energetic processes often produce transversely-heated angular distributions of the magnetized core (lowest energy) plasma. This characteristic is found in solar wind ion pickup, resulting from cometary or interstellar gas ionization, in Earths' ionosphere, and with hot ions formed around the Space Transportation System during gas releases. We investigate the thermalization of O+ ion pickup using the 2.5D hybrid simulation method (with fluid electrons and kinetic ions) of the ion pickup (ring) distributions, formed in the auroral ionosphere, with a range of ring velocities and thermal to magnetic pressure ratios. We find that in the unstable collisonless regime the anisotropy of the non-thermal distribution produces the ion-cyclotron instability, and the nonlinear relaxation is accompanied by wave-particle scattering that results in an emitted power of EMIC waves. We conclude that ionospheric pickup thermalization is slow due to the small ring speed compared to the thermal and Alfven speeds, while in the solar wind and other space plasmas regions with larger ion-ring velocity the collisionless relaxation and thermalization is rapid in terms of O+ ion gyro-period.

Ofman, Leon↗
NASA NTRS1994

Transversely accelerated ions in the topside ionosphere

Data from the rocket campaigns Mechanism in the Auroral Region for Ion Energization (MARIE) and TOpside Probe of the Auroral Zone (TOPAZ) III, within regions of low-altitude transversely accelerated ions, are interpreted to explain the acceleration of the ions. Using the Monte Carlo kinetic technique to evaluate the ion heating produced by the simultaneously observed lower hybrid waves, we find that their observed electric field amplitudes are sufficient to explain the observed ion energies in the MARIE event. Much of the uncertainty in evaluating the efficiency of a plasma wave induced particle heating process which is dependent on a velocity resonance comes from the lack of information on the phase velocities of the waves. In the case of the MARIE observations, our modeling efforts show that features in the ion velocity distribution are consistent with the wave phase velocities inferred from interferometer measurements of wavelengths. The lower hybrid waves with which low-altitude transversely accelerated ions are associated are frequently observed to be concentrated in small-scale wave packets called 'spikelets'. We demonstrate through the scaling of the size of these wave packets that they are consistent with the theory of lower hybrid collapse. Using the Monte Carlo technique, we find that if the lower hybrid field energy is concentrated in these wave packets, it is still adequate to accelerate the ionospheric ions to the observed energies.

Retterer, John M.↗
NASA NTRS1994

Effect of high-latitude ionospheric convection on Sun-aligned polar caps

A coupled magnetospheric-ionospheric (M-I) magnetohydrodynamic (MHD) model has been used to simulate the formation of Sun-aligned polar cap arcs for a variety of interplanetary magnetic field (IMF) dependent polar cap convection fields. The formation process involves launching an Alfven shear wave from the magnetosphere to the ionosphere where the ionospheric conductance can react self-consistently to changes in the upward currents. We assume that the initial Alfven shear wave is the result of solar wind-magnetosphere interactions. The simulations show how the E region density is affected by the changes in the electron precipitation that are associated with the upward currents. These changes in conductance lead to both a modified Alfven wave reflection at the ionosphere and the generation of secondary Alfven waves in the ionosphere. The ensuing bouncing of the Alfven waves between the ionosphere and magnetosphere is followed until an asymptotic solution is obtained. At the magnetosphere the Alfven waves reflect at a fixed boundary. The coupled M-I Sun-aligned polar cap arc model of Zhu et al.(1993a) is used to carry out the simulations. This study focuses on the dependence of the polar cap arc formation on the background (global) convection pattern. Since the polar cap arcs occur for northward and strong B(sub y) IMF conditions, a variety of background convection patterns can exist when the arcs are present. The study shows that polar cap arcs can be formed for all these convection patterns; however, the arc features are dramatically different for the different patterns. For weak sunward convection a relatively confined single pair of current sheets is associated with the imposed Alfven shear wave structure. However, when the electric field exceeds a threshold, the arc structure intensifies, and the conductance increases as does the local Joule heating rate. These increases are faster than a linear dependence on the background electric field strength. Furthermore, above the threshold, the single current sheet pair splits into multiple current sheet pairs. For the fixed initial ionospheric and magnetospheric conditions used in this study, the separation distance between the current pairs was found to be almost independent of the background electric field strength. For either three-cell or distorted two-cell background convection patterns the arc formation favored the positive B(sub y) case in the northern hemisphere.

Sojka, J. J.↗
NASA NTRS1994

Ionospheric response to the sustained high geomagnetic activity during the March '89 great storm

A simulation was conducted to model the high-latitude ionospheric to the sustaied level of high geomagnetic activity for the great magnetic storm period of March 13-14, 1989. The geomagnetic and solar activity indices and the Defense Meterological Satellite Program (DMSP) F8 and F9 satellite data for particle precipitation and high-latitude convection were used as inputs to a time-dependent ionospheric model (TDIM). The results of the TDIM were compared to both DMSP plasma density data and ground-based total electron content (TEC) measurements for the great storm period as well as with earlier storm observations. The comparisons show that the overall structure of the high-latitude ionosphere was dominated by an increased convection speed within the polar cap that led to increased ion temperatures. In turn, this enhanced the NO(+) density, raised the atomic-to-molecular ion transition height to over 300 km, decreased N(sub m)F(sub 2), increased h(sub m)F(sub 2), and in places either increased n(sub e) at 800 km or slightly decreased it. The morphology of the ionosphere under these extreme conditions was considerably different than that modeled for less distributed intervals. These differences included the character of the dayside tongue of ionization that no longer extended deep into the polar cap. Instead, as a result of the ion heating and consequent reduction in N(sub m)F(sub 2), a large polar hole occupied much of the polar region. This polar hole extended beyond the auroral oval and merged with the night sector midatitude trough. The limitaions associated with the applicability of the TDIM to the geomagnetic conditions present on March 13 and 14 are discussed. The primary limitations of the TDIM derive from the limited temporal resolution of the model input parameters and the lack of suitably dynamic thermospheric specification for the great storm conditions. These limitations leads to midlatitude ionospheric storm phases that do no follow those observed.

Sojka, J. J.↗
NASA NTRS1986

The distribution of atomic hydrogen in the Jovian atmosphere

An analysis is presented of the Voyager and IUE lyman alpha spectra of the Jovian equatorial emission in which was derived a zonal asymmetry in the hydrogen column abundance. Using two estimates of the fraction of Lyman alpha which is due to direct excitation by charged particle precipitation from the ionosphere, upper and lower limits were derived to the H column abundance within and without the perturbed region. That the asymmetry in H abundance may be due to localized heating near the homopause with a consequent rise in scale height is shown. The derived exospheric temperature remains fairly constant with longitude. The required additional heat input over the bulge region, 0.02 erg/cm/s, is supplied by an additional flux of magnetospheric electrons due to Jupiter's magnetic anomaly.

Killen, R. M.↗
NASA NTRS1994

Energy transfer between energetic ring current H(+) and O(+) by electromagnetic ion cyclotron waves

Electromagnetic ion cyclotron (EMIC) waves in the frequency range below the helium gyrofrequency can be excited in the equatorial region of the outer magnetosphere by cyclotron resonant instability with anisotropic ring current H(+) ions. As the unducted waves propagate to higher latitudes, the wave normal should become highly inclined to the ambient magnetic field. Under such conditions, wave energy can be absorbed by cyclotron resonant interactions with ambient O(+), leading to ion heating perpendicular to the ambient magnetic field. Resonant wave absorption peaks in the vicinity of the bi-ion frequency and the second harmonic of the O(+) gyrofrequrency. This absorption should mainly occur at latitudes between 10 deg and 30 deg along auroral field lines (L is greater than or equal to 7) in the postnoon sector. The concomitant ion heating perpendicular to the ambient magnetic field can contribute to the isotropization and geomagnetic trapping of collapsed O(+) ion conics (or beams) that originate from a low-altitude ionospheric source region. During geomagnetic storms when the O(+) content of the magnetosphere is significantly enhanced, the absorption of EMIC waves should become more efficient, and it may contribute to the observed acceleration of O(+) ions of ionospheric origin up to ring current energies.

Thorne, Richard M.↗
NASA NTRS2000

Investigating the Role of the Earth's Ionosphere in Space Weather: Modeling and Observations of High-Latitude Ionospheric Outflows

We will present what we believe to be the first reported observations of up- and downflows of topside ionospheric thermal plasmas from multiple near-simultaneous tracks through the high-latitude topside ionosphere. From several Southern polar passes, it has been possible to construct plots of field-aligned flows of 0+ observed by the Thermal Ion Dynamics Experiment(TIDE) on the POLAR spacecraft near 5000 km altitude together with vertical ion flow observations from one or more DMSP spacecraft near 800 km altitude. We also will present simultaneous observations of POLAR auroral LTVI images with field-aligned flows. We will further show new fluid-kinetic simulations of auroral ionospheric flows in the altitude range 120 km to several R-E as subjected to the synergistic auroral effects of soft electron precipitation, transverse wave-driven ion heating, and hot plasma-driven electric potentials.

Horwitz, J. L.↗
NASA NTRS1985

A theoretical study of the global F region for June solstice, solar maximum, and low magnetic activity

A time-dependent, three-dimensional, multi-ion model of the ionospheric F region at 120-800 km altitude is presented. Account is taken of field-aligned diffusion, cross-field electrodynamic drifts in equatorial and high latitude regions, interhemispheric flow, thermospheric winds, polar wind escape, energy-dependent chemical reactions and neutral composition changes. Attention is also given to the effects of ion production by solar EUV radiation and auroral precipitation, thermal conduction, diffusion-thermal heat flow, local heating and cooling processes, offsets between the geomagnetic and geographic poles, and bending of field lines near the magnetic equator. The model incorporates all phenomena described by previous models and can be applied to tracing magnetic storm and substorm disturbances from high to low latitudes on a global scale. Sample results are provided for ionospheric features during a June solstice, the solar maximum and in a period of low geomagnetic activity. The model will eventually be used to study coupled ionosphere-thermosphere activity.

Sojka, J. J.↗
NASA NTRS1994

A time-dependent gyro-kinetic model of thermal ion upflows in the high-latitude F region

Ample evidence supports the significance of the high-latitude ionospheric contribution to magnetospheric plasma. Assuming flux conservation along a flux tube, the upward field-aligned ion flows observed in the magnetosphere require high-latitude ionospheric field-aligned ion upflows of the order of 10(exp 8) to 10(exp 9)/sq cm/s. Since radar and satellite observations of high-latitude F region flows at times exceed this flux requirement by an order of magnitude, the thermal ionospheric upflows are not simply the ionospheric response to a magnetospheric flux requirement. Several ionospheric ion upflow mechanisms have been proposed, but simulations based on fluid theory do not reproduce all the observed features of ionospheric ion upflows. Certain asymmetries in the statistical morphology of high-latitude F region ion upflows suggest that the ion upflows may be generated by ion-neutral frictional heating. We developed a single-component (O(+)), time-dependent gyro-kinetic model of the high-latitude F region response to frictional heating in which the neutral exobase is a discontinuous boundary between fully collisional and collisionless plasmas. The concept of a discontinuous neutreal exobase and the assumption of a constant and uniform polarization electric field reduce the ion velocity distribution function, from which we can compute the ion density, parallel velocity, parallel and perpendicular temperature, and parallel flux. Using our model, we simulated the response of a convecting flux tube between 500 km and 2500 km to various frictional heating inputs; the results were both qualitatively and quantitatively different from fluid model results, which may indicate an inadequacy of the fluid theory approach. The gyro-kinetic frictional heating model responses to the various simulations were qualitatively similar: (1) initial perturbations of all the modeled parameters propagated rapidly up the flux tube, (2) transient values of the ion parallel velocity, temperature, and flux exceeded 3 km/s, 2 x 10(exp 4) K, and 10(exp 9)/sq cm/s, respectively, (3) a second transient regime developed wherein the parallel temperature drops to very low values (a few hundred Kelvins), and (4) well after heating ceased, large parallel temperatures and large downward parallel velocities and fluxes developed as the flux tube slowly returned to diffusive equilibrium. The ion velocity distributions during the simulation are often non-Maxwellian and are sometimes composed of two distinct ion populations.

Loranc, M.↗
NASA NTRS2010

Mechanisms of Ionospheric Mass Ejection

Ionospheric outflows are directly responsive to solar wind disturbances, particularly in the dayside auroral cusp or cleft regions. Inputs of both electromagnetic energy (Poynting flux) and kinetic energy (particle precipitation) are closely correlated with these outflows. We assess the importance of processes thought to drive ionospheric outflows. These begin with the diffuse effects of photoionization and thermal equilibrium of the ionospheric topside, enhancing Jeans' escape, with ambipolar diffusion and acceleration. Auroral outflows begin with dayside reconnexion and resultant field-aligned currents and driven convection. These produce plasmaspheric plumes, collisional heating and wave-particle interactions, centrifugal acceleration, and auroral acceleration by parallel electric fields, including enhanced ambipolar fields from electron heating by precipitation particles. Solar wind energy dissipation is concentrated by the geomagnetic field into auroral regions with an amplification factor of 10-100, enhancing heavy species plasma and gas escape from gravity, and providing more current carrying capacity. Internal plasmas thus enable electromagnetic driving via coupling to the plasma and neutral gas. We assess the importance of each of these processes in terms of local escape flux production as well as global outflow, and suggest methods for their implementation within multi-species global simulation codes. We conclude by assessing outstanding obstacles to this objective.

Moore, Thomas Earle↗
NASA NTRS2009

Physical Processes for Driving Ionospheric Outflows in Global Simulations

We review and assess the importance of processes thought to drive ionospheric outflows, linking them as appropriate to the solar wind and interplanetary magnetic field, and to the spatial and temporal distribution of their magnetospheric internal responses. These begin with the diffuse effects of photoionization and thermal equilibrium of the ionospheric topside, enhancing Jeans' escape, with ambipolar diffusion and acceleration. Auroral outflows begin with dayside reconnexion and resultant field-aligned currents and driven convection. These produce plasmaspheric plumes, collisional heating and wave-particle interactions, centrifugal acceleration, and auroral acceleration by parallel electric fields, including enhanced ambipolar fields from electron heating by precipitating particles. Observations and simulations show that solar wind energy dissipation into the atmosphere is concentrated by the geomagnetic field into auroral regions with an amplification factor of 10-100, enhancing heavy species plasma and gas escape from gravity, and providing more current carrying capacity. Internal plasmas thus enable electromagnetic driving via coupling to the plasma, neutral gas and by extension, the entire body " We assess the Importance of each of these processes in terms of local escape flux production as well as global outflow, and suggest methods for their implementation within multispecies global simulation codes. We complete 'he survey with an assessment of outstanding obstacles to this objective.

Moore, Thomas Earle↗
NASA NTRS2017

Simultaneous Observations of Traveling Convection Vortices: Ionosphere-Thermosphere Coupling

We present simultaneous observations of magnetosphere-ionosphere-thermosphere coupling over Svalbard during a traveling convection vortex (TCV) event. Various spaceborne and ground-based instruments made coordinated measurements, including magnetometers, particle detectors, an all-sky camera, European Incoherent Scatter (EISCAT) Svalbard Radar, Super Dual Auroral Radar Network (SuperDARN), and SCANning Doppler Imager (SCANDI). The instruments recorded TCVs associated with a sudden change in solar wind dynamic pressure. The data display typical features of TCVs including vortical ionospheric convection patterns seen by the ground magnetometers and SuperDARN radars and auroral precipitation near the cusp observed by the all-sky camera. Simultaneously, electron and ion temperature enhancements with corresponding density increase from soft precipitation are also observed by the EISCAT Svalbard Radar. The ground magnetometers also detected electromagnetic ion cyclotron waves at the approximate time of the TCV arrival. This implies that they were generated by a temperature anisotropy resulting from a compression on the dayside magnetosphere. SCANDI data show a divergence in thermospheric winds during the TCVs, presumably due to thermospheric heating associated with the current closure linked to a field-aligned current system generated by the TCVs. We conclude that solar wind pressure impulse-related transient phenomena can affect even the upper atmospheric dynamics via current systems established by a magnetosphere-ionosphere-thermosphere coupling process.

Kim, Hyomin↗
NASA NTRS2010

Lithosphere-Atmosphere-Ionosphere Coupling (LAIC) Model - An Unified Concept for Earthquake Precursors Validation

The paper presents a conception of complex multidisciplinary approach to the problem of clarification the nature of short-term earthquake precursors observed in atmosphere, atmospheric electricity and in ionosphere and magnetosphere. Our approach is based on the most fundamental principles of tectonics giving understanding that earthquake is an ultimate result of relative movement of tectonic plates and blocks of different sizes. Different kind of gases: methane, helium, hydrogen, and carbon dioxide leaking from the crust can serve as carrier gases for radon including underwater seismically active faults. Radon action on atmospheric gases is similar to the cosmic rays effects in upper layers of atmosphere: it is the air ionization and formation by ions the nucleus of water condensation. Condensation of water vapor is accompanied by the latent heat exhalation is the main cause for observing atmospheric thermal anomalies. Formation of large ion clusters changes the conductivity of boundary layer of atmosphere and parameters of the global electric circuit over the active tectonic faults. Variations of atmospheric electricity are the main source of ionospheric anomalies over seismically active areas. Lithosphere-Atmosphere-Ionosphere Coupling (LAIC) model can explain most of these events as a synergy between different ground surface, atmosphere and ionosphere processes and anomalous variations which are usually named as short-term earthquake precursors. A newly developed approach of Interdisciplinary Space-Terrestrial Framework (ISTF) can provide also a verification of these precursory processes in seismically active regions. The main outcome of this paper is the unified concept for systematic validation of different types of earthquake precursors united by physical basis in one common theory.

Pulinets, S.↗
NASA NTRS2013

Wavelength Dependence of Solar Irradiance Enhancement During X-class Flares and Its Influence on the Upper Atmosphere

The wavelength dependence of solar irradiance enhancement during flare events is one of the important factors in determining how the Thermosphere-Ionosphere (TI) system responds to flares. To investigate the wavelength dependence of flare enhancement, the Flare Irradiance Spectral Model (FISM) was run for 61X-class flares. The absolute and the percentage increases of solar irradiance at flare peaks, compared to pre-flare conditions, have clear wavelength dependences. The 0-4 nm irradiance increases much more ((is) approximately 680 on average) than that in the 14-25 nm waveband ((is) approximately 65 on average), except at 24 nm ( (is) approximately 220). The average percentage increases for the 25-105 nm and 122-190 nm wave bands are approximately 120 and approximately 35, respectively. The influence of 6 different wavebands (0-14 nm, 14-25 nm, 25-105 nm, 105-120 nm, 121.56 nm,and122-175 nm) on the thermosphere was examined for the October 28th, 2003 flare (X17-class) event by coupling FISM with the National Center for Atmospheric Research (NCAR) Thermosphere-Ionosphere-Electrodynamics General Circulation Model(TIE-GCM) under geomagnetically quiet conditions (Kp=1). While the enhancement in the0-14nm waveband caused the largest enhancement of the globally integrated solar heating, the impact of solar irradiance enhancement on the thermosphere at 400 km is largest for the 25-105 nm waveband (EUV), which accounts for about 33 K of the total 45 K temperature enhancement, and approximately 7.4% of the total approximately 11.5% neutral density enhancement. The effect of 122-175 nm flare radiation on the thermosphere is rather small. The study also illustrates that the high-altitude thermospheric response to the flare radiation at 0-175 nm is almost a linear combination of the responses to the individual wavebands. The upper thermospheric temperature and density enhancements peaked 3-5 h after the maximum flare radiation.

FISM↗
NASA NTRS1982

Stimulated plasma instability and nonlinear phenomena in the ionosphere

Several hundred topside ionograms were used to study simulated wave-particle interactions in the ionosphere. The study combined the benefits of high-frequency-resolution Alouette 2 analog sounder data with advanced digital graphics techniques. The study shows that the sounder phase can cause significant plasma heating when the plasma parameter is confined to specific ranges. The observations support the Harris instability generation process and the nonlinear Landau damping maintaining process for long-duration diffuse resonances. The observations also suggest that the so-called Q resonances have characteristics which imply that generation processes in a sounder-stimulated plasma turbulence may be involved.

Benson, R. F.↗
NASA NTRS1987

Orbital atmospheric physics and dynamics

There are two ways of modeling the upper atmosphere. One is the empirical model that makes use of experimental data on means and excursions from the mean and fits the data in a self-consistent manner. The other approach is to deal directly with the physical processes. This is difficult since what is happening is extremely complex. Data measured using an interferometer to give Doppler shifts of airglow lines showed 300 to 800 m/sec winds with a complex structure in the upper region of the thermosphere at high latitudes. Ionospheric electric fields, strongly influenced by interaction with the solar wind, drive the ionized component and large neutral winds result due to momentum transfer between the charged particles and the neutrals. Frictional heating results from movement of ions through the neutrals, which also influences the compositional structure. These are examples of the complex interactions involved. The NCAR General Circulation Model (tropospheric) was adapted for use at thermospheric altitudes: the Thermospheric General Circulation Model (TGCM). The model makes use partly of primitive equations and partly of empirical data for some quantities such as electron density, magnetic field, and ion drift.

Roble, Raymond↗