Annual and seasonal variations in the low-latitude topside ionosphere

Annual and seasonal variations in the low-latitude topside ionosphere are investigated using observations made by the Hinotori satellite and the Sheffield University Plasmasphere Ionosphere Model (SUPIM). The observed electron densities at 600 km altitude show a strong annual anomaly at all longitudes. The average electron densities of conjugate latitudes within the latitude range ±25° are higher at the December solstice than at the June solstice by about 100% during daytime and 30% during night-time. Model calculations show that the annual variations in the neutral gas densities play important roles. The model values obtained from calculations with inputs for the neutral densities obtained from MSIS86 reproduce the general behaviour of the observed annual anomaly. However, the differences in the modelled electron densities at the two solstices are only about 30% of that seen in the observed values. The model calculations suggest that while the differences between the solstice values of neutral wind, resulting from the coupling of the neutral gas and plasma, may also make a significant contribution to the daytime annual anomaly, the E × B drift velocity may slightly weaken the annual anomaly during daytime and strengthen the anomaly during the post-sunset period. It is suggested that energy sources, other than those arising from the 6% difference in the solar EUV fluxes at the two solstices due to the change in the Sun-Earth distance, may contribute to the annual anomaly. Observations show strong seasonal variations at the solstices, with the electron density at 600 km altitude being higher in the summer hemisphere than in the winter hemisphere, contrary to the behaviour in NmF2. Model calculations confirm that the seasonal behaviour results from effects caused by transequatorial component of the neutral wind in the direction summer hemisphere to winter hemisphere.     

Ultraviolet-B solar irradiance,the influence of Ångström turbidity,at Central Spain

UV-B solar irradiance and meteorological variables were measured at the C.I.B.A. site (Low Atmosphere Research Laboratory), University of Valladolid, Spain, between January 2003 and March 2006. Calculated Ångström turbidity and aerosol optical thickness values were evaluated from the direct downward irradiance, surface pressure, air temperature and relative humidity values. Monthly turbidity β values for an average year showed minimum values in winter and maximum values in summer. The obtained values are according to an aerosol standard atmosphere between mean and clean continental model. UV-B model calculations were performed using a radiative transfer tropospheric model, TUV 4.1a; measured and calculated UV-B irradiance values were compared in order to establish possible values of single scattering albedo. The results show that calculated single scattering albedo values are according to an aerosol standard atmosphere between mean and polluted continental model.     

Annual and semiannual variations in the ionospheric F2-layer: II. Physical discussion

The companion paper by Zou et al. shows that the annual and semiannual variations in the peak F2-layer electron density (NmF2) at midlatitudes can be reproduced by a coupled thermosphere-ionosphere computational model (CTIP), without recourse to external influences such as the solar wind, or waves and tides originating in the lower atmosphere. The present work discusses the physics in greater detail. It shows that noon NmF2 is closely related to the ambient atomic/molecular concentration ratio, and suggests that the variations of NmF2 with geographic and magnetic longitude are largely due to the geometry of the auroral ovals. It also concludes that electric fields play no important part in the dynamics of the midlatitude thermosphere. Our modelling leads to the following picture of the global three-dimensional thermospheric circulation which, as envisaged by Duncan, is the key to explaining the F2-layer variations. At solstice, the almost continuous solar input at high summer latitudes drives a prevailing summer-to-winter wind, with upwelling at low latitudes and throughout most of the summer hemisphere, and a zone of downwelling in the winter hemisphere, just equatorward of the auroral oval. These motions affect thermospheric composition more than do the alternating day/night (up-and-down) motions at equinox. As a result, the thermosphere as a whole is more molecular at solstice than at equinox. Taken in conjunction with the well-known relation of F2-layer electron density to the atomic/molecular ratio in the neutral air, this explains the F2-layer semiannual effect in NmF2 that prevails at low and middle latitudes. At higher midlatitudes, the seasonal behaviour depends on the geographic latitude of the winter downwelling zone, though the effect of the composition changes is modified by the large solar zenith angle at midwinter. The zenith angle effect is especially important in longitudes far from the magnetic poles. Here, the downwelling occurs at high geographic latitudes, where the zenith angle effect becomes overwhelming and causes a midwinter depression of electron density, despite the enhanced atomic/molecular ratio. This leads to a semiannual variation of NmF2. A different situation exists in winter at longitudes near the magnetic poles, where the downwelling occurs at relatively low geographic latitudes so that solar radiation is strong enough to produce large values of NmF2. This circulation-driven mechanism provides a reasonably complete explanation of the observed pattern of F2 layer annual and semiannual quiet-day variations.     

