Search NASA⌕ Search

Engineering topics

Emery, B. A.

Publications and source records attributed to Emery, B. A..

At least 19 records

CEDAR-GEM Challenge for Systematic Assessment of Ionosphere/Thermosphere Models in Predicting TEC During the 2006 December Storm Event

In order to assess current modeling capability of reproducing storm impacts on total electron content (TEC), we considered quantities such as TEC, TEC changes compared to quiet time values, and the maximum value of the TEC and TEC changes during a storm. We compared the quantities obtained from ionospheric models against ground-based GPS TEC measurements during the 2006 AGU storm event (14-15 December 2006) in the selected eight longitude sectors. We used 15 simulations obtained from eight ionospheric models, including empirical, physics-based, coupled ionosphere-thermosphere, and data assimilation models. To quantitatively evaluate performance of the models in TEC prediction during the storm, we calculated skill scores such as RMS error, Normalized RMS error (NRMSE), ratio of the modeled to observed maximum increase (Yield), and the difference between the modeled peak time and observed peak time. Furthermore, to investigate latitudinal dependence of the performance of the models, the skill scores were calculated for five latitude regions. Our study shows that RMSE of TEC and TEC changes of the model simulations range from about 3 TECU (total electron content unit, 1 TECU = 1016 el m2) (in high latitudes) to about 13 TECU (in low latitudes), which is larger than latitudinal average GPS TEC error of about 2 TECU. Most model simulations predict TEC better than TEC changes in terms of NRMSE and the difference in peak time, while the opposite holds true in terms of Yield. Model performance strongly depends on the quantities considered, the type of metrics used, and the latitude considered.

Shim, J. S.↗

Atmospheric response in aurora experiment: Observations of E and F region neutral winds in a region of postmidnight diffuse aurora

The goal of the Atmospheric Response in Aurora (ARIA) experiment carried out at Poker Flat, Alaska, on March 3, 1992, was to determine the response of the neutral atmosphere to the long-lived, large-scale forcing that is characteristic of the diffuse aurora in the post midnight sector. A combination of chemical release rocket wind measurements, instrumented rocket composition measurements, and ground-based optical measurements were used to characterize the response of the neutral atmosphere. The rocket measurements were made at the end of a 90-min period of strong Joule heating. We focus on the neutral wind measurements made with the rocket. The forcing was determined by running the assimilated mapping of ionospheric electrodynamics (AMIE) analysis procedure developed at the National Center for Atmospheric Research. The winds expected at the latitude and longitude of the experiment were calculated using the spectral thermospheric general circulation model developed at the Danish Meteorological Institute. Comparisons of the observations and the model suggest that the neutral winds responded strongly in two height ranges. An eastward wind perturbation of approximately 100 m/s developed between 140 and 200 km altitude with a peak near 160 km. A southwestward wind with peak magnitude of approximately 150 m/s developed near 115 km altitude. The large amplitude winds at the lower altitude are particularly surprising. They appear to be associated with the upward propagating semidiurnal tide. However, the amplitude is much larger than predicted by any of the tidal models, and the shear found just below the peak in the winds was nominally unstable with a Richardson number of approximately 0.08.

Larsen, M. F.↗

One-dimensional hybrid satellite track model for the Dynamics Explorer 2 (DE 2) satellite

A one-dimensional hybrid satellite track model has been developed to calculate the high-latitude thermospheric/ionospheric structure below the satellite altitude using Dynamics Explorer 2 (DE 2) satellite measurements and theory. This model is based on Emery et al. (1985) satellite track code but also includes elements of Roble et al. (1987b) global mean thermosphere/ionosphere model. A number of parameterizations and data handling techniques are used to input satellite data from several DE 2 instruments into this model. Profiles of neutral atmospheric densities are determined from the Mass Spectrometer Incoherent Scatter 1990 (MSIS-90) model and measured neutral temperatures. Measured electron precipitation spectra are used in an auroral model to calculate particle impact ionization rates below the satellite. These rates are combined with a solar ionization rate profile and used to solve the O(+) diffusion equation, with the measured electron density as an upper boundary condition. The calculated O(+) density distribution, as well as the ionization profiles, are then used in a photochemical equilibrium model to calculate the electron and molecular ion densities. The electron temperature is also calculated by solving the electron energy equation with an upper boundary condition determined by the DE 2 measurement. The model enables calculations of altitude profiles of conductivity and Joule heating rate along and below the satellite track. In a first application of the new model, a study is made of thermospheric and ionospheric structure below the DE 2 satellite for a single orbit which occurred on October 25, 1981. The field-aligned Poynting flux, which is independently obtained for this orbit, is compared with the model predictions of the height-integrated energy conversion rate. Good quantitative agreement between these two estimates has been reached. In addition, measurements taken at the incoherent scatter radar site at Chatanika (65.1 deg N, 147.4 deg W) during a DE 2 overflight are compared with the model calculations. A good agreement was found in lower thermospheric conductivities and Joule heating rate.

