Nitric oxide and lower ionosphere quantities during solar particle events of October 1989 after rocket and ground-based measurements

The most dramatic demonstrations of solar activity are solar proton flares. One such very strong flare, accompanied by a solar proton event (SPE) and a large ground level enhancement of cosmic rays on Earth, was observed in October 1989. During this SPE, ion density and nitric oxide concentration profiles were measured by rockets launched from the Soviet research vessel ‘Akademik Shirshov’ in the southern part of the Indian Ocean. The rocket experiment yielded the first in-situ measurement of NO concentration increased by SPE. The NO concentrations estimated from ion-pair production rates due to measured fluxes of high energy particles agree fairly well with the observed NO concentrations in the stratopause region. The results of rocket measurements are compared with measurements of the radio wave absorption in the lower ionosphere performed at similar latitudes in central Europe. Model calculations of absorption show that while the night-time enhancement of absorption can be explained by increased electron density related to the measured increase of ion density as a consequence of enhanced penetration of high energy particles, the daytime increase of absorption needs to be explained mainly in terms of the observed increase of nitric oxide concentration.     

Solar flux and seasonal variations of the mesopause temperatures at 51°N

Observations of the OH (8-3) band airglow emission, using a multichannel tilting filter type photometer, have been carried out at Calgary (51°N, 114°W), Canada, since 1981. In this paper recent measurements of the nocturnal, seasonal and solar flux variations of the mesopause temperature, obtained from the rotational temperature of the OH (8-3) band observations, are presented. The data presented span the ascending phase of the present solar cycle viz. 1987–1988 (low solar activity) and 1990 (high solar activity). Good correlations (r = 0.73) between the OH (8-3) band rotational temperature and the 10.7 cm solar flux were observed. The mean temperature for the period investigated was about 210 K. The seasonally averaged nocturnal variations show only small irregular excursions, possibly associated with solar tides and the passage of gravity waves in the mesopause region. However, the observed rotational temperatures show considerable night-to-night changes.     

Soft X-rays from the sunlit earth's atmosphere

The HEAO-1 A-2 experiment low energy proportional counters have been used to measure the X-ray spectrum of the sunlit earth in the energy range 0.2–0.8 keV. The X-rays arise by coherent scattering of, or fluorescence of atmospheric constituents by, solar coronal X-rays incident on the atmosphere. Although the relative spectral contributions of the two processes depend upon the sun-earth-satellite geometry, fluorescent oxygen and nitrogen K X-ray emission is always important. The observed spectra were compared with calculations in order to derive the coronal temperature and emission measure, parameters that characterize the incident solar spectrum. These derived parameters agree well with the expected values for the nonflaring sun, and good agreement was obtained between measurements closely spaced in time but having a wide range of geometries and counting rates. Thus X-ray observations of the sunlit earth's atmosphere can be a useful monitor of solar activity for satellite-borne instrumentation unable to view the sun directly. The total measured fluorescent line flux agreed well with calculations, but the N : O line ratio did not. This disagreement is attributed to several causes, including the relative weakness of N emission at high altitudes where fluorescence is particularly important, the presence of line emission in the solar spectrum, and possible small calibration errors. Since present detectors cannot resolve the oxygen and nitrogen K X-rays, observation of X-ray fluorescence from the sunlit atmosphere will be useful in monitoring atmospheric constituents only to the extent that total line counting rates depend upon composition. X-rays from the sunlit earth are briefly examined as a source of background in auroral X-ray observations. During nonflare periods this background should be unimportant above about 2 keV.     

Day-to-day ionospheric variations in a period of high solar activity

In November and December 1979 the solar 10.7 cm radio flux density, sunspot number, X-ray flux and EUV flux were high and very variable. The day-to-day variations of noon F2-layer height and Elayer electron density at three ionosonde stations (Slough, Port Stanley and Huancayo) are found to be well correlated with the day-to-day variations of solar activity. The short-term E- and F-layer variations are consistent with those derived from longer-term studies.     

The 27 day solar UV response of stratospheric ozone: solar cycle 21 vs. solar cycle 22

A correlative study of ozone and the solar UV flux on the time scale of a solar rotation shows an anomalous response of ozone in the upper stratosphere during solar cycle 22. The study, which is based on the analysis of ozone and solar UV flux measured by the SBUV/2 spectrometer on NOAA 11 (January 1989–December 1990), shows a sharp transition from an in-phase relation between ozone and the solar UV flux below 2 mb to an almost out-of-phase relation above 1 mb. Such a phase change is not predicted by photochemical models and was not observed during solar cycle 21. The ozone measurements from the Nimbus-7 SBUV spectrometer from 1979 to 1984 showed an almost in-phase relation between ozone and the solar UV flux at these heights (in agreement with model predictions). Similar studies of ozone and temperature relations between 30 and 1 mb did not show significant changes from the solar cycle 21 to 22. The temperature oscillations appear to be primarily of dynamical origin, with no apparent correlation with solar UV flux.     