Spectral energy contributions of quasi-periodic oscillations (2–35 days) to the variability of the f oF2

The relative contributions of quasi-periodic oscillations from 2 to 35 days to the variability of f_oF2 at middle northern latitudes between 42°N and 60°N are investigated. The f_oF2 hourly data for the whole solar cycle 21 (1976–1986) for four European ionospheric stations Rome (41.9°N, 12.5°E), Poitiers (46.5°N, 0.3°E), Kaliningrad (54.7°N, 20.6°E) and Uppsala (59.8°N, 17.6°E) are used for analysis. The relative contributions of different periodic bands due to planetary wave activity and solar flux variations are evaluated by integrated percent contributions of spectral energy for these bands. The observations suggest that a clearly expressed seasonal variation of percent contributions exists with maximum at summer solstice and minimum at winter solstice for all periodic bands. The contributions for summer increase when the latitude increases. The contributions are modulated by the solar cycle and simultaneously influenced by the long-term geomagnetic activity variations. The greater percentage of spectral energy between 2 to 35 days is contributed by the periodic bands related to the middle atmosphere planetary wave activity.     

The effects of nitric oxide cooling and the photodissociation of molecular oxygen on the thermosphere/ionosphere system over the Argentine Islands

In the past the global, fully coupled, time-dependent mathematical model of the Earths thermo-sphere/ionosphere/plasmasphere (CTIP) has been unable to reproduce accurately observed values of the maximum plasma frequency, foF2, at extreme geophysical locations such as the Argentine Islands during the summer solstice where the ionosphere remains in sunlight throughout the day. This is probably because the seasonal dependence of thermospheric cooling by 5.3 m nitric oxide has been neglected and the photodissociation of O₂ and heating rate calculations have been over-simplified. Now we have included an up-to-date calculation of the solar EUV and UV thermospheric heating rate, coupled with a new calculation of a diurnally varying O₂ photodissociation rate, in the model. Seasonally dependent 5.3 m nitric oxide cooling is also included. With these important improvements, it is found that model values of foF2 are in substantially better agreement with observation. The height of the F2-peak is reduced throughout the day, but remains within acceptable limits of values derived from observation, except at around 0600 h LT. We also carry out two studies of the sensitivity of the upper atmosphere to changes in the magnitude of nitric oxide cooling and photodissociation rates. We find that hmF2 increases with increased heating, whilst foF2 falls. The converse is true for an increase in the cooling rate. Similarly increasing the photodissociation rate increases both hmF2 and foF2. These changes are explained in terms of changes in the neutral temperature, composition and neutral wind.     

FFIM用于机载大气激光雷达 ABCT反演的可行性研究

Fernald前向积分法能否用于机载大气探测激光雷达气溶胶后向散射系数的反演一直是一个有争议的课题.本文利用青岛机载大气探测激光雷达实测数据、国外机载大气探测激光雷达实测数据和机载大气探测激光雷达模拟数据,对Fernald前向积分法应用于不同高度的机载大气探测激光雷达气溶胶后向散射系数反演的误差进行了定量分析,分析结果表明:飞机的飞行高度在3.5 km左右,标定值存在20%的误差时,离地面2 km的高度范围内反演得到的气溶胶后向散射系数的相对误差在12%以内,但在标定点附近相对误差可达20%;飞机飞行高度在7 km左右,当标定值存在100%的误差时,反演得到的气溶胶后向散射系数的相对误差大都在10%~15%之间,标定值存在400%的误差时,反演得到的气溶胶后向散射系数的相对误差大部分在15%~50%之间.本文从理论上对Fernald前向积分法应用于机载大气探测激光雷达气溶胶后向散射系数反演出现负值的原因进行了探讨.研究表明:Fernald前向积分法能够较准确地反演出中高空探测(4.5 km以上)机载大气探测激光雷达气溶胶后向散射系数,但应用于低空探测(4.5 km以下)机载大气探测激光雷达气溶胶后向散射系数反演时,反演误差较大甚至反演结果会出现负值.     

Seasonal variation of the auroral E-region neutral wind for different solar activities

