Vertical motions in the thermosphere over Mawson,Antarctica

High time resolution measurements of Doppler shift and broadening of the (OI) >1630 nm emission in the night airglow and aurora have provided determinations of vertical velocities and temperatures in the neutral thermosphere over Mawson, Antarctica. The vertical wind exhibits a large, rapid and complex response to geomagnetic energy input. Upward winds greater than 50 m s⁻¹ are frequently associated with the expansion phase of auroral substorms. Following the disturbance, prolonged periods of downward winds produce temperature enhancements of 200K outside the source region, thus providing a mechanism for the redistribution of geomagnetic energy. Oscillatory behaviour consistent with thermospheric gravity waves is observed during both quiet and disturbed conditions.     

Observations of thermospheric neutral winds by the UCL Doppler imaging system at Kiruna in northern Scandinavia

Since the 1982/1983 winter, the UCL group, in collaboration with the Swedish Institute for Space Physics (previously Kiruna Geophysical Institute), has operated a Doppler imaging system at the high latitude station of Kiruna (67°N, 22°E). The Doppler imaging system is an imaging Fabry-Perot interferometer of 13.2 cm aperture. This instrument has been operated on a ‘campaign’ basis for mapping thermospheric winds using the OI emission at 630 nm (240 km altitude) from a region up to about 400 km radius about Kiruna. In November 1986, the performance of this wide-field Doppler imaging system was augmented by improvements to the detector and all-sky optics. We present data from December 1986, obtained during periods with both clear skies and active auroral and geomagnetic conditions. Maps of the neutral wind flow within the auororal oval during disturbed conditions and near magnetic midnight show continuous and rapid changes of thermospheric winds. The typical scale sizes of eddies observed within the mean flow around magnetic midnight are 100–300 km, with fluctuations at all time scales resolved by the 10 min between successive Doppler images. The local and short period fluctuations appear to be a filtered response of the thermosphere to rapid local variations of the convection and precipitation patterns, within a background of global scale changes     

Comparison of plasma flow and thermospheric circulation over northern Scandinavia using EISCAT and a Fabry-Perot interferometer

The EISCAT incoherent scatter radar, operating in a full tristatic mode, provided data on the ionospheric plasma drift above northern Scandinavia, during the 24 h period, 11 UT 25 November to 11 UT 26 November 1982. For the hours of darkness, 14 UT until 05 UT, observations of thermospheric winds were made by means of a ground-based Fabry-Perot interferometer (FPI) operated at Kiruna Geophysical Institute (21° E, 68° N). During this period, the radar observations describe well the ebbing and flowing of regions of strong convective ion flow associated with the auroral oval. As individual geomagnetic disturbances occur, the overall ion flow pattern intensifies and moves equatorward. The zonal thermospheric wind observed by the FPI responds rapidly to surges of the local ionospheric convection, while the meridional wind response is slower and apparently to much larger-scale features of the geomagnetic input to the high latitude thermosphere. From the data base, periods of strong heating of the ionospheric ions and of the thermospheric gas can be identified, which can be compared with Joule and particle heating rates deduced from the observations of ionospheric drifts, neutral winds, electron densities and auroral emission rates. A three-dimensional, time-dependent global thermospheric model is used to distinguish local and global features of the thermospheric wind field. Meridional and zonal wind components at 312 km may be theoretically derived from the EISCAT data using an appropriate model (MSIS) for neutral temperature. The EISCAT-derived meridional wind is within about 50 m s⁻¹ of the FPI observations throughout the period of joint observations. The EISCAT-derived zonal wind is systematically larger (by about 50%) than the FPI measurement, but the two independent measurements follow closely the same fluctuations in response to geophysical events until 03 UT, when the EISCAT solution is driven away from the FPI measurement by a sharp increase in both neutral and ion temperatures. Between 03 and 05 UT the EISCAT-derived zonal wind is 200–400 m s⁻¹ westward. Allowance for the neutral temperature rise would reduce the EISCAT values towards the very small zonal winds shown by the FPI during this period. We describe the relatively straightforward analysis required to derive the meridional wind from the radar data and the limitations inherent in the derivation of zonal wind, using the ion energy equation, due to the lack of precise knowledge of the background neutral temperature from the EISCAT data alone. For analysis of EISCAT ion drift observations at 312 km, the ground-based FPI temperature measurements do not improve the accuracy of the analysis, since the median altitude of the FPI measurement is probably in the range 180–240 km throughout the observation period. This median altitude and the temperature gradient both fluctuate in response to local geomagnetic events, while the temperature gradient may be considerably greater than that predicted by standard atmospheric models. When the neutral temperature is well known, or when there is a large enhancement of the ion temperature, the EISCAT-derived zonal wind exceeds the FPI measurement, but the consistency with which they correlate and follow ion-drag accelerations suggests that the differences are purely due to the considerable altitude gradients which are predicted by theoretical models.     