Deng, Wei↗

One-Dimensional Hybrid Satellite Track Model for the Dynamics Explorer 2 (DE 2) Satellite

A one-dimensional hybrid satellite track model has been developed to calculate the high-latitude thermospheric/ionospheric structure below the satellite altitude using Dynamics Explorer 2 (DE 2) satellite measurements and theory. This model is based on Emery et al. satellite track code but also includes elements of Roble et al. global mean thermosphere/ionosphere model. A number of parameterizations and data handling techniques are used to input satellite data from several DE 2 instruments into this model. Profiles of neutral atmospheric densities are determined from the MSIS-90 model and measured neutral temperatures. Measured electron precipitation spectra are used in an auroral model to calculate particle impact ionization rates below the satellite. These rates are combined with a solar ionization rate profile and used to solve the O(+) diffusion equation, with the measured electron density as an upper boundary condition. The calculated O(+) density distribution, as well as the ionization profiles, are then used in a photochemical equilibrium model to calculate the electron and molecular ion densities. The electron temperature is also calculated by solving the electron energy equation with an upper boundary condition determined by the DE 2 measurement. The model enables calculations of altitude profiles of conductivity and Joule beating rate along and below the satellite track. In a first application of the new model, a study is made of thermospheric and ionospheric structure below the DE 2 satellite for a single orbit which occurred on October 25, 1981. The field-aligned Poynting flux, which is independently obtained for this orbit, is compared with the model predictions of the height-integrated energy conversion rate. Good quantitative agreement between these two estimates has been reached. In addition, measurements taken at the incoherent scatter radar site at Chatanika (65.1 deg N, 147.4 deg W) during a DE 2 overflight are compared with the model calculations. A good agreement was found in lower thermospheric conductivities and Joule heating rate.

Deng, Wei↗

Ground-based studies of ionospheric convection associated with substorm expansion

The instantaneous patterns of electric fields and currents in the high-latitude ionosphere are deduced by combining satellite and radar measurements of the ionospheric drift velocity, along with ground-based magnetometer observations for October 25, 1981. The period under study was characterized by a relatively stable southward interplanetary magnetic field (IMF), so that the obtained electric field patterns do reflect, in general, the state of sustained and enhanced plasma convection in the magnetosphere. During one of the satellite passes, however, an intense westward electrojet caused by a substorm intruded into the satellite (DE2) and radar (Chatanika, Alaska) field of view in the premidnight sector, providing a unique opportunity to differentiate the enhanced convection and substorm expansion fields. The distributions of the calculated electric potential for the expansion and maximum phases of the substorm show the first clear evidence of the coexistence of two physically different systems in the global convection pattern. The changes in the convection pattern during the substorm indicate that the large-scale potential distributions are indeed of general two-cell patterns representing the southward IMF status, but the night-morning cell has two positive peaks, one in the midnight sector and the other in the late morning hours, corresponding to the substorm expansion and the convection enhancement, respectively.

Kamide, Y.↗

Interhemispheric asymmetry of the high-latitude ionospheric convection pattern