Seasonal variations in the auroral E-region neutral wind for different solar activity periods are studied. This work is based on neutral wind data obtained over 56 days between 95–119 km altitude under geomagnetic quiet conditions (A_p<16) during one solar cycle by the European Incoherent Scatter radar located in northern Scandinavia. In general, the meridional mean wind shifts northward, and the zonal mean wind increases in eastward amplitude from winter to summer. The zonal mean wind blows eastward in the middle and lower E-region for each season and for each solar condition except for the equinox, where the zonal mean wind blows westward at and below 104 km. Solar activity dependence of the mean wind exists during the winter and equinox seasons, while in summer it is less prominent. Under high solar activity conditions, the altitude profiles of the horizontal mean winds in winter and the equinoxes tend to resemble those in summer. The horizontal diurnal tide is less sensitive to solar activity except during summer when the meridional amplitude increases by ∼10 m s⁻¹ and the corresponding phase shifts to a later time period (1–2 h) during high solar activity. Seasonal dependence of the semidiurnal tide is complex, but is found to vary with solar activity. Under low solar activity conditions the horizontal semidiurnal amplitude shows seasonal dependence except at upper E-region heights, while under high solar activity conditions it becomes less sensitive to seasonal effects (except for the meridional component above 107 km). Comparisons of mean winds with LF and UARS observations are made, and the driving forces for the horizontal mean winds are discussed for various conditions.     

The role of water-vapour photodissociation on the formation of a deep minimum in mesopause ozone

A one-dimensional atmospheric photochemical model with an altitude grid of about 1.5 km was used to examine the structure of the global mean vertical ozone profile and its night-time-to-daytime variation in the upper atmosphere. Two distinct ozone layers are predicted, separated by a sharp drop in the ozone concentration near the mesopause. This naturally occurring mesopause ozone deep minimum is primarily produced by the rapid increase in the destruction of water vapour, and hence increase in HO_x, at altitudes between 80 and 85 km, a region where water-vapour photodissociation by ultraviolet radiation of the solar Lyman-alpha line is significant, and where the supply of water vapour is maintained by methane oxidation even for very dry conditions at the tropospheric-stratospheric exchange region. The model indicates that the depth of the mesopause ozone minimum is limited by the efficiency with which inactive molecular hydrogen is produced, either by the conversion of atomic hydrogen to molecular hydrogen via one of the reaction channels of H with HO₂, or by Lyman-alpha photodissociation of water vapour via the channel that leads to the production of molecular hydrogen. The ozone concentration rapidly recovers above 85 km due to the rapid increase in O produced by the photodissociation of O₂ by absorption of ultraviolet solar radiation in the Schumann-Runge bands and continuum. Above 90 km, there is a decrease in ozone due to photolysis as the production of ozone through the three-body recombination of O₂ and O becomes slower with decreasing pressure. The model also predicts two peaks in the night-time/daytime ozone ratio, one near 75 km and the other near 110 km, plus a strong peak in the night-time/daytime ratio of OH near 110 km. Recent observational evidence supports the predictions of the model.     

基于激光雷达手段的海南地区重力波与其波谱的季节分布特性研究

利用子午工程海南激光雷达对我国海南地区上空进行持续观测,通过3年的累积观测数据对我国低纬度地区重力波活动的季节分布特性进行研究,依据重力波线性理论对海南地区上空的大气密度扰动规律、空间功率谱及时间频率谱进行分析,并通过选择波长在1km至8km范围内具有特定波长以及具有波动周期为60 min至25min的特定频率的重力波辅助研究大气密度扰动的季节变化规律,总结得出海南地区重力波活动具有夏季大、春秋季小、而冬季依然频繁的季节性分布规律.结合海南地区特殊的地理位置与当地季节性气候特征分析得出海南地区上空重力波活动季节性变化的可能原因为青藏高原地形及我国南海地区存在的热带强对流与赤道潜流共同作用的结果.     

Gravity wave intensity and momentum fluxes in the mesosphere over Shigaraki, Japan (35°N, 136°E) during 1986–1997

Averaged seasonal variations of wind perturbation intensities and vertical flux of horizontal momentum produced by internal gravity waves (IGWs) with periods 0.2/1 h and 1/6 h are studied at the altitudes 65/80 km using the MU radar measurement data from the middle and upper atmosphere during 1986/1997 at Shigaraki, Japan (35°N, 136°E). IGW intensity has maxima in winter and summer, winter values having substantial interannual variations. Mean wave momentum flux is directed to the west in winter and to the east in summer, opposite to the mean wind in the middle atmosphere. Major IGW momentum fluxes come to the mesosphere over Shigaraki from the Pacific direction in winter and continental Asia in summer.     

