Engineering PapersSearch

Engineering topics

Roble, R. G.

Publications and source records attributed to Roble, R. G..

At least 19 records

Northern Hemisphere Atmospheric Influence of the Solar Proton Events and Ground Level Enhancement in January 2005

Solar eruptions in early 2005 led substantial barrage of charged particles on the Earth's atmosphere during the January 16-21 period. Proton fluxes were greatly increased during these several days and led to the production ofHO(x)(H, OH, BO2)and NO(x)(N, NO, NO2), which then caused the destruction of ozone. We focus on the Northern polar region, where satellite measurements and simulations with the Whole Atmosphere Community Climate Model (WACCM3) showed large enhancements in mesospheric HO(x) and NO(x) constituents, and associated ozone reductions, due 10 these solar proton events (SPEs). The WACCM3 simulations show enhanced short-lived OH throughout the mesosphere in the 60-82.5degN latitude band due to the SPEs for most days in the Jan.16-2l,2005 period, in reasonable agreement with the Aura Microwave Limb Sounder (MLS) measurements. Mesospheric HO2 is also predicted to be increased by the SPEs, however, the modeled HO2 results are somewhat larger than the MLS measurements. These HO(x) enhancements led to huge predicted and MLS-measured ozone decreases of greater than 40% throughout most of the Northern polar mesosphere during the SPE period. Envisat Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) measurements of hydrogen peroxide (H2O2) show increases throughout the stratosphere with highest enhancements of about 60 ppt y in the lowermost mesosphere over the Jan. 16-18, 2005 period due to the solar protons. WACCM3 predictions indicate H2O2 enhancements over the same time period of more than twice that amount. Measurements of nitric acid (HNO3) by both MLS and MIPAS show an increase of about 1 ppbv above background levels in the upper stratosphere during January 16-29, 2005. WACCM3 simulations show only minuscule HNO3 changes in the upper stratosphere during this time period. However due to the small loss rates during winter, polar mesospheric enhancements of NO(x) are computed to be greater than 50 ppbv during the SPE period. Computed NO(x)increases, which were statistically significant at the 95% level, lasted about a month past the SPEs. The SCISAT-I Atmospheric Chemistry Experiment Fourier Transform Spectrometer (ACE-FTS) NO(x) measurements and MIPAS NO, measurements for the polar Northern Hemisphere are in reasonable agreement with these predictions. An extremely large ground level enhancement (GLE) occurred during the SPE period on January 20, 2005. We find that protons of energies 300 to 20,000 MeV, not normally included in our computations, led to enhanced lower stratospheric odd nitrogen concentrations of less than 0.1% as a result of this GLE.

Jackman, C. H.

Geomagnetic Storm Effects in the Low- to Middle-Latitude Upper Thermosphere

In this paper, we use data from the Dynamics Explorer 2 (DE 2) satellite and a theoretical simulation made by using the National Center for Atmospheric Research thermosphere/ionosphere general circulation model (NCAR-TIGCM) to study storm-induced changes in the structure of the upper thermosphere in the low- to middle-latitude (20 deg-40 deg N) region of the winter hemisphere. Our principal results are as follows: (1) The winds associated with the diurnal tide weaken during geomagnetic storms, causing primarily zonally oriented changes in the evening sector, few changes in the middle of the afternoon, a combination of zonal and meridional changes in the late morning region, and mainly meridional changes early in the morning; (2) Decreases in the magnitudes of the horizontal winds associated with the diurnal tide lead to a net downward tendency in the vertical winds blowing through a constant pressure surface; (3) Because of these changes in the vertical wind, there is an increase in compressional heating (or a decrease in cooling through expansion), and thus temperatures in the low- to middle-latitudes of the winter hemisphere increase; (4) Densities of all neutral species increase on a constant height surface, but the pattern of changes in the O/N2 ratio is not well ordered on these surfaces; (5) The pattern of changes in the O/N2 ratio is better ordered on constant pressure surfaces. The increases in this ratio on constant pressure surfaces in the low- to middle-latitude, winter hemisphere are caused by a more downward tendency in the vertical winds that blow through the constant pressure surfaces. Nitrogen-poor air is then advected downward through the pressure surface, increasing the O/N2 ratio; (6) The daytime geographical distribution of the modeled increases in the O/N2 ratio on a constant pressure surface in the low- to middle-latitudes of the winter hemisphere correspond very closely with those of increases in the modeled electron densities at the F2 peak.

Burns, A. G.

Large Enhancements in the O/N2 Ratio in the Evening Sector of the Winter Hemisphere During Geomagnetic Storms

