Search NASA⌕ Search

Engineering topics

Fuller-Rowell, T. J.

Publications and source records attributed to Fuller-Rowell, T. J..

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.↗

Response of the thermosphere and ionosphere to geomagnetic storms

Four numerical simulations have been performed, at equinox, using a coupled thermosphere-ionosphere model, to illustrate the response of the upper atmosphere to geomagnetic storms. The storms are characterized by an increase in magnetospheric energy input at high latitude for a 12-hour period; each storm commences at a different universal time (UT). The initial response at high latitude is that Joule heating raises the temperature of the upper thermosphere and ion drag drives high-velocity neutral winds. The heat source drives a global wind surge, from both polar regions, which propagates to low latitudes and into the opposite hemisphere. The surge has the character of a large-scale gravity wave with a phase speed of about 600 m/s. Behind the surge a global circulation of magnitude 100 m/s is established at middle latitudes, indicating that the wave and the onset of global circulation are manifestations of the same phenomena. A dominant feature of the response is the penetration of the surge into the opposite hemisphere where it drives poleward winds for a few hours. The global wind surge has a preference for the night sector and for the longitude of the magnetic pole and therefore depends on the UT start time of the storm. A second phase of the meridional circulation develops after the wave interaction but is also restricted, in this case by the buildup of zonal winds via the Coriolis interaction. Conservation of angular momentum may limit the buildup of zonal wind in extreme cases. The divergent wind field drives upwelling and composition change on both height and pressure surfaces. The composition bulge responds to both the background and the storm-induced horizontal winds; it does not simply rotate with Earth. During the storm the disturbance wind modulates the location of the bulge; during the recovery the background winds induce a diurnal variation in its position. Equatorward winds in sunlight produce positive ionospheric changes during the main driving phase of the storm. Negative ionospheric phases are caused by increases of molecular nitrogen in regions of sunlight, the strength of which depends on longitude and the local time of the sector during the storm input. Regions of positive phase in the ionosphere persist in the recovery period due to decreases in mean molecular mass in regions of previous downwelling. Ion density changes, expressed as a ratio of disturbed to quiet values, exhibit a diurnal variation that is driven by the location of the composition bulge; this variation explains the ac component of the local time variation of the observed negative storm phase.

Fuller-Rowell, T. J.↗

Neutral winds in the lower thermosphere from Dynamics Explorer 2

Doppler line profile measurements of the OI lambda 557.7 nm 'green line' emission, made by the Fabry-Perot interferometer on Dynamics Explorer 2, have provided altitude profiles of the meridional component of the lower thermospheric neutral wind. The wind inversion technique of Nardi has been used to extract the neutral wind profiles from the line-of-sight measurements. Individual lambda 557.7 nm Doppler line profiles and inverted volume-emission-rate and neutral-wind profiles are presented. Neutral wind measurements from the about 120 km altitude level, obtained on multiple orbital passes over the summer hemisphere polar region, have been merged to produce a synthesized averaged 'vector' wind field in geomagnetic coordinates. The wind pattern exhibits a region of anticyclonic vorticity in the daytime sector of the magnetic polar cap. Averaged winds of about 300 m/sec in the equatorial direction are observed in the early morning sector. The satellite winds are in reasonable agreement with the predictions of the NCAR and UCL thermosphere general circulation models, with significant regional discrepancies evident in both magnitude and direction.

Killeen, T. L.↗

Comparison of the USU ionospheric model with the UCL-Sheffield coupled thermospheric-ionospheric model

Several physical models of the high-latitude ionosphere have been developed that describe the time-dependent evolution of the E- and F-region plasma density. The models require a variety of inputs, including solar EUV fluxes, magnetospheric convection, auroral precipitation, and neutral atmosphere. Of specific relevance to this study is how the neutral atmosphere is incorporated into the ionospheric models. For the USU ionospheric model, the neutral atmosphere is the MSIS 1986 empirical model, while for the UCL-Sheffield coupled thermospheric-ionospheric model the neutral atmosphere is computed simultaneously with the ionosphere. Both models were run for similar solar and magnetospheric conditions (solar maximum, moderate geomagnetic activity, and winter solstice). Solar maximum conditions ensured a strong coupling between the ionosphere and thermosphere, which provided the possibility of a large ionospheric difference between the two physical models. This was further enhanced by choosing winter conditions so that the densities were not dominated by sunlight. The comparison of the two models indicated that both models predict the same morphological features with similar ionospheric densities, generally within about 30 percent.