A numerical model of the zonal mean circulation of the middle atmosphere

The annual cycle of the zonally averaged circulation in the middle atmosphere (16–96 km) is simulated using a numerical model based on the primitive equations in log pressure coordinates. The circulation is driven radiatively by heating due to solar ultraviolet absorption by ozone and infrared cooling due to carbon dioxide and ozone (parameterized as a Newtonian cooling). Since eddy fluxes due to planetary waves are neglected in the model, the computed mean meridional circulation must be interpreted as thediabatic circulation, not as the total eulerian mean. Rayleigh friction with a short (2–4 day) time constant above 70 km is included to simulate the strong mechanical dissipation which is hypothesized to exist in the vicinity of the mesopause due to turbulence associated with gravity waves and tides near the mesopause.Computed mean winds and temperatures are in general agreement with observations for both equinox and solstice conditions. In particular, the strong mechanical damping specified near the mesopause makes it possible to simulate the cold summer and warm winter mesopause temperatures without generating excessive mean zonal winds. In addition, the model exhibits a strong semiannual cycle in the mean zonal wind at the equator, with both amplitude and vertical structure in agreement with the easterly phase of the observed equatorial semiannual oscillation.Contribution No. 497, Department of Atmospheric Sciences, University of Washington, Seattle.     

Persistent, Widespread, and Strongly Absorbing Haze Over the Himalayan Foothills and the Indo-Gangetic Plains

We examine the impact of the Atmospheric Brown Clouds on the direct radiative forcing of the Himalayan foothills and the Indo-Gangetic Plains (IGP) regions, home for over 500 million S. Asians. The NASA-Terra MODIS satellite data reveal an extensive layer of aerosols covering the entire IGP and Himalayan foothills region with seasonal mean AODs of about 0.4 to 0.5 in the visible wavelengths (0.55 micron), which fall among the largest seasonal mean dry season AODs for the tropics. We show new surface data which reveal the presence of strongly absorbing aerosols that lead to a large reduction in solar radiation fluxes at the surface during the October to May period. The three-year mean (2001 to 2003) October to May seasonal and diurnal average reduction in surface solar radiation for the IGP region is about 32 (±5) W m⁻² (about 10% of TOA insolation or 20% of surface insolation). The forcing efficiency (forcing per unit optical depth) is as large as −27% (note that the forcing is negative) of top-of-atmosphere (TOA) solar insolation, and exceeds the forcing efficiency that has been observed for other polluted regions in America, Africa, East Asia, and Europe. General circulation model sensitivity studies suggest that both the local and remote influence of the aerosol induced radiative forcing is to strengthen the lower atmosphere inversion, stabilize the boundary layer, amplify the climatological tendency for a drier troposphere, and decrease evaporation. These aerosol-induced changes could potentially increase the life times of aerosols, make them more persistent, and decrease their single scattering albedos, thus potentially leading to a detrimental positive feedback between aerosol concentrations, aerosol forcing, and aerosol persistence. In addition, both the model studies and observations of pan evaporation suggest that the reduction in surface solar radiation may have led to a reduction in surface evaporation of moisture. These results suggest the vulnerability of this vital region to air pollution related direct and indirect (through climate changes) impacts on agricultural productivity of the region.     

Enhancement of stratospheric aerosols after solar proton event

The lidar measurements at Verhnetulomski observatory (68.6°N, 31.8°E) at Kola peninsula detected a considerable increase of stratospheric aerosol concentration after the solar proton event of GLE (ground level event) type on the 16/02/84. This increase was located at precisely the same altitude range where the energetic solar protons lost their energy in the atmosphere. The aerosol layer formed precipitated quickly (1–2 km per day) during 18, 19, and 20 February 1984, and the increase of R(H) (backscattering ratio) at 17 km altitude reached 40% on 20/02/84. We present the model calculation of CN (condensation nuclei) altitude distribution on the basis of an ion-nucleation mechanism, taking into account the experimental energy distribution of incident solar protons. The meteorological situation during the event was also investigated.deceased     

Model results for the ionospheric E region: solar and seasonal changes