Splitting of atomic oxygen height profiles by eddy diffusion

A time-dependent model of upper atmospheric composition has been used to explain the splitting of the atomic oxygen concentration height profile into two maxima. The major feature of this model is multicomponent molecular diffusion. This model of the major constituents makes it possible to obtain the wave height profiles of the atomic oxygen as a result of the smooth eddy diffusion, multicomponent molecular diffusion and photochemistry. The position of maxima of the split profiles depends very strongly on the gradient of the eddy diffusion coefficient above its maximum. The depth of the splitting increases with the gradient increase, but the splitting does not depend on the maximum value of the eddy diffusion coefficient itself. It is important to note that, in this case, the atomic oxygen concentration may be varied by a factor of 2 in the thermosphere due only to the variation of the eddy diffusion coefficient gradient above the maximum. The splitting of the atomic oxygen height profiles is reflected in the height profiles of the Herzberg band and also of the green line emissions which are a result of atomic oxygen three-body recombination.     

The middle atmospheric response to short and long term solar UV variations: analysis of observations and 2D model results

We have investigated the middle atmospheric response to the 27-day and 11-yr solar UV flux variations at low to middle latitudes using a two-dimensional photochemical model. The model reproduced most features of the observed 27-day sensitivity and phase lag of the profile ozone response in the upper stratosphere and lower mesosphere, with a maximum sensitivity of +0.51% per 1% change in 205 nm flux. The model also reproduced the observed transition to a negative phase lag above 2 mb, reflecting the increasing importance with height of the solar modulated HO_x chemistry on the ozone response above 45 km. The rnodel revealed the general anti-correlation of ozone and solar UV at 65–75 km, and simulated strong UV responses of water vapor and HO_x species in the mesosphere. Consistent with previous 1D model studies, the observed upper mesospheric positive ozone response averaged over ±40° was simulated only when the model water vapor concentrations above 75 km were significantly reduced relative to current observations. Including the observed temperature-UV response in the model to account for temperature-chemistry feedback improved the model agreement with observations in the middle mesosphere, but did not improve the overall agreement above 75 km or in the stratosphere for all time periods considered. Consistent with the short photochemical time scales in the upper stratosphere, the model computed ozone-UV sensitivity was similar for the 27-day and 11-yr variations in this region. However, unlike the 27-day variation, the model simulation of the 11-yr solar cycle revealed a positive ozone-UV response throughout the mesosphere due to the large depletion of water vapor and reduced HO_x-UV sensitivity. A small negative ozone response at 65–75 km was obtained in the 11-yr simulation when temperature-chemistry feedback was included,In agreement with observations, the model computed a low to middle latitude total ozone phase lag of +3 days and a sensitivity of +0.077% per 1% change in 205 nm flux for the 27-day solar variation, and a total ozone sensitivity of +0.27% for the 11-yr solar cycle. This factor of 3 sensitivity difference is indicative of the photochemical time constant for ozone in the lower stratosphere which is comparable to the 27-day solar rotation period but is much shorter than the 11-yr solar cycle.     

EISCAT incoherent scatter radar observations and model studies of day to twilight variations in the D-region during the PCA event of August 1989

An intense solar proton event causing enhanced ionization in the ionospheric D-region occurred on 12 August 1989. The event was partially observed during three successive nights by the EISCAT UHF incoherent scatter radar at Ramfjordmoen near Tromsa, Norway. Ion production rates calculated from GOES-7 satellite measurements of proton flux and a detailed ion chemistry model of the D-region are used together with the radar data to deduce electron concentration, negative ion to electron concentration ratio, mean ion mass and neutral temperature in the height region from 70 to 90 km, at selected times which correspond to the maximum and minimum solar elevations occurring during the radar observations. The quantitative interpretation of EISCAT data as physical parameters is discussed. The obtained temperature values are compared with nearly simultaneous temperature measurements at Andøya based on lidar technique.     

A model of the lower ionosphere above Red Lake,Canada, during the 26 February 1979 solar eclipse

