Trends in the F2 ionospheric layer due to long-term variations in the Earth's magnetic field

The Earth's magnetic field presents long-term variations with changes in strength and orientation. Particularly, changes in the dip angle (I) and, consequently, in the sin(I)cos(I) factor, affect the thermospheric neutral winds that move the conducting plasma of the ionosphere. In this way, a lowering or lifting of the F2-peak (hmF2) is induced together with changes in foF2, depending on season, time and location. A simple theoretical approximation, developed in a previous work, is extended to a worldwide latitude–longitude grid to assess hmF2 and foF2 trends due to Earth's magnetic field secular variations. Compared to the greenhouse gases effects over the ionosphere, the Earth's magnetic field may be able to produce stronger trends which vary with season, time and location. However, to elucidate the origin of F2-region trends, long-term variations in the three possible known mechanisms should be considered altogether—greenhouse gases, geomagnetic activity and Earth's magnetic field.     

Ionospheric long-term trends for South American mid-latitudes

Using ionosonde made observations at Concepción (36.8°S; 73.0°W) for the 1958–1994 interval, long-term trends of critical frequency (foF2) and peak height (hmF2) of the ionospheric F2-layer are analysed. The trends found for different times-of-day and all seasons are consistent with an increasing diurnal-variation amplitude of both foF2 and hmF2. An increasing hmF2 trend of up to 1.5 km/year found between midnight and dawn during winter has no precedent. It is suggested that these long-term amplitude changes may be associated with changes in the prevailing thermospheric meridional neutral winds.     

The zonal E×B plasma drift effects on the low latitude ionosphere electron density at solar minimum near equinox

The F2-layer peak density, NmF2, and peak altitude, hmF2, which were observed by 12 ionospheric sounders during the 20 September 1964 geomagnetically quiet time period at solar minimum are compared with those calculated by the three-dimensional time-dependent theoretical model of the Earth's low and middle latitude ionosphere and plasmasphere. The modeled NmF2 are also compared with those measured during the geomagnetically quiet time periods of 12–15, 18–21, and 26 September 1964 to take into account observed day-to-day ionospheric variability. Major features of the data are reproduced by the model if the corrected HWM90 neutral wind is used. The changes in NmF2 due to the zonal E×B plasma drift are found to be less than 20% in the daytime low latitude ionosphere. The model, which does not take into account the zonal E×B plasma drift, underestimates night-time NmF2 up to the maximum factor of 2 at low geomagnetic latitudes. The night-time increase of NmF2 caused by the zonal E×B plasma drift is less pronounced at −20° and 20° geomagnetic latitudes in comparison with that between −10° and 10° geomagnetic latitude. The longitude dependence of the calculated night-time low latitude influence of the zonal E×B plasma drift on NmF2 is explained in terms of the longitudinal asymmetry in B (the eccentric magnetic dipole is displaced from the Earth's center and the Earth's eccentric tilted magnetic dipole moment is inclined with respect to the Earth's rotational axis), and the variations of the wind induced plasma drift and the meridional E×B plasma drift in geomagnetic longitude. The difference between the hmF2 values calculated by including the effect of zonal E×B drift and that obtained when it is excluded does not exceed 19 km in the low latitude ionosphere. Over the geomagnetic equator the zonal E×B plasma drift produces the maximum increase in the electron density by a factor of 1.06–1.48 and 1.05–1.30 at 700 and 1000 km altitude, respectively, and this increase is not significant above about 1500 km. Changes in the vertical electron content, VEC, caused by the zonal E×B plasma do not exceed 16% during the day, while the value of the night-time VEC is increased up to a factor of 1.4 due to this drift. The maximum effects of the zonal E×B plasma drift on the night-time electron density derived from the model results corresponding to solar minimum and maximum are quite comparable.     

The ionospheric F2-region at low geomagnetic latitudes during the geomagnetic storms of 22–26 April 1990: Comparison of observed and modeled response