A new, empirical model for NO densities is developed, to include physically reasonable variations with local time, season, latitude and solar cycle. Model calculations making full allowance for secondary production, and ionising radiations at wavelengths down to 25 Å, then give values for the peak density N _mE that are only 6% below the empirical IRI values for summer conditions at solar minimum. At solar maximum the difference increases to 16%. Solar-cycle changes in the EUVAC radiation model seem insufficient to explain the observed changes in N _mE, with any reasonable modifications to current atmospheric constants. Hinteregger radiations give the correct change, with results that are just 2% below the IRI values throughout the solar cycle, but give too little ionisation in the E-F valley region. To match the observed solar increase in N _mE, the high-flux reference spectrum in the EUVAC model needs an overall increase of about 20% (or 33% if the change is confined to the less well defined radiations at <150 Å). Observed values of N _mE show a seasonal anomaly, at mid-latitudes, with densities about 10% higher in winter than in summer (for a constant solar zenith angle). Composition changes in the MSIS86 atmospheric model produce a summer-to-winter change in N _mE of about–2% in the northern hemisphere, and +3% in the southern hemisphere. Seasonal changes in NO produce an additional increase of about 5% in winter, near solar minimum, to give an overall seasonal anomaly of 8% in the southern hemisphere. Near solar maximum, reported NO densities suggest a much smaller seasonal change that is insufficient to produce any winter increase in N _mE. Other mechanisms, such as the effects of winds or electric fields, seem inadequate to explain the observed change in N _mE. It therefore seems possible that current satellite data may underestimate the mean seasonal variation in NO near solar maximum. A not unreasonable change in the data, to give the same 2:1 variation as at solar minimum, can produce a seasonal anomaly in NmE that accounts for 35–70% of the observed effect at all times.     

The equatorial ionospheric anomaly in electron content from solar minimum to solar maximum for South East Asia

Median hourly, electron content-latitude profiles obtained in South East Asia under solar minimum and maximum conditions have been used to establish seasonal and solar differences in the diurnal variations of the ionospheric equatorial anomaly (EIA). The seasonal changes have been mainly accounted for from a consideration of the daytime meridional wind, affecting the EIA diffusion of ionization from the magnetic equator down the magnetic field lines towards the crests. Depending upon the seasonal location of the subsolar point in relation to the magnetic equator diffusion rates were increased or decreased. This led to crest asymmetries at the solstices with (1) the winter crest enhanced in the morning (increased diffusion rate) and (2) the same crest decaying most rapidly in the late afternoon (faster recombination rate at lower ionospheric levels). Such asymmetries were also observed, to a lesser extent, at the equinoxes since the magnetic equator (located at about 9○N lat) does not coincide with the geographic equator. Another factor affecting the magnitude of a particular electron content crest was the proximity of the subsolar point, since this increased the local ionization production rate. Enhancements of the EIA took place around sunset, mainly during the equinoxes and more frequently at solar maximum, and also there was evidence of apparent EIA crest resurgences around 0300 LST for all seasons at solar maximum. The latter are thought to be associated with the commonly observed, post-midnight, ionization enhancements at midlatitudes, ionization being transported to low latitudes by an equatorward wind. The ratio increases in crest peak electron contents from solar minimum to maximum of 2.7 at the equinoxes, 2.0 at the northern summer solstice and 1.7 at northern winter solstice can be explained, only partly, by increases in the magnitude of the eastward electric field E overhead the magnetic equator affecting the [E×B] vertical drifts. The most important factor is the corresponding increase in ionization production rate due to the increase in solar radiation flux. The EIA crest asymmetries observed at solar maximum were less significant, and this is probably due to the corresponding increase in ionization densities leading to an increase of the retarding effect of ion-drag on the daytime meridional winds.     

Short-period planetary waves in the Antarctic middle atmosphere