The assimilative mapping of ionospheric electrodynamics technique has been used to derive the large-scale high-latitude ionospheric convection patterns simultaneously in both northern and southern hemispheres during the period of January 27 to 29, 1992. When the interplanetary magnetic field (IMF) B(sub z) component is negative, the convection patterns in the southern hemisphere are basically the mirror images of those in the northern hemisphere. The total cross-polar cap potential drops in the two hemispheres are similar. When B(sub z) is positive and absolute value of B(sub y) greater than B(sub z), the convection configurations are mainly determined by B(sub y) and they may appear as normal 'two-cell' patterns in both hemispheres much as one would expect under southward IMF conditions. However, there is a significant difference in the cross-polar-cap potential drop between the two hemispheres, with the potential drop in the southern (summer) hemisphere over 50% larger than that in the northern (winter) hemisphere. As the ratio of absolute value of B(sub y)/B(sub z) decreases (less than one), the convection configuration in the two hemispheres may be significantly different, with reverse convection in the southern hemisphere and weak but disturbed convection in the northern hemisphere. By comparing the convection patterns with the corresponding spectrograms of precipitating particles, we interpret the convection patterns in terms of the concept of merging cells, lobe cells, and viscous cells. Estimates of the ' merging cell' potential drops, that is, the potential ascribed to the opening of the dayside field lines, are usually comparable between the two hemispheres, as they should be. The 'lobe cell' provides a potential between 8.5 and 26 kV and can differ greatly between hemispheres, as predicted. Lobe cells can be significant even for southward IMF, if absolute value of B(sub y) greater than the absolute value of B(sub z). To estimate the potential drop of the 'viscous cells,' we assume that the low-latitude boundary layer is on closed field lines. We find that this potential drop varies from case to case, with a typical value of 10 kV. If the source of these cells is truly a viscous interaction at the flank of the magnetopause, the process is likely spatially and temporally varying rather than steady state.

Lu, G.↗

Mapping ionospheric convection response to IMF B(sub y) negative and B(sub z) positive conditions

We have used the Assimilative Mapping of Ionospheric Electrodynamics (AMIE) procedure to derive electric potential patterns for a period of strong and slowly varying northward Interplanetary Magnetic Field (IMF). Using these patterns we offer some insight into high-latitude convection response to differing IMF conditions between the hours of 0850 and 1720 UT on 14 January 1988. We introduce a newly assimilated data set, Defense Meteorological Satellite Program (DMSP) ion drift measurements, as a means of mapping electric fields in polar regions. Also, for the first time we map both southern and northern hemisphere convection flows. We show that significantly different convection patterns are simultaneously possible in the two polar caps when B(sub z) is positive and that asymmetries in these patterns are influenced by the strength of B(sub y) and B(sub x). When the ratio absolute value of B(sub y)/B(sub z) exceeds unity the polar cap convection increases and changes from sunward to anti-sunward for B(sub z) positive conditions.

Knipp, D. J.↗

Ionospheric convection response to slow, strong variations in a Northward interplanetary magnetic field: A case study for January 14, 1988

We analyze ionospheric convection patterns over the polar regions during the passage of an interplanetary magnetic cloud on January 14, 1988, when the interplanetary magnetic field (IMF) rotated slowly in direction and had a large amplitude. Using the assimilative mapping of ionospheric electrodynamics (AMIE) procedure, we combine simultaneous observations of ionspheric drifts and magnetic perturbations from many different instruments into consistent patterns of high-latitude electrodynamics, focusing on the period of northward IMF. By combining satellite data with ground-based observations, we have generated one of the most comprehensive data sets yet assembled and used it to produce convection maps for both hemispheres. We present evidence that a lobe convection cell was embedded within normal merging convection during a period when the IMF B(sub y) and B(sub z) components were large and positive. As the IMF became predominantly northward, a strong reversed convection pattern (afternoon-to-morning potential drop of around 100 kV) appeared in the southern (summer) polar cap, while convection in the northern (winter) hemisphere became weak and disordered with a dawn-to-dust potential drop of the order of 30 kV. These patterns persisted for about 3 hours, until the IMF rotated significantly toward the west. We interpret this behavior in terms of a recently proposed merging model for northward IMF under solstice conditions, for which lobe field lines from the hemisphere tilted toward the Sun (summer hemisphere) drape over the dayside magnetosphere, producing reverse convection in the summer hemisphere and impeding direct contact between the solar wind and field lines connected to the winter polar cap. The positive IMF B(sub x) component present at this time could have contributed to the observed hemispheric asymmetry. Reverse convection in the summer hemisphere broke down rapidly after the ratio absolute value of B(sub y)/B(sub z) exceeded unity, while convection in the winter hemisphere strengthened. A dominant dawn-to-dusk potential drop was established in both hemispheres when the magnitude of B(sub y) exceeded that of B(sub z) with potential drops of the order of 100 kV, even while B(sub z) remained northward. The latter transition to southward B(sub z) produced a gradual intensification of the convection, but a greater qualitative change occurred at the transition through absolute value of B(sub y)/B(sub z) = 1 than at at the transition through B(sub z) = 0. The various convection patterns we derive under northward IMF conditions illustrate all possibilities previously discussed in the literature: nearly single-cell and multicell, distorted and symmetric, ordered and unordered, and sunward and antisunward.