First summer results on winds in the upper mesosphere derived from the 843 nm hydroxyl emissions measured from the Bear Lake Observatory,Utah

The Imaging Fabry-Perot Interferometer (IFPI) at the Bear Lake Observatory (BLO), Utah (41.9°N, 111.4°W) is used for studies of the aeronomy of the middle and upper atmosphere. Wind and temperature structure can be determined from observations of the Doppler shift and Doppler broadening of the airglow and auroral emissions from the mesosphere and thermosphere. The mesospheric winds recorded at the end of August, September and early October 1992 are consistent with a semi-diurnal tidal variation. The amplitude of this variation is approximately 30 ms⁻¹ at the end of August and early September and approximately 20 ms⁻¹ at the end of September and early October. However, during June and July, the semi-diurnal tidal variation, if present, is weak, with amplitude < 5 ms⁻¹. No consistent semi-diurnal tidal variation is observed during late October 1992. During the solstice period, antisymmetric tidal components may be preferentially generated in such a way that they can result in destructive interference with the normally dominant symmetric modes, resulting in a decrease of tidal variation. This is consistent with the observed decrease in tides during the June, July and late October periods. Near the equinoxes, however, the excitation of these antisymmetric modes is expected to be weaker, possibly explaining why a pronounced and consistent semi-diurnal tidal variation has been observed during the August, September and early October periods. In contrast, the mesospheric winds derived from the Sheffield Meteor Wind Radar (53.4°N, 1.5°W) reveal a clear semi-diurnal tidal variation throughout the year, with an amplitude that may vary between 15 ms⁻¹ and 50 ms⁻¹, being about 25 ms⁻¹ on average. The IFPI records winds from a region of the atmosphere centred at 87 km, whereas the Sheffield Meteor Wind Radar measures winds centred at 95 km. Therefore, the two regions may experience different tidal modes due to the different latitude, longitude and altitude of the observed regions and/or the different topography of the observing sites. Some proposed reasons for these differences are presented.     

Simulations of the behaviour of the polar thermosphere and ionosphere for northward IMF

The dynamics and structure of the polar thermosphere and ionosphere within the polar regions are strongly influenced by the magnetospheric electric field. The convection of ionospheric plasma imposed by this electric field generates a large-scale thermospheric circulation which tends to follow the pattern of the ionospheric circulation itself. The magnetospheric electric field pattern is strongly influenced by the magnitude and direction of the interplanetary magnetic field (IMF), and by the dynamic pressure of the solar wind. Previous numerical simulations of the thermospheric response to magnetospheric activity have used available models of auroral precipitation and magnetospheric electric fields appropriate for a southward-directed IMF. In this study, the UCL/Sheffield coupled thermosphere/ionosphere model has been used, including convection electric field models for a northward IMF configuration. During periods of persistent strong northward IMF B_z, regions of sunward thermospheric winds (up to 200 m s⁻¹) may occur deep within the polar cap, reversing the generally anti-sunward polar cap winds driven by low-latitude solar EUV heating and enhanced by geomagnetic forcing under all conditions of southward IMF B_z. The development of sunward polar cap winds requires persistent northward IMF and enhanced solar wind dynamic pressure for at least 2–4 h, and the magnitude of the northward IMF component should exceed approximately 5 nT. Sunward winds will occur preferentially on the dawn (dusk) side of the polar cap for IMF B_y negative (positive) in the northern hemisphere (reverse in the southern hemisphere). The magnitude of sunward polar cap winds will be significantly modulated by UT and season, reflecting E-and F-region plasma densities. For example, in northern mid-winter, sunward polar cap winds will tend to be a factor of two stronger around 1800 UT, when the geomagnetic polar cusp is sunlit, then at 0600 UT, when the entire polar cap is in darkness.     