Measurements of ionization sources, ionization profiles and minor atmospheric constituents were conducted during the 26 February 1979 solar eclipse above Red Lake, Canada. A model of the lower thermosphere has developed to describe the D- and E-regions of the ionosphere for this case with the model being guided by the measurements. During the eclipse a rather intense particle precipitation event was in progress. For this reason, an auroral deposition code was coupled to a chemical-kinetics code to calculate degraded primary and secondary electron fluxes, ionization rates, positive ion and electron densities. The model was calibrated with the experimental measurements of electron flux below 100 km and electron density between 70 and 150 km. This calculation not only satisfactorily described the ionization in the E-region but also the gross electron density characteristics of the D-region. Bursts in the observed electron flux were also simulated with the model to give electron density profiles that were remarkably consistent with small perturbations seen in the electron density measurements.     

Short term variabilities of ionospheric electron content (IEC) and peak electron density (NP) during solar cycles 20 and 21 for a low latitude station

Hourly values of IEC and of f₀F₂ (critical frequency) for a low latitude station, Hawaii (21.2°N, 157.7°W), during the solar maxima (1969 and 1981) and minima (1965 and 1985) years of two consecutive solar cycles, 20 and 21, are used to study the day to day variabilities of the ionospheric parameters IEC and NP. It is found that there is good correspondence in the day to day variations of IEC and NP from one solar cycle to the other for both solar maximum and minimum years in the two solar cycles. Depending on solar phase and season, while the mean daytime IEC and NP variations range from about 20% to 35%, the mean night time values vary from about 25% to 60%. The mean daytime variations in NP for the solar minimum phase are remarkably higher in all the three seasons compared to the solar maximum phase. However, no such increase is observed in the mean daytime IEC variations, indicating the highly variable nature of the daytime ionospheric F region compared to the topside during solar minimum for this low latitude station. The winter night time IEC also seems to be a relatively stable parameter during the solar minimum. The short term day to day variabilities of the day time peak values of IEC and NP (ie IEC_max and NP_max) are not closely associated with the variations in F_10.7 solar flux. Contrary to the common expectation, the variabilities in both the parameters, particularly in NP_max, are somewhat reduced during the solar maximum (when the variability in F_10.7 solar flux is much higher compared to the solar minimum) which is more evident in the stronger 21 solar cycle. A larger number of significant components are seen in the spectra of the percentage variation of both IEC_max and NP_max during both solar phases of the two solar cycles compared to the corresponding F_10.7 solar flux spectra. The number of additional components for both the parameters with periods less than 15 days are more for the low solar activity years than for the solar maximum years.     

Can the high latitude ionosphere support large field-aligned ion drifts?

Ion velocities perpendicular and parallel to the geomagnetic field have recently been deduced by Smith et al. from bistatic measurements at 71° geomagnetic latitude in the afternoon sector. The results of this experiment include large (>400 m s⁻¹) downward ion velocities parallel to the magnetic field that persist for hours, small (100 m s⁻¹) ion velocities perpendicular to the magnetic field and electron density profiles with extremely narrow full-width at half-maximum. The explanation of these results was that the ionospheric flux tubes observed were near the terminator, and thus, sunlit at the top and in darkness at the bottom. The difference in production between the top and bottom of the flux tube creates an excess of ions at the top, which rapidly diffuse downwards. A three-dimensional, time-dependent model of the ionosphere has been used to test this explanation. Numerical experiments were performed to determine upper limits for the downward ion velocity. Assuming reasonable vertically-induced ion drifts due to either neutral winds or plasma convection, these upper limits were substantially smaller than the measurements. The location of the terminator was found to contribute a maximum of about 60 m s⁻¹ to the vertical ion velocity due to diffusion in a partially illuminated flux tube. In an attempt to explain the narrow density profiles without invoking an additional ionization source, the downward force in the model was arbitrarily increased, as would occur due to parallel electric fields in the ionosphere. Since the interpretation of these measurements as large field-aligned flows seems untenable by a model thought to be consistent with the currently accepted physics of the atmosphere, an alternate hypothesis is presented. 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. Therefore, the field-aligned velocities would be small, and the narrow width of the profiles would be due to enhanced electron densities in an O⁺ arc.     

Oscillations in D-region absorption at periods of one to two months

One to two month oscillations in D-region absorption are found in seven years of daily ƒ-min data from low latitude stations at Singapore (1°N, 104°E) and Rarotonga (21°S, 160°W). Coherency (cross-spectral) analyses reveal that solar flux variations account for much of the ƒ-min variance at these periods. Over the range of periods from 10 to 200 days, statistically significant linear correlation is found between the ƒ-min time series and contemporaneous 10.7cm solar flux measurements at periods of 16–19 days, the 26–29 day solar rotation band, and a broad band covering 43–80 day periods.     