Knipp, D. J.↗

Electric potential patterns deduced for the SUNDIAL period of September 23-26, 1986

High-latitude electric potential and ionospheric conductance patterns are presented and discussed for the Northern Hemisphere during the SUNDIAL period of September 23-26, 1986 using the Assimilative Mapping of Ionospheric Electrodynamics (AMIE) technique of Richmond and Kamide (1988). Data sources used for the model are satellite and ground magnetometers, electron precipitation instruments, incoherent scatter radars, and ionospheric coherent radars. The period was characterized by many substorms, and a wide variety of instantaneous patterns of electrodynamic parameters on a hemispheric scale are derived throughout this period, of which only a few examples are displayed. The entire set of electric potentials and conductances are being made available through the NCAR CEDAR Data Base for further analysis and utilization in simulation models.

Emery, B. A.↗

Thermospheric dynamics during September 18-19, 1984. I - Model simulations

The September 18-19, 1984 Equinox Transient Study (ETS) interval was simulated using the National Center for Atmospheric Research thermospheric GCM with inputs guided by observations. High-latitude forcings for the model were defined with precision for this period. The importance of upward-propagating semidiurnal tides is emphasized, the first-realistic time-dependent simulation to incorporate such tides. A description of the simulated quiet time thermospheric structure for September 18 is contrasted with the storm time behavior of September 19. Several features of the storm response are described, including the generation of large-scale equatorward propagating disturbances. Longer-lived perturbations of the neutral temperature and composition are also discussed. A new feature discovered in the model predictions for altitudes around 200 km comprises a quasi-fixed, four-cell pattern of high-and low-density regions in the polar cap.

Crowley, G.↗

Thermospheric dynamics during September 18-19, 1984. II - Validation of the NCAR thermospheric general circulation model

The winds, temperatures, and densities predicted by the thermospheric GCM are compared with measurements from the Equinox Transition Study of September 17-24, 1984. Agreement between predictions and observation is good in many respects. The quiet day observations contain a strong semidiurnal wind variation which is mainly due to upward-propagating tides. The storm day wind behavior is significantly different and includes a surge of equatorward winds due to a global propagating disturbance associated with the storm onset. A quantitative statistical comparison of the predicted and measured winds indicates that the equatorward winds in the model are weaker than the observed winds, particularly during storm times. A quiet day phase anomaly in the measured F region winds which is not reproduced by the model suggests the occurrence of an important unmodeled interaction between upward propagating semidiurnal tides and high-latitude effects.

Crowley, G.↗

On the mechanisms responsible for high-latitude thermospheric composition variations during the recovery phase of a geomagnetic storm

The causal mechanisms for the recovery of the perturbed high-latitude thermospheric composition to the unperturbed state in the period following a geomagnetic storm are investigated. Model runs of the NCAR thermosphere/ionosphere GCM (TIGCM) and thermosphere GCM are used to calculate the averaged mass mixing ratio variations and the forcing terms responsible for these variations during the recovery phase of a geomagnetic storm. High latitude compositional recovery is found to occur in the NCAR TIGCM on a time scale of about 12 hr to 1 day. This time scale is in agreement with previously observed time scales for typical poststorm F region electron density recoveries. Neither molecular diffusion nor large-scale horizontal advection is the dominant process in determining the compositional state during the recovery period. Thermospheric compositional recovery at high geomagnetic latitudes is driven primarily by vertical advection.

Burns, A. G.↗

Joule heating in the mesosphere and thermosphere during the July 13, 1982, solar proton event

The solar proton event of July 13, 1982 produced considerable ionization in the polar-cap mesosphere. Energetic solar proton fluxes were measured by the NOAA-6 satellite. The DE-2 satellite measured the low-energy electrons, the ion drift velocity, and other atmospheric and ionospheric properties during the event in the region of the measured maximum electric field (189 mV/m at 2215 UT near 60 deg N), a Joule heating rate of 1-3 K/day is calculated between 70 and 80 km, exceeding the heating due to ozone absorption at noon in the summer hemisphere in that altitude range. The Joule heating rate above 90 km greatly exceeded 20 K/day. The calculated height-integrated Joule heating rate above 100 km in the same region exceeded 400 ergs/sq cm sec, and DE-2 near 350 km measured neutral winds of nearly 1000 m/s and neutral gas temperatures of over 2000 K. The overall ionospheric structure calculated below the DE-2 satellite is described.