Fabry-Perot spectrometer observations of the auroral oval/polar cap boundary above Mawson,Antarctica

Zenith observations of the oxygen λ1630 nm auroral/airglow emission (produced at an altitude of ∼220 to ∼250 km) were obtained with the Mawson Fabry-Perot Spectrometer (FPS) during three ‘zenith direction only’ observing campaigns in 1993. The data show many instances of strong (50 to 100 m s⁻¹) upwellings in the vertical wind, when the auroral oval is located equatorward of the zenith. Our data appear consistent with the existence of a region of upwelling up to ∼ 4° poleward of the poleward boundary of the visible auroral oval, rather than short duration, explosive heating events. The upwellings are probably the vertical component of wind shear produced by reversal of the zonal thermospheric winds, which occurs near the poleward boundary of the visible auroral oval. Zenith temperature was also seen to increase when the oval was equatorward of Mawson, showing rises of up to 300 K or more. However, this increase is at times unrelated to the upwellings, and seems to be caused by the expansion of the warm polar cap over the observing site.On a number of nights the boundary between the polar cap and the auroral oval was observed to pass over our site several times, occasionally showing a quasi-periodic expansion and contraction. We speculate that this quasi-periodic movement may be related to periodic auroral activity that is known to generate large-scale gravity waves.     

Mesospheric wind studies during AIDA Act '89: morphology and comparison of various techniques

The Arecibo Initiative in Dynamics of the Atmosphere (AIDA) '89 was a multi-instrument campaign designed to compare various mesospheric wind measurement techniques. Our emphasis here is the comparison of the incoherent scatter radar (ISR) measurements with those of a 3.175 MHz radar operating a s an imaging Doppler interferometer (1131). We have performed further analyses in order to justify the interpretation of the long term IDI measurements in terms of prevailing winds and tides. Initial comparison of 14 profiles by Hines et al., 1993, J. atmos. terr. Phys. 55, 241–288, showed good agreement between the ISR and IDI measurements up to about 80 km, with fair to poor agreement above that altitude. We have compiled statistics from 208 profiles which show that the prevailing wind and diurnal and semidiurnal tides deduced from the IDI data provide a background wind about which both the IDI and ISR winds are normally distributed over the height range from 70 to 97 km. The 3.175 MHz radar data have also been processed using an interferometry (INT) technique [Van Baelen and Richmond 1991, Radio Sts. 26, 1209–1218] and two spaced antenna (SA) techniques [Meek, 1980, J. atmos. terr. Phys. 42, 837–839; Briggs. 1984, MAP Handbook, Vol. 13, pp. 166–186] to determine the three dimensional wind vector. These are then compared with the IDI results. Tidal amplitudes and phases were calculated using the generalized analysis of Groves, 1959, S. atmos. terr. Phys. 16, 344–356, historically used on meteor wind radar data. Results show a predominance of the diurnal S₁¹ tidal mode in the altitude range 70–110 km, reaching a maximum amplitude 45 ms⁻¹ at 95 km, with semidiurnal amplitudes being about 10–15 ms⁻¹ throughout the height range considered. There is evidence of the two day wave in data from 86–120 km, with amplitudes on the order of 20 ms⁻¹.     

Recent work on mid-latitude and equatorial sporadic-E

Theoretical and experimental work since 1970 is summarized. Mid-latitude sporadic-E is most likely due to a vertical shear in the horizontal east-west wind and this theory accounts for the detailed observations of the wind and electron density profiles. Preferred heights of sporadic-E are separated by about 6km and descending layers are often seen moving down with velocities in the range 0.6–4 ms⁻. Sometimes sporadic-E layers are very flat and uniform, and at other times form clouds of electrons 2–100km in size moving horizontally at 20–130 ms⁻¹. Sporadic-E is probably not correlated with meteor showers; this is a rather surprising result since the ions are meteor debris.The major problems with windshear theory are to account for the dramatic seasonal variation and, to a lesser extent, for the geographical and diurnal distributions.The Q-type equatorial sporadic-E appears to be due to the gradient instability. There is a very much smaller amount of new experimental data available in this area.     