A comparison between the modeled NmF2 and hmF2 and NmF2 and hmF2, which were observed by the Kokubunji, Okinawa, Manila, Vanimo, and Darwin ionospheric sounders and by the middle and upper (MU) atmosphere radar, have been used to study the time-dependent response of the low-latitude ionosphere to geomagnetic forcing during a time series of geomagnetic storms from 22 to 26 April 1990. The reasonable agreement between the model results and data requires the modified equatorial meridional E×B plasma drift, the modified HWM90 wind, and the modified NRLMSISE-00 neutral densities. We found that changes in a flux of plasma into the nighttime equatorial F2-region from higher L-shells to lower L-shells caused by the meridional component of the E×B plasma drift lead to enhancements in NmF2 close to the geomagnetic equator. The equatorward wind-induced plasma drift along magnetic field lines, which cross the Earth equatorward of about 20° geomagnetic latitude in the northern hemisphere and about −19° geomagnetic latitude in the southern hemisphere, contributes to the maintenance of the F2-layer close to the geomagnetic equator. The nighttime weakening of the equatorial zonal electric field (in comparison with that produced by the empirical model of Fejer and Scherliess [Fejer, B.G., Scherliess, L., 1997. Empirical models of storm time equatorial zonal electric fields. J. Geophys. Res. 102, 24047–24056] or Scherliess and Fejer [Scherliess, L., Fejer, B.G., 1999. Radar and satellite global equatorial F region vertical drift model. J. Geophys. Res. 104, 6829–6842) in combination with corrected equatorward nighttime wind-induced plasma drift along magnetic field lines in the both geomagnetic hemispheres are found to be the physical mechanism of the nighttime NmF2 enhancement formation close to the geomagnetic equator over Manila during 22–26 April 1990. The model crest-to-trough ratios of the equatorial anomaly are used to study the relative role of the main mechanisms of the equatorial anomaly suppression for the 22–26 April 1990 geomagnetic storms. During the most part of the studied time period, a total contribution from geomagnetic storm disturbances in the neutral temperature and densities to the equatorial anomaly changes is less than that from meridional neutral winds and variations in the E×B plasma drift. It is shown that the latitudinal positions of the crests are determined by the E×B drift velocity and the neutral wind velocity.     

Variations in the altitude of the F2 peak associated with trough-formation processes

A novel approach is described which can help to determine, from ground-based data, which of the possible production mechanisms for the mid-latitude F-region ionospheric trough is dominant during a particular event. This approach involves numerically modelling the possible causal mechanisms of the mid-latitude trough to see how each will affect the altitude of the F2-layer electron-concentration peak (hmF2), and then comparing these predictions with the observed variation of hmF2 during trough formation. The modelling work predicts that, if the neutral-wind velocity does not vary, hmF2 will remain almost constant if the trough is formed via stagnation, but will rise if it is formed as a result of high ion velocities or neutral upwelling. Observations made at Halley (76°S, 27°W, L=4.2), Antarctica, show that most frequently the only changes in hmF2 during trough formation are those expected due to variations in the neutral wind, which suggests that stagnation is the most common production mechanism. During the most geomagnetically active night studied, on which Ap varied between 18 and 32, there was a rise in hmF2 that cannot be explained by changes in the neutral wind. On this night the plasma also decayed faster, and the poleward edge of the trough was seen earlier than on other nights. These differences, together with the fact that the ion velocities remained relatively low, suggest the trough was caused by a change in neutral composition, possibly advected into the observing area.     

Critical frequency foF2 as an indicator of trends in thermospheric dynamics

Variations with time during recent decades of three parameters are considered. R(foF2) is the correlation coefficient between the nighttime and daytime values of foF2 within the same day. Stable trends are found for minimal (R(foF2)(min)) and maximal (R(foF2)(max)) values of R(foF2) over the year. The foF2(day)/foF2(night) ratio demonstrates both negative and positive trends; the sign of the trend being governed by the inclination I and declination D of the magnetic field. The correlation coefficient r(h,fo) between foF2 and the 100-hPa level in the stratosphere demonstrates a decrease (both, for the years of maximum and minimum solar activity) from the 1980s to the 1990s. The trends in all three groups of data are considered in the scope of an assumption that there is a long-term change in the circulation in the upper atmosphere. The data considered in the paper provide an indirect confirmation of the existence of this change and show the possibility that further studies of the thermospheric dynamics can be undertaken using ground-based ionospheric observations.     

Global pattern of trends in the upper atmosphere and ionosphere: Recent progress

The global pattern of long-term trends and changes in the upper atmosphere and ionosphere has been presented by Laštovička et al. [2006a. Global change in the upper atmosphere. Science 314 (5803), 1253–1254]. Trends in the mesospheric temperature, electron concentration in the lower ionosphere, electron concentration and height of its maximum in the E-region, electron concentration in the F1-region maximum, thermospheric neutral density and F-region ion temperature qualitatively agree with consequences of the enhanced greenhouse effect and form a consistent pattern of global change in the upper atmosphere. Three groups of parameters were identified as not-fitting this global pattern, the F2-region ionosphere, mesospheric water vapour, and the mesosphere/upper thermosphere dynamics. The paper reports progress in development of the global pattern of trends with emphasis to these three open problems. There are several other factors contributing to long-term trends, namely the stratospheric ozone depletion, mesospheric water vapour concentration changes, long-term changes of geomagnetic activity and of the Earth's magnetic field.     

Comparisons of observed ionospheric F2 peak parameters with IRI-2001 predictions over South Africa

The monthly means of the ionospheric F2 peak parameters (foF2 and hmF2) over three stations in South Africa (Grahamstown, 33.3°S, 26.5°E, Madimbo, 22.4°S, 26.5°E, and Louisvale, 28.5°S, 21.2°E) were analyzed and compared with IRI-2001, using CCIR (Comité Consultatif International des Radio communications) and URSI (Union Radio-Scientifique Internationale coefficients) options. The analysis covers a few selected quiet and disturbed days during various seasons represented by the months of January, April, July and October 2003. IRI-2001 generally overestimates hmF2 for both quiet and disturbed days and it overestimates and underestimates foF2 at different times for all the stations. In general, foF2 is predicted more accurately by IRI-2001 than hmF2, and on average, the CCIR option performed better than the URSI option when predicting both foF2 and hmF2.In general, the model generates good results, although some improvements are still necessary to be implemented in order to obtain better predictions. There are no significant differences in the model predictions of hmF2 and foF2 for quiet and disturbed days.     

An inspection of the long-term behaviour of the range of the daily geomagnetic field variation from comprehensive modelling

This paper attempts to reveal whether long-term trends in the ionosphere are reflected in the amplitude range of the geomagnetic daily variation recorded at ground level. The smooth and regular variation observed in the magnetograms on magnetically quiet days is induced by the ionospheric currents flowing in the dynamo region. So it is likely that trends in the conductivity or in the dynamics of this region could produce changes in the current densities, and consequently in the range of the geomagnetic variation. The crucial aspect is how to separate the changes produced by the geomagnetic activity itself, or by secular changes of the Earth's magnetic field, from the part of the variation produced by factors affecting trends in the ionosphere, which could have an anthropogenic origin. To investigate this, we synthesized for several geomagnetic observatories the daily ranges of the geomagnetic field components with a comprehensive model of the quiet-time, near-Earth magnetic field, and finally we removed the synthetic values from the observed ranges at those observatories. This comprehensive model accounts for contributions from Earth's core, lithosphere, ionosphere, magnetosphere and coupling currents, and, additionally, accounts for influences of main field and solar activity variations on the ionosphere. Therefore, any trend remaining in the residuals, assuming that all the contributions mentioned above are properly described and thus removed by the comprehensive model, should reflect the influence of other sources. Results, based on series of magnetic data from observatories worldwide distributed, are presented. Trends in the X and Z components are misleading, since the current system changes in form as well as in intensity, producing changes of the focus latitude in the course of a solar cycle and from one cycle to another. Some differences exist between the long-term trends in the Y component between the real and modelled ranges, suggesting that other non-direct solar causes to the amplitude changes of the solar quiet geomagnetic variation should not be ruled out. Nevertheless, the results also reflect some short-comings in the way that the comprehensive modelling accounts for the influence of the solar activity on the range of the daily geomagnetic variation.     

On the geomagnetic and ionospheric responses to an intense storm associated with weak IMF Bz and high solar wind dynamic pressure

A study of the geomagnetic storm of July 13–14, 1982, and its ionospheric response is presented using the low-latitude magnetic index, Dst, and interpreted using solar wind interplanetary data: proton number density, solar wind flow speed, interplanetary magnetic field southward component B _Z , and solar wind dynamic pressure. The F2 region structure response to the geomagnetic storm was studied using foF2 data obtained during the storm from a network of various ionosonde stations. Our results appear to show simultaneous abrupt depletion of foF2 that occurred at all latitudes in both the East Asian and African/European longitudinal zone during the period: 18:00–19:00 UT on July 13 and is as result of an abrupt increase in the dynamic pressure between 16:00 and 17:00 UT. The dynamic pressure increased from 3.21 to 28.07 nPa within an hour. The aforementioned abrupt depletion of foF2 simultaneously resulted in an intense negative storm with peak depletion of foF2 at about 19:00 at all the stations in the East Asian longitudinal zone. In the African/European longitudinal zone, this simultaneous abrupt depletion of foF2 resulted in intense negative storm that occurred simultaneously at the low latitude stations with peak depletion at about 20:00 UT on July 13, while the resulting negative storm at the mid latitude stations recorded peak depletion of foF2 simultaneously at about 2:00 UT on July 14. The present results indicate that most of the stations in the three longitudinal zones showed some level of simultaneity in the depletion of foF2 between 18:00 UT on July 13 and 2:00 UT on July 14. The depletion of foF2 during the main phase of the storm was especially strongly dependent on the solar wind dynamic pressure.     

Errors in Estimating of the F2-Layer Peak Parameters in Automatic Systems for Processing the Ionograms in the Vertical Radio Sounding of the Ionosphere under Low Solar Activity Conditions

Geomagnetism and Aeronomy - The article analyzes the errors in estimating the parameters of the main ionospheric maximum, plasma frequency foF2 and its height hmF2, by automated systems for...     

Structural changes in trend parameters of the MLT winds based on wind measurements at Obninsk (55°N, 37°E) and Collm (52°N, 15°E)

Recent analysis of the long-term behavior of different geophysical data has demonstrated that trend parameters can change during a period of observation. Sophisticated general methods for an objective analysis of structural changes in linear trends have been developed during the last 10 years. Such methods are applied for an analysis of changes in trend parameters of the mesosphere/lower thermosphere wind observed over Obninsk (55°N, 37°E) from 1964 to 2007 and Collm (52°N, 15°E) from 1979 to 2008, respectively. We found that trend models with breakpoints are generally preferred against straight lines. At Obninsk, there are break-years in trends of the winter prevailing winds close to 1977, when a climatic regime shift was observed. The break-years in trends of the semidiurnal tides for both stations are close to years of possible changes in stratospheric ozone. Correlations of the Obninsk and Collm winds with atmospheric indices are also considered.     

A comparison of results derived from scaling VS chirp-ionosonde ionograms with the international reference ionosphere (IRI)

The results derived from processing vertical-incidence ionograms obtained with the chirp-ionosonde at Irkutsk for different winter time intervals (February) and at equinox are presented. The peak height hmF2 was determined by Dudeney's formula based on ionogram parameters, including the coefficient M(3000). The algorithm is suggested for determining the coefficient M(3000) in the automatic mode using the conventional form of the transfer curve method without invoking a standard transparency called the “transfer curve”. The parameters foF2 and hmF2 are compared with the international reference ionosphere (IRI-95) model. It is found that in most cases the values of the foF2 and hmF2 parameters, calculated in the IRI-95 model, are similar to the median ones. It is confirmed that for practical purposes where it is necessary to know the radio wave propagation conditions along the propagation path, the IRI model is convenient and attractive.     

On the long-term change of ionospheric parameters

Independent of the possible sources (solar activity, geomagnetic activity, greenhouse effect, etc.) of a global change in the upper atmosphere, it is the sign of a long-term trend of temperature that might reveal the cause of a global change.Long-term change of temperature in the F region of the ionosphere has been studied and is assumed to be expressed in terms of thickness of the bottomside F2 layer characterized by the difference between height of the maximum electron density of the F2 layer hmF2 and altitude of the lower boundary of the F region represented by h′F. Using the difference of two ionospheric parameters has the advantage that it reduces the effect of changes resulting from alteration of equipment and scaling personnel. In this study, in summer only night values of the difference hmF2−h′F and in winter both day and night values have been taken into account considering that h′F might indicate the lower boundary of the F region in these periods. The study of the behaviour of hmF2−h′F taking separately the stations and determining yearly the mean measure (trend) of the variation of hmF2−h′F with solar and geomagnetic activities found that this difference increases significantly with enhanced solar activity, but trends of the solar activity effect exerted on this difference themselves do not practically change with increasing sunspot number. Further, hmF2−h′F decreases only insignificantly with growing geomagnetic activity. Trends of the geomagnetic activity effect related to hmF2−h′F change only insignificantly with increasing Ap; however, trends of the geomagnetic activity effect decreased with increasing latitude.As a result of this investigation it has been found that hmF2−h′F regarded as thickness of the bottomside F2 layer shows an effect of the change of solar activity during the last three solar cycles, indicating temperature change in the upper atmosphere to be expected on the basis of changing solar activity. Furthermore, though a long-term variation of solar activity considering only years around solar activity minima is relatively small, the difference hmF2−h′F indicates a trend opposing the change of solar activity; that is, it decreases slightly during the first three 20, 21, 22 solar cycle minima (1964–1986), but decreases more abruptly according to the change of solar activity towards the minimum of solar cycle 23 (1986–1996), thus also indicating variation of temperature in the F region. However, this variation cannot be explained by the change of solar and geomagnetic activities alone, but assumes some other source (e.g. greenhouse gases) too.     

Upper ionosphere variability over Alma-Ata and Observatorio Del Ebro using the ΔfoF2 data obtained during the winter/spring period of 2003–2004

The foF2 data obtained at Alma-Ata and Observatorio Del Ebro during the winter/spring of 2003–2004 are analyzed to compare and investigate the upper ionosphere variability at the two selected sites. The geomagnetic activity and the middle stratosphere dynamics, involving planetary wave (PW) activity, are analyzed for understanding the physical conditions and processes that can explain the observed ionospheric variability. By applying the same method of wavelet analysis to the data sets and doing a direct comparison of the results, two types of foF2 disturbances were found. The first type is 2–7-day oscillations, which appeared during periods of increased geomagnetic activity. The second type is oscillations arising from PW activity in the lower atmosphere. These consist of (1) 6–11-day oscillations arising from PW activity in lower atmospheric regions developed during the final stratosphere warming and indicating the timing of the transition from the winter to the summer circulation and (2) 9–13-day and 8–10-day oscillations mostly during the quiet level of geomagnetic activity, indicating a likely close relation with those in the geopotential height at the 1 hPa level for westward-propagating waves at 40°N, which strengthened during stratosphere warming events in January 2004. The time delay of the oscillations in the ΔfoF2 with respect to those in the geopotential height is about 10 days and it seems that the assumed ionosphere response can occur under weakened eastward zonal wind or relatively weak westward zonal wind (V<30 m s⁻¹).     

A modelling study of the latitudinal variations in the nighttime plasma temperatures of the equatorial topside ionosphere during northern winter at solar maximum

Latitudinal variations in the nighttime plasma temperatures of the equatorial topside ionosphere during northern winter at solar maximum have been examined by using values modelled by SUPIM (Sheffield University Plasmasphere Ionosphere Model) and observations made by the DMSP F10 satellite at 21.00 LT near 800 km altitude. The modelled values confirm that the crests observed near 15° latitude in the winter hemisphere are due to adiabatic heating and the troughs observed near the magnetic equator are due to adiabatic cooling as plasma is transported along the magnetic field lines from the summer hemisphere to the winter hemisphere. The modelled values also confirm that the interhemispheric plasma transport needed to produce the required adiabatic heating/cooling can be induced by F-region neutral winds. It is shown that the longitudinal variations in the observed troughs and crests arise mainly from the longitudinal variations in the magnetic meridional wind. At longitudes where the magnetic declination angle is positive the eastward geographic zonal wind combines with the northward (summer hemisphere to winter hemisphere) geographic meridional wind to enhance the northward magnetic meridional wind. This leads to deeper troughs and enhanced crests. At longitudes where the magnetic declination angle is negative the eastward geographic zonal wind opposes the northward geographic meridional wind and the trough depth and crest values are reduced. The characteristic features of the troughs and crests depend, in a complicated manner, on the field-aligned flow of plasma, thermal conduction, and inter-gas heat transfer. At the latitudes of the troughs/crests, the low/high plasma temperatures lead to increased/decreased plasma concentrations.     

Variations in the Intensity of the Earth's Magnetic Field

The Earth's main magnetic field can be approximated by an axial, geocentric dipole. The remaining non-dipole field is much smaller and is a regional rather than a global feature – quite large changes can occur in a few ka. This review is concerned with changes in the dipole component of the geomagnetic field, and one of the problems is in separating the non-dipole from the dipole contributions to the field. Unlike the many determinations of the direction of the Earth's magnetic field in the past (which have led to fundamental contributions to our understanding of plate tectonics and shown that the field can on occasion reverse its polarity), estimates of the intensity of the field are comparatively few, especially before the Holocene. This is mainly the result of experimental difficulties in obtaining reliable measurements of the field. These problems are discussed in some detail and are followed by a short account of archaeomagnetic intensities and results from Hawaii where many of the first determinations were obtained. Measurements for the last 100 ka from both lavas and lacustrine and oceanic sediments are reviewed and results from different areas compared. An asymmetric saw-tooth pattern has been observed in some of the records over the last few Ma, and this rather controversial question is discussed. Finally an account is given of the far more limited data on palaeointensities in earlier times.A short discussion is given of the interpretation of coherent short wavelength variations which are observed in many marine magnetic profiles. Although short reversals of the field may be responsible for some of these tiny wiggles, it is more likely that in general they are the result of changes in the strength of the Earth's magnetic field.     

Quasi-biennial oscillation in foF2 at the south crest of the equatorial anomaly

The aim of this paper is to report some periodicities observed in the ionospheric parameter foF2 measured at Tucuman (26.9°S; 65.4°W), station placed near the southern crest of the equatorial anomaly. For that, monthly medians of foF2 at several hours of LT for the period 1958–1987 are used. The data are run with Fast Fourier Transform (FFT). Data gaps (4–5 months) are filled by means of linear interpolation. Several periodicities are present. Besides the solar cycle dominant dependence (11 years), semi-annual, annual, five years and quasi-biennial periodicities are also observed. A marked quasi-biennial periodicity is observed at daytime and nighttime hours being their greater amplitude at local noon and midnight. Different mechanisms or combined effects possibly cause them. It is suggested that the solar activity by means of extreme ultraviolet radiation (EUV), which present a quasi-biennial oscillation (QBO) and it is responsible for the ionization, could be the dominant mechanism for the diurnal quasi-biennial periodicity of foF2. At night, since the photoionization by extreme ultraviolet radiation is not significant and the F2 layer is lower than during daytime (100 km) other mechanism may be operative for the quasi-biennial periodicity observed. Possibly the stratospheric QBO contributes to the modulation of the observed behaviour in foF2 at night. This result is preliminary because it needs to be extended to other stations so as to extract definite conclusions. Moreover, we cannot dismiss the possibility of a combined effect of both these mechanisms mainly at daytime and/or QBO influence of geomagnetic parameters.     

Scatter of hmF2 values as an indicator of trends in thermospheric dynamics

Analysis of the variability of the @F2-layer height, @hmF2@, from the 1950s to 1960s to the 1990s is performed on the basis of vertical sounding data for a series of midlatitude ionospheric stations. It is found that the scatter of the @hmF2@ values (standard deviation) increases sharply from the earlier decades to the later ones. This increase is better pronounced in the spring period of the year and does not depend on geomagnetic activity. The increase in the scatter of @hmF2@ indicates, evidently, the presence of systematic changes (trends) in thermospheric dynamics, the existence of these trends having been suggested in the recent publications of the author on the basis of the analysis of the foF2(night)/foF2(day) and R(foF2) parameters.     

Longitudinal variation in E- and F-region ionospheric trends

A novel technique is used to examine northern hemisphere midlatitude longitudinal variations in ionospheric long-term trends. Differences in hour-by-hour monthly median ionospheric parameters between equilatitudinal observatory pairs are analysed for long-term trends, thus eliminating at source the large solar cycle and geomagnetic variability that normally hinders ionospheric trend calculations. The results confirm the finding of Bremer [1998. Trends in the ionsopheric E- and F-regions over Europe. Annales Geophysicae 16, 698–996] that there are longitudinal variations in the F-region altitude trend across Europe, but suggest the influence of a stationary wave-like feature between 3°W and 104°E. Possible causes such as scaling errors, insufficient account of changes in ionisation underlying the F-region, varying gravity wave fluxes, and secular change in the geomagnetic field are ruled out. The data suggest that the longitudinal variation may reflect long-term changes in a large-scale stationary feature induced via non-migrating tides induced by latent heat release in the troposphere.Significant differences in the long-term trend of E-region peak plasma frequency between observatories were also found. These E-region differential trends varied with solar zenith angle reaching over 0.3 MHz per decade between Juliusruh and Moscow at midday in summer.