In this paper, we have looked for enhancements of the O/N2 ratio in data measured by the Dynamics Explorer 2 (DE 2) satellite in the middle latitudes of the winter hemisphere, based on a prediction that was made by the National Center for Atmospheric Research thermosphere/tonosphere general circulation model (NCAR-TIGCM) that such increases occur. The NCAR-TIGCM predicts that these enhancements should be seen throughout the low latitude region and in many middle latitude locations, but that the enhancements in O/N2 are particularly strong in the middle-latitude, evening-to-midnight sector of the winter hemisphere. When this prediction was used to look for these effects in DE 2 NACS (neutral atmosphere composition spectrometer) data, large enhancements in the O/N2 ratio (approx. 50 to 90%) were seen. These enhancements were observed during the main phase of a storm that occurred on November 24, 1982, and were seen in the same region of the winter hemisphere predicted by the NCAR-TIGCM. They are partially the result of the depletion of N2 and, as electron loss is dependent on dissociative recombination at F(sub 2) altitudes, they have implications for electron densities in this area. Parcel trajectories, which have been followed through the NCAR-TIGCM history file for this event, show that large O/N2 enhancements occur in this limited region in the winter hemisphere for two reasons. First, these parcels of air are decelerated by the antisunward edge of the ion convection pattern; individual parcels converge and subsidence occurs. Thus molecular-nitrogen-poor air is brought from higher to lower heights. Because neutral parcels that are found a little poleward of the equatorial edge of the eveningside convection pattern are swept inward toward the center of the auroral oval, the enhancements occur only in a very limited range of latitudes. Second, nitrogen-poor air is transported from regions close to the magnetic pole in the winter hemisphere. During geomagnetic storms, enhanced meridional winds are driven by the increased pressure-gradient force that is associated with intensified Joule heating in the auroral oval. These pressure-driven winds decrease rapidly on the dayside beyond the auroral oval where the parcels originate, limiting the region into which the parcels can be transported. Thus these two processes drive values of O/N2 in a limited region of the winter hemisphere, and reinforce only in the evening sector, causing large changes in this region.

Burns, A. G.

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

Database of theoretical thermosphere/ionosphere predictions

The National Center for Atmospheric Research (NCAR) thermosphere-ionosphere general circulation model has performed a series of self-consistent calculations of the thermosphere/ionosphere structure for a variety of geophysical conditions. The model simulates March, June, and December periods for solar cycle minimum and maximum for three levels of geomagnetic activity representing quiet, moderate, and disturbed conditions. The calculations include predictions of N2, O2, O, N(sup 4)S, N(sup 2)D, NO, N2(+), O2(+), O(+), NO(+), and N(+), along with the ion, electron, and neutral temperatures and three components of the neutral winds. The model simulations are available from the authors or through the Coupling, Energetic, and Dynamics of Atmospheric Regions database at NCAR.

Fesen, C. G.

CO2 cooling in terrestrial planet thermospheres

We examine the recent progress in the debate on the CO2-O relaxation rate, its temperature dependence, and its corresponding impact on the thermospheric heat budgets of Venus, Earth, and Mars. This comparative approach provides the broadest range of conditions under which a common CO2-O relaxation rate should provide consistent results. New global mean calculations are presented for the heat budgets of these three planets using large CO2-O relaxation rates that have been inferred recently from Earth CO2 radiance measurements and laboratory studies. Results indicate that available Venus and Mars data constrain the CO2-O relaxation rate to be 2-4 x 10(exp -12)/cu cm/s at 300 K. For Venus, this strong cooling serves as an effective thermostat that gives rise to a small variation of thermospheric temperatures over the solar cycle, just as observed. Conversely, CO2 cooling does not appear to be dominant in the dayside heat budget of the Mars thermosphere over most of the solar cycle. For the Earth, this strong cooling implies that the lower thermosphere does not typically require significant eddy diffusion or heat conduction. However, global-scale dynamics or an additional heating mechanism may be needed to restore calculated temperatures to observed values when relaxation rates exceeding 2 x 10(exp -12)/cu cm/s are employed.

Bougher, S. W.

A thermosphere-ionosphere-mesosphere-electrodynamic general circulation model (time-GCM): Equinox solar cycle minimum simulations (30-500 km)