The effect of the mid-latitude ionospheric trough on whistler mode ducting during magnetic storms

Whistler-mode signals observed at Faraday, Antarctica (65° S, 64° W, Λ=50.8°) show anomalous changes in group delay and Doppler shift with time during the main phase of intense geomagnetic activity. These changes are interpreted as the effect of refracting signals into and out of ducts near L=2.5 by electron concentration gradients associated with edges of the mid-latitude ionospheric trough. The refraction region is observed to propagate equatorwards at velocities in the range 20–85 ms⁻¹ during periods of high geomagnetic activity (K_p ≥ 5), which is in good agreement with typical trough velocities. Model estimates of the time that the trough edges come into view from Faraday show a good correlation with the observed start times of the anomalous features. Whistler-mode signals observed at Dunedin, New Zealand (46° S, 171° E, Λ=52.5°) that have propagated at an average L-shell of 2.2 (Λ=47.6°) do not show such trough-related changes in group delay. These observations are consistent with a lower occurrence of the trough at lower invariant latitudes.     

Thermosphere dynamics

Some recent investigations of thermosphere dynamics, carried out in the U.S.S.R., are reviewed briefly. The global empirical models of thermospheric motions are obtained on the basis of ground-based HF and meteor radar measurements of ionospheric irregularities drifts. Numerical modelling of large scale thermospheric electrodynamics for the low and mid-latitudes for quiet geomagnetic conditions, is presented. Disturbances of thermospheric wind systems from high latitude heat sources are considered. The response of lower thermosphere dynamics due to variations of solar and geomagnetic activity are discussed.     

First results with the Adelaide VHF radar: spaced antenna studies of tropospheric winds

The first results from a VHF radar of the ST type located at Buckland Park near Adelaide, Australia (35°S, 138°E), are presented. The radar is designed to be versatile and can be used to measure velocities in the lower atmosphere using both the spaced antenna (SA) and Doppler beam-swinging (DBS) techniques. Here studies of irregularities and motions made with the spaced antenna technique are discussed. It is shown that the scale of the diffraction pattern formed by the backscattered radiation varies with altitude, with the mean pattern scale being smaller in the troposphere than in the stratosphere. The observations are consistent with the backscattered energy decreasing as a function of off-vertical angle by 1.5 dB per degree in the troposphere and by about 2.8 dB per degree in the lower stratosphere. An intercomparison of zonal velocities measured with the SA and DBS methods shows good agreement. In May and August 1984 an extensive comparison was made between the velocities measured by the SA method and winds determined from over 80 balloon-borne radiosondes released from Adelaide Airport, situated some 36 km to the south of the radar. The velocities were compared on a statistical basis and showed excellent agreement, although the SA speeds tended to be 1–2 m s⁻¹ smaller in magnitude than the radiosonde velocities. Overall, the rms differences between the two sets of measurements was only 3–4ms⁻¹ throughout the troposphere, a result which is consistent with the random errors inherent in each technique, as well as the spatial separation between the radar and balloon observations. The utility of the SA method for meteorological observations is illustrated by a study of both the horizontal and vertical wind fields during the passage of a cold front made in November 1984. The high time resolution available with the radar allows detailed studies of the development of the pre-frontal jet, the wind convergence into the front and associated vertical motions.     

Optical interferometer measurements of thermospheric temperature at Kavalur (12.5°N, 78.5°E), India