Sojka, J. J.↗

Modelling of E-region auroral winds

A numerical code, describing the interactions between chemistry and dynamics in the coupled thermosphere and ionosphere systems, has been used to investigate the morphology of wind systems of the auroral E-region and, in particular, their response to geomagnetic forcing. The structure and dynamics of both the thermosphere and the ionosphere show very strong interactions at all levels of geomagnetic disturbance. By comparison, other perturbations, particularly those resulting from tidal and gravity wave sources within the lower atmosphere, which are important and easily identified during quiet periods, are always overwhelmed by the magnitude of the effects of large geomagnetic disturbances above about 110 km altitude. The model simulations will be used to characterize the geomagnetic response of the E-region neutral wind, temperature, and composition in different parts of the polar cap and the auroral oval, and the major changes resulting from seasonal variations.

Rees, David↗

Numerical simulations of the seasonal/latitudinal variations of atomic oxygen and nitric oxide in the lower thermosphere and mesosphere

A 2-Dimensional zonally-averaged thermospheric model and the global University College London (UCL) thermospheric model have been used to investigate the seasonal, solar activity and geomagnetic variation of atomic oxygen and nitric oxide. The 2-dimensional model includes detailed oxygen and nitrogen chemistry, with appropriate completion of the energy equation, by adding the thermal infrared cooling by O and NO. This solution includes solar and auroral production of odd nitrogen compounds and metastable species. This model has been used for three investigations; firstly, to study the interactions between atmospheric dynamics and minor species transport and density; secondly, to examine the seasonal variations of atomic oxygen and nitric oxide within the upper mesosphere and thermosphere and their response to solar and geomagnetic activity variations; thirdly, to study the factor of 7 to 8 peak nitric oxide density increase as solar F sub 10.7 cm flux increases from 70 to 240 reported from the Solar Mesospheric Explorer. Auroral production of NO is shown to be the dominant source at high latitudes, generating peak NO densities a factor of 10 greater than typical number densities at low latitudes. At low latitudes, the predicted variation of the peak NO density, near 110 km, with the solar F sub 10.7 cm flux is rather smaller than is observed. This is most likely due to an overestimate of the soft X-ray flux at low solar activity, for times of extremely low support number, as occurred in June 1986. As observed on pressure levels, the variation of O density is small. The global circulation during solstice and periods of elevated geomagnetic activity causes depletion of O in regions of upwelling, and enhancements in regions of downwelling.

Rees, D.↗

A theoretical and empirical study of the response of the high latitude thermosphere to the sense of the 'Y' component of the interplanetary magnetic field

Patterns of magnetospheric energetic plasma precipitation as a function of the Y component of the Interplanetary Magnetic Field (IMF) are studied. The development of a three-dimensional, time-dependent global thermospheric model using a polar conversion electric field with a dependence on the Y component of the IMF to evaluate thermospheric wind circulation is examined. Thermospheric wind data from the ISEE-3 satellite, Dynamics Explorer-2 satellite, and a ground-based Fabry-Perot interferometer in Kiruna, Sweden, collected on December 1, 2, 6, 25, 1981 and February 12, 13, 1982 are described. The observed data and simulations of polar thermospheric winds are compared. In the Northern Hemisphere a strong antisunward ion flow on the dawn side of the geomagnetic polar cap is observed when the BY is positive, and the flow is detected on the dusk side when the BY is negative. It is concluded that the strength and direction of the IMF directly control the transfer of solar wind momentum and energy to the high latitude thermosphere.

Rees, D.↗

The composition, structure, temperature and dynamics of the upper thermosphere in the polar regions during October to December 1981