A new simulation model of the mesosphere, thermosphere, and ionosphere with coupled electrodynamics has been developed and used to calculate the global circulation, temperature and compositional structure between 30-500 km for equinox, solar cycle minimum, geomagnetic quiet conditions. The model incorporates all of the features of the National Center for Atmospheric Research (NCAR) thermosphere-ionosphere- electrodynamics general circulation model (TIE-GCM) but the lower boundary has been extended downward from 97 to 30 km (10 mb) and it includes the physical and chemical processes appropriate for the mesosphere and upper stratosphere. The first simulation used Rayleigh friction to represent gravity wave drag in the middle atmosphere and although it was able to close the mesospheric jets it severely damped the diurnal tide. Reduced Rayleigh friction allowed the tide to penetrate to thermospheric heights but did not close the jets. A gravity wave parameterization developed by Fritts and Lu (1993) allows both features to exist simultaneously with the structure of tides and mean flow dependent upon the strength of the gravity wave source. The model calculates a changing dynamic structure with the mean flow and diurnal tide dominant in the mesosphere, the in-situ generated semi-diurnal tide dominating the lower thermosphere and an in-situ generated diurnal tide in the upper thermosphere. The results also show considerable interaction between dynamics and composition, especially atomic oxygen between 85 and 120 km.

Roble, R. G.

The effects of neutral inertia on ionospheric currents in the high-latitude thermosphere following a geomagnetic storm

Neutral flywheel effects are investigated in NCAR-TIGCM simulation of geomagnetic storms that occurred in November 23, 1982 and December 7-8, 1982. Theoretical calculations from the latter storm are compared with measurements of currents form instruments on the Dynamics Explorer 2 satellite. It is concluded that neutral flywheel effects can make a contribution to high latitude electrodynamics for a few hours after the main phase of a geomagnetic storm. The Hall currents that are driven by neutral winds during B(Z) northward conditions are generally in the opposite direction to those that occur during B(Z) southward conditions, when they are driven primarily by ion winds. The morphology of the field-aligned current system calculated by the NCAR-TIGCM during southward B(Z) conditions is in general agreement with observations.

Deng, W.

The Effects of Neutral Inertia on Ionospheric Currents in the High-Latitude Thermosphere Following a Geomagnetic Storm