Planetary waves with periods between two and four days in the middle atmosphere over Antarctica are characterized using one year of data from the medium-frequency spaced antenna (MFSA) radars at Scott Base, Rothera, and Davis. In order to investigate the origin of the observed waves, the ground-based data are complemented by temperature measurements from the Earth Observing System Microwave Limb Sounder (EOS MLS) instrument on the Aura satellite as well as wind velocity data from the United Kingdom Met. Office (UKMO) stratospheric assimilation. Observed characteristics of waves with a period of approximately two days in summer are consistent with the quasi-two-day wave (QTDW) generally found after the summer solstice at low- and mid-latitudes. The Scott Base observations of the QTDW presented here are the highest-latitude ground-based observations of this wave to date. Waves with preferred periods of two and four days occur in bursts throughout the winter with maximum activity in June, July, and August. The mean of the two- and four-day wave amplitudes is relatively constant, suggesting constant wave forcing. When several waves with different periods occur at the same time, they often have similar phase velocities, supporting suggestions that they are quasi-non-dispersive. In 2005, a “warmpool” lasts from late July to late August. An alternative interpretation of this phenomenon is the presence of a structure propagating with the background wind. Consideration of the role of vertical shear (baroclinic instabilities) and horizontal shear (barotropic instabilities) of the zonal wind suggests that instabilities are likely to play a role in the forcing of the two- and four-day waves, which are near-resonant modes and thus supported by the atmosphere.     

Solar cycle changes to planetary wave propagation and their influence on the middle atmosphere circulation

Recent observations suggest that there may be a causal relationship between solar activity and the strength of the winter Northern Hemisphere circulation in the stratosphere. A three-dimensional model of the atmosphere between 10–140 km was developed to assess the influence of solar minimum and solar maximum conditions on the propagation of planetary waves and the subsequent changes to the circulation of the stratosphere. Ultraviolet heating in the middle atmosphere was kept constant in order to emphasise the importance of non-linear dynamical coupling. A realistic thermo-sphere was achieved by relaxing the upper layers to the MSIS-90 empirical temperature model. In the summer hemisphere, strong radiative damping prevents significant dynamical coupling from taking place. Within the dynamically controlled winter hemisphere, small perturbations are reinforced over long periods of time, resulting in systematic changes to the stratospheric circulation. The winter vortex was significantly weakened during solar maximum and western phase of the quasi-biennial oscillation, in accordance with reported 30 mb geopotential height and total ozone measurements.     

Ionospheric Scintillation and Dynamics of Fresnel-Scale Irregularities in the Inner Region of the Equatorial Ionization Anomaly

This paper reports differences in the occurrence statistics of global positioning system (GPS) L-band scintillations at observational sites located in the inner regions of the northern and southern crests of the equatorial ionization anomaly. Ground-based GPS data acquired at the closed magnetically aligned stations of Manaus (3.1°S; 59.9°W; dip lat. 6.2°N) and Cuiabá (15.5°S; 56.1°W; dip. lat. 6.2°S), Brazil, from December 2001 to February 2007 are used in the analysis. The drift dynamics of Fresnel-scale ionospheric irregularities at the southern station of Cuiabá are also investigated. Only geomagnetically quiet days with the sum of daily Kp < 24 were used in the analysis statistics and in the irregularity drift studies. The results reveal a clear dependence of the scintillation occurrence with the solar activity, but there exists an asymmetry in the percentage of scintillation occurrence between the two stations throughout the period analyzed. The nocturnal occurrence of the scintillations over Cuiabá is predominantly larger than over Manaus, but this scenario seems to change with the decline in the solar activity (mainly during local post-midnight hours). A broad minimum and maximum in the scintillation occurrence appears to occur over both the stations, respectively, during the June solstice (winter) and December solstice (summer) months. The dynamics of the Fresnel-scale irregularities, as investigated from the estimations of the mean zonal drift velocities, reveals that the amplitude of the eastward drifts tends to reduce with the decline in the solar activity. The magnitude of the zonal drift velocities during the December solstice months is larger than during the equinoxes, with the differences being more pronounced at solar maximum years. Other relevant aspects of the observations, with complementary data from a low-latitude ionospheric model, are highlighted and discussed.     

Horizontal eddy diffusivities in the northern hemispheric stratosphere and lower mesosphere

Meterological rocket soundings, launched between 1969–74 at six locations representative of low, middle, and high altitudes, are employed with the use of the statistical theory of diffusion, to determine the zonal and meridional component of eddy diffusivity between 30 and 55 km as a function of season, latitude, and altitude. A comparison is also made between annually-averaged eddy diffusivities above and below 30 km.It is shown that the zonal component of eddy diffusivity is approximately three to five times as large as the meridional component, in most cases. Both components of eddy diffusivity vary greatly with season, latitude, and altitude. Highest eddy diffusivities, found in the vicinity of the winter westerly jet, are approximately one order of magnitude higher than those present during the summer. Tropical eddy diffusivities, however, remain relatively small throughout the year. Annually, a minimum is indicated near 25 km between maximums located at the stratopause and tropopause.