Diurnal,semi-diurnal and terdiurnal Hough components of surface pressure

Diurnal and semi-diurnal Hough components of surface pressure are evaluated by classical tidal theory for previously presented profiles of water vapour and ozone heating. Values are compared with the observational results of Haurwitz and Cowley (1973) and show significant discrepancies which are considered to indicate the need for additional heating to be identified. The present calculations are in satisfactory agreement with those for semi-diurnal modes evaluated by Walterscheid et al. (1980).Profiles of terdiurnal heating due respectively to ozone and water vapour absorption of solar radiation are calculated for the (3,3,3) to (3,3,6) Hough modes and corresponding surface pressure oscillations are evaluated by classical tidal theory. Comparisons are made with earlier evaluations and with observationally derived results. The calculated solstitial (3, 3, 4) mode, which is the mode of the largest surface pressure amplitude, shows good agreement in phase with observation but underpredicts observed amplitudes by about 36% in contrast to earlier evaluations which were based on a now unacceptable basic temperature profile.     

Effects of solar tides on the zonal mean circulation in the lower thermosphere: solstice condition

Effects of momentum deposition due to solar diurnal and semi-diurnal tidal waves on the zonal mean circulation in the mesosphere and lower thermosphere for a solstice condition are discussed. In the present model, the system of zonally averaged equations and the system of perturbation equations are integrated simultaneously, so that the propagation of tidal waves is affected not only by the basic mean fields but also by the induced zonal mean fields due to the momentum deposition. Results for two different vertical eddy diffusion profiles are presented. It is shown that the solar tides make a significant contribution to the generation of the mean zonal winds in the upper mesosphere and the lower thermosphere. Below 120 km the main contribution is due to propagating diurnal tides, while above 120 km it is due to semidiurnal tides.     

On the global atmospheric electrical circuit

A model of global atmospheric electric circuit has been developed, which incorporates the pollution due to aerosol particles of anthropogenic and volcanic origins and the ionization caused by the coronal discharges due to intense electric field beneath thunderclouds in the atmospheric layers close to the earth's surface. Also the effects of solar activity and of Stratospheric Aerosol Particles (SAP) on the parameters of the circuit have been studied. The global distribution of Aerosol Particles (AP) of man-made origin has been assumed on the basis of world population density and that of volcanic origin on the basis of global distribution of volcanic activities. The thunderstorms are assumed to be the current generators of the global circuit. The latitudinal, longitudinal and height variations of atmospheric conductivity have been calculated from the known variations of ionization caused by cosmic rays, the radio-active emanations and intense electric fields beneath thunderclouds (over continents), along with the assumed variation of AP concentration.On increasing the AP concentration over the northern hemisphere by about 5 times that of the southern hemisphere, the ionospheric potential increases to about 6% and electric field increases in non-mountainous regions more than that in the high mountain regions. A large increase in SAP serves to increase the global resistance while both global current and ionospheric potential decrease. However, for low SAP concentration, the ionospheric potential increases. The SAP affect the electrical structure of the stratosphere without much influencing the troposphere except in the volcanically active regions, where conductivity is low due to high AP concentration. The major influence on the global atmospheric electric circuit occurs in response to solar activity. In the absence of local effects, the calculated variation of ionospheric potential is 14%. However, in real atmosphere (where both local and global processes act simultaneously), the ionospheric potential is found to vary by only about 3%. This has little effect on the ground electrical properties where more than 30% of the variations have been found to be caused by the local effect.     

Changes of pulsation activity during two solar cycles

Monthly and yearly averages of the pulsation activity in the mid-latitude station Nagycenk are compared to solar wind velocity and ionospheric-plasmaspheric electron concentration data. It is found that pulsation amplitudes are correlated with solar wind velocities with the exception of some month around December in solar maximum years, when they are significantly lower than computed from the corresponding solar wind velocities. This decrease can be caused either by a cutoff of the magnetospheric shell resonances or by local ionospheric damping. In addition to these effects, pulsations amplitudes slightly depend also on the geomagnetic activity and have a semi-annual activity change with maxima around equinoxes.     

On the IRI model predictions for the low-latitude ionosphere