Observational data obtained by the Dynamics Explorer-2 (DE-2) spacecraft were compared with global model simulations in order to study the composition, structure, temperature, and dynamics of the upper thermosphere in polar regions during the period October-December 1981. A UCL three-dimensional model was used to simulate the seasonal, diurnal, and geomagnetic response of the neutral thermosphere and to follow the major features of the solar and geomagnetic inputs during the late 1981 period. Overall agreement was obtained between the simulations and the DE-2 data for thermospheric wind structure at high latitudes, and for the combined thermal and compositional structure in both hemispheres. Computer-generated line drawings of the variations in thermospheric structure are given, as well as a series of color graphic illustrations of the DE-2 data.

Rees, D.↗

The westward thermospheric jet-stream of the evening auroral oval

Ground-based and satellite measurements of the thermospheric wind in jet-streams during the evening auroral oval are analyzed, in order to study the geophysical mechanisms of thermospheric wind generation. Numerical simulations using a global, three-dimensional, time-dependent model of thermospheric dynamics were compared with the satellite data, and the results are discussed in detail. The wind distribution during the storm is shown in a series of color plates.

Rees, D.↗

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 effect of realistic conductivities on the high-latitude neutral thermospheric circulation

The dynamics of the high latitude thermosphere are dominated by the ion circulation pattern driven by magnetospheric convection. The reaction of the neutral thermosphere is influenced by both the magnitude of the ion convection velocity and by the conductivity of the thermosphere. Using a three-dimensional, time-dependent, thermospheric, neutral model together with different ionospheric models, the effect of changes in conductivity can be assessed. The ion density is described by two models: the first is the empirical model of Chiu (1975) appropriate for very quiet geomagnetic conditions, and the second is a modified version of the theoretical model of Quegan et al. (1982). The differences in the neutral circulation resulting from the use of these two ionospheric models emphasizes the need for realistic high latitude conductivities when attempting to model average or disturbed geomagnetic conditions, and a requirement that models should couple realistically the ionosphere and the neutral thermosphere. An attempt is made to qualitatively interpret some of the features of the neutral circulation produced at high latitudes by magnetospheric processes.

Fuller-Rowell, T. J.↗

A comparison of wind observations of the upper thermosphere from the Dynamics Explorer satellite with the predictions of a global time-dependent model

Experimental and theoretical investigations of high-latitude circulation in the upper thermosphere are reported. Vector-wind plots constructed by combining remotely sensed meridional-wind and in situ zonal-wind data obtained during four south-polar and three north-polar passes of the NASA Dynamics Explorer satellite in October and December, 1981, are presented and analyzed, taking the geomagnetic activity level and local solar time into account. The results are then compared with the predictions of 3D time-dependent global models of thermospheric neutral winds (Fuller-Rowell and Rees, 1980, 1981, 1983) adjusted to account for both solar-UV/EUV heating and ionization effects and quiet (Kp from 1 to 2) or moderately disturbed (Kp from 3 to 4) geomagnetic conditions (models Q and MD). An MD model incorporating a self-consistent description of the high-latitude ionosphere and a Q model excluding high-latitude effects are found to give the most accurate predictions for the respective geomagnetic conditions.

Rees, D.↗

Observations of neutral wind profiles between 115- and 175-km altitude in the dayside auroral oval

Lithium trail neutral wind measurements have been made in the morning dayside auroral oval as a function of altitude. During all four experiments the winds displayed a large altitude shear. In the most extreme case the wind at 135-km altitude was toward the south at 210 m/s, and at 185-km altitude toward the northwest at 170 m/s. A one-dimensional numerical solution to the Navier Stokes equation including an initial steady-state wind profile, viscosity, pressure gradient, Lorenz force and the Coriolis force was made in order to study the effect of auroral processes on the neutral atmosphere. Since electric fields and plasma density were measured, the Lorenz force could be modeled with reasonable accuracy. The pressure gradient was determined in an iterative way as a best fit to the observations. Within the limitations of a one-dimensional model the numerical solution is in reasonable agreement with one of the observations, which is applicable to such a modeling, and indicates that both Lorenz forces and pressure gradients due to auroral processes are required to explain the observation.

Pereira, E.↗