First results on the behaviour of thermospheric temperature over Kavalur (12.5°N, 78.5°E geographic; 2.8°N geomagnetic latitude) located close to the geomagnetic equator in the Indian zone are presented. The results are based on measurements of the Doppler width of O(¹D) night airglow emission at 630 nm made with a pressure-scanned Fabry-Perot interferometer (FPI) on 16 nights during March April 1992. The average nighttime (2130-0430 IST) thermospheric temperature is found to be consistently higher than the MSIS-86 predictions on all but one of the nights. The mean difference between the observed nightly temperatures and model values is 269 K with a standard error of 91 K. On one of the nights (9/10 April 1992, Ap = 6) the temperature is found to increase by ~250 K around 2330 IST and is accompanied by a ‘midnight collapse’ of the F-region over Ahmedabad (23°N, 72°E, dip 26.3°N). This relationship between the temperature increase at Kavalur and F-region height decrease at Ahmedabad is also seen in the average behaviour of the two parameters. The temperature enhancement at Kavalur is interpreted as the signature of the equatorial midnight temperature maximum (MTM) and the descent of the F-region over Ahmedabad as the effect of the poleward neutral winds associated with the MTM.     

Role of vertical winds on the Rayleigh-Taylor mode instabilities of the night-time equatorial ionosphere

In view of the recent observations on the presence of vertical winds in the equatorial ionosphere in the evening and night-time, the role of vertical winds in the Rayleigh-Taylor (R-T) mode instability has been re-examined. The mathematical treatment of Chiu and Straus, earlier developd for a case of horizontal winds, is extended to evaluate the role of vertical winds in causing the R-T mode instability. It is shown that the vertical (downward) winds of small magnitude have a very significant effect on the instability growth rate in the. F-region. A downward wind of l m s⁻¹ can cause the same growth rate as a 200 m s⁻¹ eastward wind at 260 km altitude. Furthermore, a downward wind of 16m s⁻¹ at 300 km can be as effective as that due to the gravitational drift itself. Similarly, an upward wind can inhibit the instability on the bottomside of the F-region. It appears that the polarity of the vertical winds (upward or downward) at the base of the F-layer plays an important role in the growth of the R-T mode plasma instability in the equatorial ionosphere.     

Equatorial penetration of magnetic disturbance effects in the thermosphere and ionosphere

The neutral dynamic and electrodynamic coupling between high and low latitudes, and the mutual interactions between these two processes, are investigated. For 22 March 1979, when a sudden increase in magnetic activity occurred, we have analyzed the following experimental data: (a) neutral densities and cross-track neutral winds as a function of latitude (0°–80°) near 200 km from a satellite-borne accelerometer; (b) hourly mean H-component magnetic data from the Huancayo Observatory (0.72°S, 4.78°E; dipole geomagnetic coordinates) magnetometer; and (c) hourly mean foF2 measurements from the ionosonde at Huancayo. Comparisons are also made with a self-consistent thermosphere-ionosphere general circulation model and with observationally-based empirical models of winds and density.In concert with the increase in magnetic activity to K_p levels of 5–7, a nighttime (2230 LT) westward intensification of the neutral wind approaching 400 ± 100 ms⁻¹ occurred near the magnetic equator on 22 March 1979, accompanied by a 35% increase in neutral mass density. About 2 h after each of two substorm commencements associated with periods of southward IMF, ∼100γ and ∼200γ reductions in the daytime Huancayo H-component (corrected for ring current effects) are interpreted in terms of ∼0.5 and ∼1.0 mVm⁻¹ westward perturbation electric fields, respectively. An intervening 2-hour period of northward IMF preceded a positive equatorial magnetic perturbation of about 200γ. Time scales for field variations are a few hours, suggesting that processes other than Alfven shielding are involved. Variations in f₀F2 (∼ ± 1.0 MHz) over Huancayo are consistent with the inferred electric fields and magnetic variations. Similar equatorial perturbations are found through examination of other magnetic disturbances during 1979.     

Vertical neutral wind in the equatorial F-region deduced from electric field and ion density measurements

DC electric field and ion density measurements near density depletion regions (that is, equatorial plasma bubbles) are used to estimate the vertical neutral wind speed. The measured zonal electric field in a series of density depletions crossed by the San Marco D satellite at 01.47-01.52 UT on 25 October 1988, can be explained if a downward neutral wind of 15–30 m s⁻¹ exists. Simultaneously, the F-region plasma was moving downward at a speed of 30–50 m s⁻¹ These events appear in the local time sector of 23.002&#x0304;23.15 in which strong downward neutral winds may occur. Indeed, airglow measurements suggest that downward neutral velocities of 25–50 m s⁻¹ are possible at times near midnight in the equatorial F-region.     