Measurements of ionospheric electron density vertical profiles, carried out at a magnetic equatorial station located at Fortaleza (4°S, 38°W; dip latitude 2°S) in Brazil, are analyzed and compared with low-latitude electron density profiles predicted by the International Reference Ionosphere (IRI) model. The analysis performed here covers periods of high (1979/1980) and low (1986) solar activities, considering data obtained under magnetically quiet conditions representative of the summer, winter and equinox seasons. Some discrepancies are found to exist between the observed and the IRI model-predicted ionospheric electron density profiles. For high solar activity conditions the most remarkable one is the observed fast upward motion of the F-layer just after sunset, not considered in the IRI model and which precedes the occurrence of nighttime ionospheric plasma irregularities. These discrepancies are attributed mainly to dynamical effects associated with the low latitude E × B electromagnetic plasma drifts and the thermospheric neutral winds, which are not satisfactorily reproduced either in the CCIR numerical maps or in the IRI profile shapes. In particular, the pre-reversal enhancement in the vertical E × B plasma drifts around sunset hours has a great influence on the nighttime spatial distribution of the low-latitude ionospheric plasma. Also, the dynamical control exerted by the electromagnetic plasma drifts and by the thermospheric neutral winds on the low-latitude ionospheric plasma is strongly dependent on the magnetic declination angle at a given longitude. These important longitudinal and latitudinal dependences must be considered for improvement of IRI model predictions at low latitudes.     

Recipe for predicting the IMF Bz polarity's change of direction following solar disturbances and at the onset of geomagnetic storms

A three-dimensional, time-dependent, MHD model of solar-disturbance-caused storms (Wu, 1993; Wu et al., 1996a) is used to predict the turning direction of the interplanetary magnetic field (IMF) at Earth. More explicitly, we examine the polarity of B_z caused by solar disturbances on the Sun. Three manifestations of solar disturbances, as studied by previous workers, are examined. Firstly, twenty-nine kilometric Type II events, associated (Cane, 1985) with geomagnetic storms, are studied within the context of our three-dimensional model. Then, an additional eleven long-duration X-ray events (LDEs) with radio fluxes greater than 100 solar flux units were examined; these events were not associated with interplanetary Type II events but were also associated (Cane, 1985) with geomagnetic storms. Finally, in situ interplanetary phenomena that caused ten large (Dst < −100 nT, the intensification of the storm) geomagnetic storm episodes (Tsurutani et al., 1988) near solar maximum are also studied via the B_z predictions of our 3D MHD model. The accuracy of these B_z turning-direction-predictions is found to be as follows: (1) for the kilometric Type II events, the model's prediction was successful for 26 of the 29 events studied; (2) 10/11 for the LDE events; and (3) 7/9 for the major geomagnetic storm events. The overall prediction accuracy of these three independent data sets is 43/49. Thus, consideration of these three independent data sets strongly suggests that the recipe proposed by the basic 3D MHD model may be valid for a zero-th order prediction scheme.     

Vertical plasma drifts in the post-sunset F-region at the magnetic equator

HF doppler observations of vertical plasma drifts in the post-sunset equatorial F-region at Trivandrum (dip 0.9°S), conducted over a range of solar and geomagnetic conditions, are presented. The observations show that under magnetically quiet conditions, the characteristic post-sunset enhancement in the vertical plasma drift is quite sensitive to solar activity; the peak velocity drops by about a factor of 3 as the solar flux index (S_10.7) changes from about 125 to 70. It is found that the drift velocity enhancement has strong magnetic activity dependence only during high solar activity; the drift velocity drops by more than a factor of 2 from quiet to moderate activity, but builds back to the quiet day level for high magnetic activity. The occurrence of equatorial spread-F (ESF) is seen to be closely linked to the post-sunset enhancement in the vertical drift velocity, both showing essentially the same dependence on solar and magnetic activities. A comparison with Jicamarca observations shows that while the gross characteristics of the drift velocity pattern are about the same for the two stations, there are significant differences in the detailed variations, particularly for magnetically disturbed conditions.     

Neutral density cells in the high latitude thermosphere—1. Solar maximum cell morphology and data analysis

A new high latitude thermospheric neutral density structure has been revealed in NCAR-TIGCM simulations at 120–350 km altitude. The structure consists of density cells above 50° latitude with radii of approximately 1000 km. There are between two to four cells present depending on the altitude and magnetic activity. For example, at 200 km under magnetically active conditions, the density structure consists of four cells: low density cells are located near dawn and dusk and high density cells are located near noon and midnight. Density variations among the cells range from 5 to 50% for magnetically quiet and active conditions respectively. The cells are present at all seasons, for a wide range of magnetic activity levels, and at solar minimum and solar maximum. The density cell morphology is established for equinox solar maximum as a function of altitude and magnetic activity. Departures of the cell structure from this morphology due to seasonal and solar cycles are discussed. The cell morphology provides a new framework in which to interpret lower thermospheric density data. Data to test and confirm the model predictions were provided by the SETA-1 satellite.