Results of an experimental and theoretical investigation into the effects of the time dependent neutral wind flywheel on high-latitude ionospheric electrodynamics are presented. The results extend our previous work which used the National Center for Atmospheric Research Thermosphere/Ionosphere General Circulation Model (NCAR TIGCM) to theoretically simulate flywheel effects in the aftermath of a geomagnetic storm. The previous results indicated that the neutral circulation, set up by ion-neutral momentum coupling in the main phase of a geomagnetic storm, is maintained for several hours after the main phase has ended and may dominate height-integrated Hall currents and field-aligned currents for up to 4-5 hours. We extend the work of Deng et al. to include comparisons between the calculated time-dependent ionospheric Hall current system in the storm-time recovery period and that measured by instruments on board the Dynamics Explorer 2 (DE 2) satellite. Also, comparisons are made between calculated field-aligned currents and those derived from DE 2 magnetometer measurements. These calculations also allow us to calculate the power transfer rate (sometimes called the Poynting flux) between the magnetosphere and ionosphere. The following conclusions have been drawn: (1) Neutral winds can contribute significantly to the horizontal ionospheric current system in the period immediately following the main phase of a geomagnetic storm, especially over the magnetic polar cap and in regions of ion drift shear. (2) Neutral winds drive Hall currents that flow in the opposite direction to those driven by ion drifts. (3) The overall morphology of the calculated field-aligned current system agrees with previously published observations for the interplanetary magnetic field (IMF) B(sub Z) southward conditions, although the region I and region 2 currents are smeared by the TI(ICM model grid resolution. (4) Neutral winds can make significant contributions to the field-aligned current system when B(sub Z) northward conditions prevail following the main phase of a storm, but can account for only a fraction of the observed currents. (5) DE 2 measurements provide a demonstration of "local" (satellite-altitude) flywheel effects. (6) On the assumption that the magnetosphere acts as an insulator, we calculate neutral-wind-induced polarization electric fields of approx. 20-30 kV in the period immediately following the geomagnetic storm.

Deng, W.

On the coupling between the lower and the upper thermosphere during the first lower thermosphere coupling study

The study reports simultaneous observations of the neutral winds and temperatures and electron densities in the upper and lower thermosphere. The measurements were made at Arecibo and Millstone Hill during the first Lower Thermosphere Coupling Study (LTCS 1) campaign in September 1987. The observations show much day-to-day variability, possibly due to the geomagnetic activity which occurred during the period; storm effects are discernible in the data at low altitudes even at low latitudes. Simulations of the low and midlatitudes did not substantially improve with the use of lower boundary conditions derived from the measurements. The model predicts that diurnal waves are not negligible in analysis of lower thermosphere data; near 100 km the diurnal wave amplitudes of the meridional waves were predicted to be greater than the semidiurnal.

Fesen, C. G.

Thermospheric composition changes seen during a geomagnetic storm

The largest magnitude winds observed using the instruments on board the Dynamics Explorer 2 (DE-2) satellite were measured during the large geomagnetic storm that occurred on November 24, 1982. Neutral temperatures exceeded 2000 K during this storm, and these high temperatures, combined with the very large observed winds and the very full instrumental coverage available in both hemispheres, make it a unique event to study. This paper presents results obtained using these DE-2 data and a time dependent simulation of the event made using the NCAR Thermosphere/Ionosphere General Circulation Model (NCAR-TIGCM). In general, the agreement between model calculations and the data is very good, implying that most of the important physical processes controlling the energetics and dynamics of the thermosphere are reasonably well represented in the model. The modeled summer hemisphere changes in the mass mixing ratio of N2 (Psi(N2)) are in very good agreement with the averaged data in both hemispheres. This paper describes model-experiment comparisons for the November 24, 1982 geomagnetic storm changes.

Burns, A. G.

Thermospheric heating away from the auroral oval during geomagnetic storms

The data on the morphology of the changes in the heating of thermosphere in the middle and high latitudes during geomagnetic storms, collected by the DE-2 satellite between July 1981 and February 1983 are compared with the NCAR-Thermosphere-Ionosphere General Circulation Model simulation of the November 24, 1982 storm on a one-on-one basis for an individual orbit in the middle of this storm. Good agreement was found for the winter hemisphere. A thermodynamic diagnostic processor was then used to investigate the mechanism by which geomagnetic storms cause temperature increases at lower latitudes. It was found that the pattern of heating during geomagnetic storms is complex, indicative of the complex nature of the physical processes that alter the thermal structure of the thermosphere at these times. Model predictions indicated that the greatest temperature increase at 40 deg N occurs in the morning, in agreement with results of Proelss (1984).

Burns, A. G.

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.

A thermosphere/ionosphere general circulation model with coupled electrodynamics

A new simulation model of upper atmospheric dynamics is presented that includes self-consistent electrodynamic interactions between the thermosphere and ionosphere. This model calculates the dynamo effects of thermospheric winds, and uses the resultant electric fields and currents in calculating the neutral and plasma dynamics. A realistic geomagnetic field geometry is used. Sample simulations for solar maximum equinox conditions illustrate two previously predicted effects of the feedback. Near the magnetic equator, the afternoon uplift of the ionosphere by an eastward electric field reduces ion drag on the neutral wind, so that relatively strong eastward winds can occur in the evening. In addition, a vertical electric field is generated by the low-latitude wind, which produces east-west plasma drifts in the same direction as the wind, further reducing the ion drag and resulting in stronger zonal winds.

Richmond, A. D.

The effect of solar cycle variations on the thermosphere

A review of current knowledge of the effect of solar activity on the thermosphere is presented. Data obtained by the Dynamics Explorer 2 spacecraft and model results from the NCAR Thermospheric GCM, which demonstrate the effect of increased solar and geomagnetic activity on composition, temperature and the dynamic structure of the upper atmosphere are discussed. It is shown how subtle change in the IMF can affect the form of the neutral gas circulation at high latitudes, with important implications for polar orbiting spacecraft.

Mccormac, F. G.

Thermospheric dynamics, energetics, and composition at auroral latitudes

Recent experimental results from the Dynamics Explorer-2 (DE-2) spacecraft and theoretical calculations from the NCAR-TGCM (thermospheric general circulation models) are discussed. DE-2 observations show that the close coupling exists between the magnetosphere and the thermosphere, with magnetospheric energy and momentum sources playing a key role in establishing the thermospheric dynamical and thermal structure via ion-neural collisions. The tight ion-neutral momentum coupling results in the neutral-wind pattern which is strongly dependent on the level of geomagnetic activity and the orientation of the interplanetary magnetic field. Thermospheric temperatures are directly controlled by the Juole and particle heat sources associated with auroral latitudes. It is concluded that the various dynamical, energetic, and compositional variations associated with auroral processes are highly coupled with each other.

Killeen, T. L.

Thermospheric dynamics - A system theory approach

A system theory approach to thermospheric modeling is developed, based upon a linearization method which is capable of preserving nonlinear features of a dynamical system. The method is tested using a large, nonlinear, time-varying system, namely the thermospheric general circulation model (TGCM) of the National Center for Atmospheric Research. In the linearized version an equivalent system, defined for one of the desired TGCM output variables, is characterized by a set of response functions that is constructed from corresponding quasi-steady state and unit sample response functions. The linearized version of the system runs on a personal computer and produces an approximation of the desired TGCM output field height profile at a given geographic location.

Codrescu, M.