Large simultaneous disturbances (LSDs) in the Antarctic ionosphere

A high frequency radio Doppler experiment was deployed in the Antarctic Peninsula region, centred on Argentine Islands (65°15′S, 64°16′W; L = 2.3), to investigate the morphology and sources of ionospheric disturbances. The experiment consisted of a three-transmitter dual frequency network which permits horizontal and vertical propagation velocities to be estimated over a north-south baseline of 200 km and an east-west baseline of 100 km.A new class of ionospheric disturbance has been observed, in the period range 10 min−1 h. These disturbances are characterised by unusually good correlation between perturbations on all available Doppler signals, but are apparently non—propagating and occur simultaneously at each reflection point. Several of these events display large (2 Hz at about 5 MHz transmitted frequency) Doppler shifts, thus we have labelled them Large Simultaneous Disturbances (LSD_s).Criteria for identification of LSD_s are established and the analysis of one event is described in detail. The occurrence statistics of the LSD_s are presented, including their seasonal and diurnal distributions.There is no clear general relationship between LSD_s and local geomagnetic field perturbations. However, examination of the magnetic indices AE and I_RC indicates that there is a loose association between the occurrence and amplitude of LSD_s and magnetic activity.Several possible mechanisms for the generation of LSD_s at middle latitudes are reviewed. The most likely explanation is that high latitude electric fields penetrate to magnetic middle latitudes and drive the ionospheric plasma via the E × B drift.     

Spectral characteristics of high frequency waves backscattered by small scale F-region irregularities: evidence of strong sub-auroral ion flow

The spectra of high frequency waves backscattered at night by small scale (10–20 m) sub-auroral F-region irregularities often exhibit large Doppler shifts and widths in the local time sector 2000–2400. After local midnight the Doppler shifts and the widths of the spectra decrease rapidly. We present examples of experimental data, obtained with the two coherent backscatter radars of the EDIA¹ experiment, showing the spectral characteristics just mentioned. From the Doppler shift measured at the two sites we deduced the perpendicular velocity of the irregularities, which can reach values as high as 2000 ms ⁻¹. These observations are interpreted using results of theoretical models which predict strong sub-auroral ion flow in the trough region.     

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.     

Summertime stratospheric wind measurements above the South Pole

Observations of the mean wind flow and wave motions in the stratosphere at the South Pole are presented. The atmospheric motions are determined from the tracking of a high altitude, zero-pressure balloon launched from Amundsen-Scott Station during the austral summer of 1985–1986. The balloon position was precisely monitored by an optical theodolite for a large portion of the flight so that small scale motions could be resolved. The mean flow above the pole was approximately 3ms⁻¹. Atmospheric motions characteristic of internal gravity waves were observed with an intrinsic period of approximately 4.5 h and vertical and horizontal wavelengths of approximately 2.5km and 125km, respectively. The horizontal perturbation velocity of the observed waves was large compared to the mean horizontal flow velocity. The implication is that wave motions play a dominant role in the transport of stratospheric constituents in regions where the mean winds are light, such as over the South Pole during austral summer.     

Estimation of initial auroral proton energy fluxes from Doppler profiles

An energetic auroral proton entering the atmosphere will, by charge exchange in collisions with atmospheric constituents, alternate between being a proton H⁺ and a neutral hydrogen atom H. This study provides a procedure to evaluate the auroral Doppler shifted and broadened hydrogen Balmer profile as a function of initial energy, flux, pitch angle and view angle relative to the geomagnetic field. The differential proton energy flux entering the atmosphere is deduced using ground-based measurements of H_α and H_β from Nordlysstasjonen in Adventdalen, Longyearbyen. The main assumptions are that the geomagnetic field lines are: parallel and vertical, and that the pitch angle of the H/H⁺-particle is preserved in collisions with atmospheric constituents before being thermalized. This numerical method estimates the fate of the auroral H/H⁺-particle in the atmosphere, and from measured Doppler profiles the corresponding incoming particle flux can be deduced. Optimization of the method will continue through extensive use of observational data.