Roble, R. G.↗

Millstone Hill measurements and TGCM simulation for the 30 May 1984 annular solar eclipse

On 30 May 1984, the Millstone Hill incoherent scatter radar was operated to gather data on the effects of the annular solar eclipse on the structure and dynamics of the ionosphere and thermosphere. The eclipse path was about 3 deg. south of Millstone which experienced a maximum obscuration of 86% at about 1705 UT. Both the zenith steerable antennas at Millstone were used in the experiment to collect data on the temporal evolution of the eclipse effects. This experiment represented the first opportunity at Millstone to collect data during an eclipse in the absence of a major magnetospheric disturbance which had previously made the unravelling of eclipse effects difficult. In addition, the configuration of the experiment and analysis of the data included a detailed examination of the effects on the neutral atmosphere. A major catalyst for this study was the opportunity to compare the results with the predictions made from the Thermospheric General Circulation Model (TGCM) at NCAR, as a calibration point for the model. The Arecibo and Sondrestrom radars also participated as part of a radar chain experiment. The analysis of the data from these stations is in progress; the initial results from Millstone Hill are presented.

Salah, J. E.↗

Thermospheric and ionospheric structure of the Southern Hemisphere polar cap on October 21, 1981, as determined from Dynamics Explorer 2 satellite data

For a number of years, satellites have been employed to measure auroral particles and fields within the high-latitude thermosphere. In the present paper, data from orbit 1174 of the Dynamics Explorer 2 satellite on October 21, 1981, are utilized to calculate the structure of the ionosphere and thermosphere below the satellite altitude down to about 80 km. Attention is given to details regarding the DE 2 measurements, a satellite track model, the calculated ionospheric structure, and the calculated neutral gas heating rates. The investigation demonstrates that the technique of deriving characteristics of the ionosphere and atmosphere below a satellite track promises to be very fruitful for defining the characteristics of the lower thermosphere for use in large numerical models of the thermosphere and ionosphere.

Emery, B. A.↗

Thermospheric circulation, temperature, and compositional structure of the southern hemisphere polar cap during October-November 1981

Results of NCAR thermospheric general circulation model (TGCM) simulations were compared with Dynamics Explorer (DE2) data on the neutral wind, temperature and O and N2 number variations over the south polar cap in the F region in October-November 1981. Attention was focused on variations in the thermospheric characteristics. All variations reached a maximum over the magnetic polar cap, with the largest O and N2 changes being present in the dawn sector of the auroral oval, where Joule heating was maximized. The TGCM simulations for conditions similar to the DE2 passes agreed well with the satellite data, indicating that the variations in the wind, temperature and compositional characteristics are dependent on the magnitude of magnetospheric convection and influenced by the separation between the geomagnetic and geographic poles.

Roble, R. G.↗

Observations of the dynamics of the polar thermosphere

The dynamics of the polar thermosphere are examined by using observations made from the Dynamics Explorer 2 satellite. The results used in this study were obtained primarily from the Fabry-Perot interferometer (FPI) and the wind and temperature spectrometer (WATS) during the time period from September 1981 through January 1982. Two primary geophysical conditions were examined: these were the southern summer and the northern winter polar regions. The results support the conclusion that above 60 degrees of latitude the neutral winds are strongly controlled by ion/neutral frictional momentum transfer resulting from magnetospheric convection. This implies that the natural coordinate system within which to display the neutral winds in the high polar thermosphere is magnetic. The collected observations of this study were used to assess the validity of two of the large thermospheric general circulation models. The result of this assessment was that the models reasonably represent the vector winds at high altitudes but do not, at present, accurately simulate the thermodynamics of that regime.

Hays, P. B.↗

The high latitude circulation and temperature structure of the thermosphere near solstice

NCAR thermospheric-general-circulation-model (TGCM) computations of solar-maximum thermospheric neutral-gas temperature and circulation around the December solstice are presented and discussed. The TGCM uses a 5 x 5-deg grid and 24 constant-pressure layers, corresponding to altitudes of about 97-500 km. The results are mapped as electron-density contours, polar plots, cylindrical equidistant projections, meridional cross sections, and F-region polar plots comparing the TGCM predictions with DE-2 satellite observations. The significant differences between summer and winter high-latitude F-region winds are attributed to the ion drag momentum associated with magnetospheric convection. The TGCM wind predictions follow the same pattern as the satellite measurements but are too small; possible model corrections are considered.

Roble, R. G.↗