Components, the common component architecture, and the climate/ocean/weather community
Explore the source record for details and available documents.
common component architecture↗
Engineering topics
Rasmussen, C. E.
Publications and source records attributed to Rasmussen, C. E..
Two-stream modeling of plasmaspheric refilling
Plasmaspheric refilling on an L = 4 flux tube was studied by using a time-dependent, hydrodynamic plasmaspheric flow model in which the ion streams from the two hemispheres are treated as distinct fluids. In the model the continuity, momentum, and energy equations of a two-ion (O(+) and H(+)), quasi-neutral, currentless plasma are solved along a closed geomagnetic field line; diffusive equilibrium is not assumed. collisions between all stream pairs and with neutral species are included. The model includes a corotating, tilted dipole magnetic field and neutral winds. Ionospheric sources and sinks are accounted for in a self-consistent manner. Electrons are assumed to be heated by photoelectrons. The model flux tube extends from a 200-km altitude in one hemisphere to a 200-km altitude in the other hemisphere. Initially, the upwelling streams pass through each other practically unimpeded. When the streams approach the boundary in the conjugate ionosphere, a shock develops there, which moves upward and dissipates slowly; at about the same time a reverse shock develops in the hemisphere of origin, which moves upward. After about 1 hour, large shocks develop in each stream near the equator; these shocks move toward the equator and downward after crossing the equator. However, these shocks are probably artificial, because counterstreaming flows occur in each H(+) fluid, which the model can only handle by creating shocks.
Guiter, S. M.↗
What is the source of observed annual variations in plasmaspheric density?
Plasmaspheric densities have been observed previously to be higher in December than in June, with the ratio varying between 1.5 and 3.0 and with larger variations at lower L shells. In order to search for the cause of the observed annual variations, we have modeled plasmaspheric density, using a time-dependent hydrodynamic model. On an L = 2 field line with geomagnetic longitude equal to 300 deg, the modeled plasmaspheric densities were a factor of 1.5 times higher in December than in June. The modeled December to June density ratio was found to increase slightly with L shell, in contrast to observations; this discrepancy may be due to the fact that outer plasmaspheric flux tubes are never completely full. In addition, for an L = 2 field line with geomagnetic longitude equal to 120 deg, the modeled plasmaspheric density was higher in June than in December by a factor of about 1.2. Various numerical tests were also performed in order to examine the sensitivity of plasmaspheric density to various parameters. In particular, a large vertical neutral wind was applied in order to raise the O(+) profile which had the effect of raising plasmaspheric density by a factor of 6. This in conjunction with a theoretical analysis suggests that plasmaspheric density levels are very sensitive to O(+) levels in the upper ionosphere. We conclude that annual variations in plasmaspheric density are due to similar variations in ionospheric O(+).
Guiter, S. M.↗
A prototype Upper Atmospheric Research Collaboratory (UARC)
The National Collaboratory concept has great potential for enabling 'critical mass' working groups and highly interdisciplinary research projects. We report here on a new program to build a prototype collaboratory using the Sondrestrom Upper Atmospheric Research Facility in Kangerlussuaq, Greenland and a group of associated scientists. The Upper Atmospheric Research Collaboratory (UARC) is a joint venture of researchers in upper atmospheric and space science, computer science, and behavioral science to develop a testbed for collaborative remote research. We define the 'collaboratory' as an advanced information technology environment which enables teams to work together over distance and time on a wide variety of intellectual tasks. It provides: (1) human-to-human communications using shared computer tools and work spaces; (2) group access and use of a network of information, data, and knowledge sources; and (3) remote access and control of instruments for data acquisition. The UARC testbed is being implemented to support a distributed community of space scientists so that they have network access to the remote instrument facility in Kangerlussuaq and are able to interact among geographically distributed locations. The goal is to enable them to use the UARC rather than physical travel to Greenland to conduct team research campaigns. Even on short notice through the collaboratory from their home institutions, participants will be able to meet together to operate a battery of remote interactive observations and to acquire, process, and interpret the data.
Clauer, C. R.↗
A bounce-averaged kinetic model of the ring current ion population
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.
Jordanova, V. K.↗
Effect of magnetospheric convection on thermal plasma in the inner magnetosphere
The effects of E x B convection on the distribution of plasma parameters in the inner magnetosphere have been examined. Analytical solutions describing density distributions along convective trajectories in the equatorial plane have been found. These solutions suggest the following dependence of plasma concentration n on the magnetic field B along convective trajectories: n proportional to B(exp alpha), where the parameter alpha varies between alpha is less than or = between 4/3 and 2. The alpha = 2 case corresponds to disregarding transport parallel to B. The lower bound alpha = 4/3 describes the situation when parallel transport dominates over convective motions perpendicular to B. A solution has also been obtained describing the effect of convection on diffusive equilibrium profiles in the dayside plasmasphere. In addition, ion temperature variations due to adiabatic effects associated with plasma convection have been analyzed in detail. Convective drifts lead to ion temperature anisotropies with the value and sign of the anisotropy contingent on density and temperature variations, local time, and the location of the convective trajectory with respect to the location of the plasmapause. We have also found that convective motions lead to a substantial exchange of energy between the ionosphere and the plasmasphere through electron heat fluxes.
Khazanov, G. V.↗
Interaction of ring current and radiation belt protons with ducted plasmaspheric hiss. 1: Diffusion coefficients and timescales
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.
Kozyra, J. U.↗
Decay of equatorial ring current ions and associated aeronomical consequences
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.
Fok, M.-C.↗
Plasma transport in the equatorial ionosphere during the great magnetic storm of March 1989
We have modeled plasma transport in the low-latitude and equatorial ionosphere during the great magnetic storm of March 1989. Our goal was to provide a consistent explanation for the DMSP (Defense Meteorological Satellite Program) observations of dramatic decreases in ion density and rapid ion drifts in the low latitude ionosphere over South America during the storm. The modeling effort supports the hypothesis that abnormally large upward drifts lifted F region plasma above the satellite's altitude and created the density depletions observed by DMSP. Modeled O(+) densities at the satellite's altitude have a strong qualitative resemblance to DMSP observations. Both the model and the observations indicate a deep density trough with extremely sharp boundaries surrounding the equator. The widths of both the modeled and the observed equatorial troughs increase with time. Vertical ion drifts predicted by the model also have been compared with DMSP measurements. Like the observed vertical drifts, the modeled drifts reversed sign near the trough boundaries. The modeled vertical drifts are of the same order and direction as the vertical component of E x B convection near the equator, but of opposite direction (downward) near the trough boundaries and outside of the trough.
Rasmussen, C. E.↗
Anisotropic ion heating and parallel O(+) acceleration in regions of rapid E x B convection
A numerical solution to the 20-moment set of transport equations has been found in order to study subauroral ionospheric outflows during periods of enhanced perpendicular ion drifts. The numerical model solves the time-dependent O(+) density, momentum, and both the parallel and perpendicular energy and heat flow equations in the 200-6000 km altitude range. Assuming perpendicular drifts of 3 km/s relative to the neutral atmosphere, we have found that anisotropic heating of O(+) (a result of ion-neutral collisions) leads to a temperature anisotropy, with perpendicular temperatures exceeding 8000 K and parallel temperatures greater than 5000 K (near 200 km altitude). Above approximately 2000 km, transport processes dominate the effects of collisions and wavelike oscillations in O(+) velocity, temperature and heat flux were noted.
Korosmezey, A.↗
Diurnal variations on a plasmaspheric flux tube - Light ion flows and F region temperature enhancements
The paper concentrates on the diurnal variations on a plasmaspheric flux tube modeled using a time-dependent multispecies one-stream interhemispheric model for plasma flows. The model takes into account the effects of ionization, charge exchange, recombination, collisions, heat conduction, and allows for external heat sources. The simulation is done for June solstice conditions during solar minimum. Focus is placed on the presence of large downward H(+) velocities at about 320-km altitude in the winter (southern) hemisphere, in early morning when the summer hemisphere is sunlit but the winter hemisphere is dark. In addition, an upward H(+) flux is seen in the Southern Hemisphere at altitudes above 2000 km when the sun rises in the northern end.
Guiter, S. M.↗
The flow of plasma in the solar terrestrial environment
The overall goal of our NASA Theory Program was to study the coupling, time delays, and feedback mechanisms between the various regions of the solar-terrestrial system in a self-consistent, quantitative manner. To accomplish this goal, it will eventually be necessary to have time-dependent macroscopic models of the different regions of the solar-terrestrial system and we are continually working toward this goal. However, with the funding from this NASA program, we concentrated on the near-earth plasma environment, including the ionosphere, the plasmasphere, and the polar wind. In this area, we developed unique global models that allowed us to study the coupling between the different regions. These results are highlighted in the next section. Another important aspect of our NASA Theory Program concerned the effect that localized 'structure' had on the macroscopic flow in the ionosphere, plasmasphere, thermosphere, and polar wind. The localized structure can be created by structured magnetospheric inputs (i.e., structured plasma convection, particle precipitation or Birkland current patterns) or time variations in these input due to storms and substorms. Also, some of the plasma flows that we predicted with our macroscopic models could be unstable, and another one of our goals was to examine the stability of our predicted flows. Because time-dependent, three-dimensional numerical models of the solar-terrestrial environment generally require extensive computer resources, they are usually based on relatively simple mathematical formulations (i.e., simple MHD or hydrodynamic formulations). Therefore, another goal of our NASA Theory Program was to study the conditions under which various mathematical formulations can be applied to specific solar-terrestrial regions. This could involve a detailed comparison of kinetic, semi-kinetic, and hydrodynamic predictions for a given polar wind scenario or it could involve the comparison of a small-scale particle-in-cell (PIC) simulation of a plasma expansion event with a similar macroscopic expansion event. The different mathematical formulations have different strengths and weaknesses and a careful comparison of model predictions for similar geophysical situations provides insight into when the various models can be used with confidence.
Schunk, Robert W.↗
A three-dimensional time-dependent model of the plasmasphere
A three-dimensional, time-dependent, nonlinear, hydrodynamic model of the plasmasphere has been developed which includes the self-consistent coupling of conjugate hemispheres within the plasmasphere as well as the effects of cross-L drifts due to convection electric fields. The inner region of the plasmasphere was modeled from an L of 1.5 to an L of 3.5 to 4.5, depending on local time. In this first study of the global plasmasphere, the model was run for solar minimum conditions until diurnally reproducible results were obtained, indicating that the plasmasphere was fully filled. The results of the model were then compared with applicable measurements of the plasmasphere. It was found that the model tends to overestimate densities somewhat, although it is within a factor of 2 of an average of June and December whistler observations. The model predicts densities within a factor of 2 of most of the satellite observations as well. As anticipated, the diurnal variation of the plasmasphere was found to depend on magnetic latitude. At low latitudes the diurnal variation in density was relatively small, with the largest densities occurring in the afternoon time sector. However, near the plasmapause, the effects of changes in volume of drifting tubes of plasma due to cross-L drifts led to a factor of 3 variation in equatorial density, with the highest densities occurring near local midnight where the volume was lowest. Various limitations of the current model and suggested improvements are also discussed.
Rasmussen, C. E.↗
Multistream hydrodynamic modeling of interhemispheric plasma flow
Interhemispheric plasma flow was simulated using one-stream and two-stream hydrodymic models in order to test the suggestion of Banks et al. (1971) and others that the collision of high-speed flows originating from the conjugate hemispheres will cause the formation of a pair of shocks. The single-fluid hydrodynamic equations were modified to include multiple ion streams, allowing for the possibility of counterstreaming flow. It was found that a counterstreaming of ion streams from conjugate hemispheres does occur during the early stages of the refilling of plamaspheric flux tubes, and that a pair of reverse shocks does form. These shocks form away from the equator, and their subsequent motion creates conditions similar to those predicted by the single-stream hydrodynamic models. The findings support the conclusion of earlier studies that the refilling of the plasmasphere occurs from the equatorial region downward.
Rasmussen, C. E.↗
A photochemical equilibrium model for ionospheric conductivity
A photochemical equilibrium model of the high-latitude ionosphere has been developed. This model provides densities of the ionospheric constituents, N2(+), O2(+), O(+), and NO(+), from 85 km to approximately 220 km. These densities are then used to calculate Pedersen and Hall conductivities. A comparison of the model results with Arecibo and Chatanika radar observations was made, covering periods of solar minimum and solar maximum. The comparison showed the model to predict ionospheric densities to within 50 percent and conductivities to within 40 percent in the illuminated portion of the ionosphere. In regions of electron precipitation, the model showed good agreement with measurements. Results of this study indicate the following: (1) ionospheric conductivity increases by a factor of about 1.6 from solar minimum to solar maximum conditions, (2) the portion of the ionosphere above 170 km can contribute as much as 40 percent during daylight and 80 percent during nighttime to the total height-integrated Pedersen conductivity, and (3) the ratio of the height-integrated Hall to Pedersen conductivities is approximately 1.1-1.3 for sunlit conditions; this is appreciably lower than the value of 2 found in previous studies. These and other factors indicate that, under certain conditions, the height-integrated Pedersen conductivity may be as much as 2-3 times larger than previously reported.
Rasmussen, C. E.↗
Ionospheric convection inferred from interplanetary magnetic field-dependent Birkeland currents
Computer simulations of ionospheric convection have been performed, combining empirical models of Birkeland currents with a model of ionospheric conductivity in order to investigate IMF-dependent convection characteristics. Birkeland currents representing conditions in the northern polar cap of the negative IMF By component are used. Two possibilities are considered: (1) the morning cell shifting into the polar cap as the IMF turns northward, and this cell and a distorted evening cell providing for sunward flow in the polar cap; and (2) the existence of a three-cell pattern when the IMF is strongly northward.
Rasmussen, C. E.↗
Can the high latitude ionosphere support large field-aligned ion drifts?
A three-dimensional time-dependent model of the ionosphere is used to examine recent results on vertical drift velocities and electron densities in the high latitude ionosphere. Upper limits for the downward ion velocity were found to be smaller than those obtained from previous measurements. The downward force in the model was arbitrarily increased to try to account for the narrow electron density profiles. It is noted that if the common volume measurement is made in a region of O(+) precipitation then the line profile would not be Doppler shifted when viewed off-zenith, and small field-aligned velocities and narrow profile widths would result.
Sica, R. J.↗
Ionospheric convection driven by NBZ currents
Computer simulations of Birkeland currents and electric fields in the polar ionosphere during periods of northward IMF were conducted. When the IMF z component is northward, an additional current system, called the NBZ current system, is present in the polar cap. These simulations show the effect of the addition of NBZ currents on ionospheric convection, particularly in the polar cap. When the total current in the NBZ system is roughly 25 to 50 percent of the net region 1 and 2 currents, convection in the central portion of the polar cap reverses direction and turns sunward. This creates a pattern of four-cell convection with two small cells located in the polar cap, rotating in an opposite direction from the larger cells. When the Birkeland currents are fixed (constant current source), the electric field is reduced in regions of relatively high conductivity, which affects the pattern of ionospheric convection. Day-night asymmetries in conductivity change convection in such a way that the two polar-cap cells are located within the large dusk cell. When ionospheric convection is fixed (constant voltage source), Birkeland currents are increased in regions of relatively high conductivity. Ionospheric currents, which flow horizontally to close the Birkeland currents, are changed appreciably by the NBZ current system. The principal effect is an increase in ionospheric current in the polar cap.
Rasmussen, C. E.↗