GEO 1009 ANNALS OF GEOPHYSICS, SUPPLEMENT TO VOL. 47, N. 2/3, 2004 5 Long-term trends in the ionosphere and upper atmosphere parameters JÜRGEN BREMER (1), LUCILLA ALFONSI (2), PAL BENCZE (3), JAN LAŠTOVIČKA (4), ANDREI V. MIKHAILOV (5) and NEIL ROGERS (6) (1) Leibniz-Institute of Atmospheric Physics, Kühlungsborn, Germany (2) Istituto Nazionale di Geofisica e Vulcanologia, Roma, Italy (3) Geodetic and Geophysical Research Institute, Hungarian Academy of Sciences, Sopron, Hungary (4) Institute of Atmospheric Physics, Academy of Sciences of Czech Republic, Prague, Czech Republic (5) Institute of Terrestrial Magnetism, Ionosphere and Radio Wave Propagation (IZMIRAN), Russian Academy of Sciences, Troitsk (Moscow Region), Russia (6) Centre for RF Propagation and Atmospheric Research, QinetiQ, Malvern, U.K. The first part of the paper is directed to the investigation of the practical importance of possible long- term trends in the F2-layer for ionospheric prediction models. Using observations of about 50 different ionosonde stations with more than 30 years data series of foF2 and hmF2, trends have been derived with the solar sunspot number R12 as index of the solar activity. The final result of this trend analysis is that the differences between the trends derived from the data of the individual stations are relatively large, the calculated global mean values of the foF2 and hmF2 trends, however, are relatively small. There- fore, these small global trends can be neglected for practical purposes and must not be considered in ionospheric prediction models. This conclusion is in agreement with the results of other investigations analyzing data of globally distributed stations. As shown with the data of the ionosonde station Trom- sø, however, at individual stations the ionospheric trends may be markedly stronger and lead to essen- tial effects in ionospheric radio propagation. The second part of the paper deals with the reasons for pos- sible trends in the Earth’s atmo- and ionosphere as investigated by different methods using characteris- tic parameters of the ionospheric D-, E-, and F-regions. Mainly in the F2-region different analyses have been carried out. The derived trends are mainly discussed in connection with an increasing greenhouse effect or by long-term changes in geomagnetic activity. In the F1-layer the derived mean global trend in foF1 is in good agreement with model predictions of an increasing greenhouse effect. In the E-region the derived trends in foE and h´E are compared with model results of an atmospheric greenhouse effect, or explained by geomagnetic effects or other anthropogenic disturbances. The trend results in the D-re- gion derived from ionospheric reflection height and absorption measurements in the LF, MF and HF ranges can at least partly be explained by an increasing atmospheric greenhouse effect. 5.1. INTRODUCTION Trend investigations of different ionospheric and atmospheric parameters became more important during recent years as such trends could be caused by anthropogenic activities. In this paper we will concentrate on two main questions: 1) Are there trends in the ionosphere which have to be considered in radio wave prediction models? 2) What are the reasons for possible trends in the Earth’s atmosphere and ionosphere? 1010 Jürgen Bremer, Lucilla Alfonsi, Pal Bencze, Jan Laštovička, Andrei V. Mikhailov and Neil Rogers The results related to these two topics are presented in Sections 5.2 and 5.3, discussed in Section 5.4 and briefly summarized in Section 5.5. 5.2. PRACTICAL ASPECTS OF IONOSPHERIC TRENDS In this section some practical aspects of ionospheric radio propagation are discussed. We restrict our investigations to trends in the ionospheric F2-region using their height and maximum electron density data from different ionosonde stations. In Section 5.2.1 mean global trends are derived, whereas in Section 5.2.2 some practical aspects of the ionospheric radio propagation are presented using data from an individual ionosonde station. 5.2.1. Derivation of mean global trends As is well-known from many investigations the behavior of the ionosphere is markedly controlled by solar activity. In most of the ionospheric prediction models the solar activity is described by the solar sunspot number R12. Therefore, in the analyses to investigate the practical aspects of possible ionospheric trends this parameter has been used. As the ionospheric radio propagation in the HF range is mainly dependent on the behavior of the ionospheric F2-region (height and maximum elec- tron density), in the following the results of trend analyses of hmF2 and foF2 observations all around the world are presented. The analyses have been carried out for data series mainly starting at 1957 with a duration of at least 30 years. Altogether the results of 53 stations with foF2 data and 50 sta- tions with hmF2 data have been investigated. The hmF 2 values were derived from M(3000)F2 data using the simple formula of Shimazaki (1955). The trend method used starts with calculations of sim- ple regression equations using monthly mean values the ionospheric parameter X (= hmF2 or foF2) at each full hour in dependence on the solar activity assuming a linear dependence X a b R12th $= + (5.1a) or a quadratic equation X a b cR R12 12 2 th $ $= + + . (5.1b) Then the differences between the observed Xobs and the calculated values Xth are estimated X X X thobs= -∆ (5.2a) or relative differences are calculated X X X Xth thobs= -∆ ^ h . (5.2b) Using mean yearly ∆X values derived by averaging the monthly ∆X values at each full hour (averages of 12 × 24 data series), mean linear tends have been estimated according to X d e year$= +∆ . (5.3) In fig. 5.1a,b two histograms are shown for the trend results of foF2 observations at 53 different ionosonde stations. In the upper part the solar activity has been eliminated by a linear relation ac- cording to eq. (5.1a), in the lower part by the quadratic relation according to eq. (5.1b). In both cas- es differences according to eq. (5.2a) have been used in the trend estimations. The distributions of the 1011 Long-term trends in the ionosphere and upper atmosphere parameters individual foF2 trends are relatively broad, the median trends however are relatively near to zero. As to be expected both data sets are well correlated with a correlation coefficient r = 0.92. In fig. 5.2a,b the trend results of hmF2 observations at 50 different stations are shown in a simi- lar way as for foF2 in fig. 5.1a,b. Also here relatively broad distributions of the individual trends can be seen, whereas the median values are again very near to zero. The correlation between both trend ensembles is very good with r = 0.995. In table 5.I the mean trends of foF2 and hmF2 (median as well as mean values) have been com- piled. These data can partly be found in figs. 5.1a,b and 5.2a,b and from trend results using a linear so- lar activity influence (eq. (5.1a)) and relative differences ∆X (eq. (5.2b)) not shown before. For the mean trends the error values with a significance level of 95% are also presented (for calculation see Taubenheim, 1969). All mean trends are smaller than the corresponding error limits, indicating that the mean trends are not significantly different from zero. Nevertheless from all mean and median values in table 5.I the expected changes for ∆foF2 and ∆hmF2 have been estimated for a time period of 100 years. As to be seen these values are very much smaller compared with changes induced by the vari- able solar activity. Therefore, possible long-term changes have not to be introduced in ionospheric pre- diction models. Until now global mean trends in foF2 and hmF2 have no essential practical meaning. Fig. 5.1a,b. Histograms of 53 individual foF2 trends derived after a linear (a) and quadratic (b) elimination of the solar influence. The corresponding median values are marked by arrows. b a 1012 Jürgen Bremer, Lucilla Alfonsi, Pal Bencze, Jan Laštovička, Andrei V. Mikhailov and Neil Rogers Fig. 5.2a,b. Histograms of 50 individual hmF2 trends derived after a linear (a) and quadratic (b) elimination of the solar influence. The corresponding median values are marked by arrows. Table 5.I. Mean trends of foF2 and hmF2 derived from observations at about 50 ionosonde stations (beginning 1957, more than 30 years) by different methods. ∆foF2 and ∆hmF2 are the expected changes of foF2 and hmF2 during 100 years. Trend (R12) Trend (R12 2 ) Trend (R12) rel foF2 Mean trend + 0.0003 MHz/yr − 0.0013 MHz/yr − 0.016%/yr 95% error ± 0.0027 MHz/yr ± 0.0023 MHz/yr ± 0.041%/yr ∆∆foF2mean ++ 0.03 MHz −− 0.13 MHz −− 0.13 MHz Median trend + 0.0009 MHz/yr − 0.0010 MHz/yr + 0.001%/yr ∆∆foF2med [MHz] ++ 0.09 MHz −− 0.10 MHz ++ 0.008 MHz hmF2 Mean trend + 0.058 km/yr + 0.063 km/yr + 0.021%/yr 95% error ± 0.096 km/yr ± 0.092 km/yr ± 0.030%/yr ∆∆hmF2mean ++ 5.8 km ++ 6.3 km ++ 6.3 km Median trend + 0.012 km/yr + 0.009 km/yr + 0.004%/yr ∆∆hmF2med [km] ++ 1.2 km ++ 0.9 km ++ 1.2 km a b 1013 Long-term trends in the ionosphere and upper atmosphere parameters 5.2.2. Trends at an individual station The ionosonde at Tromsø in Northern Norway (69°N, 19°E) (Hall and Hansen, 2003) is a unique resource providing continuous observations of ionospheric critical frequencies since 1935 and layer heights since 1952. Hall and Cannon (2001, 2002) have examined long-term trends in F-region parameters at Tromsø using a method similar to that developed by Bremer (1992) and used in Section 5.3.1 below, eq. (5.4). Using daily ionosonde measurements taken at or near noon, linear components of the background variations in R12 and geomagnetic index (Ap) are subtract- ed to leave a residual trend. Based on a recently extended set of ionograms, the peak height of the F2 layer, hmF2 is found to fall at − 0.106 ± 0.087 km per year and its critical frequency, foF2 falls at − 0.013 ± 0.002 MHz per year, where the errors indicate one standard deviation. Ulich et al. (2003) presented an hmF2 data set for a nearby ionosonde at Sodankylä, Finland for the period 1958-2003 with a trend of − 0.41 ± 0.04 km per year – a significantly steeper trend than that seen in the Tromsø data. Modeling of ionospheric radio system performance requires a propagation tool that contains a three-dimensional map of ionospheric electron density. This may be used to find the skip zones of High Frequency (HF) radio transmissions, or the group delay, phase advance, etc., in trans- ionospheric systems. Parameterized «monthly median» models are often employed, and these use worldwide maps of foF2 and hmF2 and similar parameters, often referred to as CCIR coefficients (CCIR, 1966; Jones et al., 1969) or URSI coefficients (Rush et al., 1989). The former were based on ionosonde data for the years 1954-1958 and 1964 only and foF2 was described by a set of spe- cial mapping functions. The worldwide RMS error in foF2 is ∼ 0.5 MHz. The URSI coefficients are derived from a much larger database of 45 000 station-months and include data from 1975- 1979. URSI coefficients showed improved accuracy over the oceans but at the expense of accu- racy over land, so that overall the CCIR and URSI maps are of similar accuracy (Bilitza, 2002). Fig. 5.3. Impact of 16 km reduction in hmF2 on ground ranges for propagation due south of Tromsø at 8 MHz for December, 12 UT, R12 = 50. 1014 Jürgen Bremer, Lucilla Alfonsi, Pal Bencze, Jan Laštovička, Andrei V. Mikhailov and Neil Rogers Five decades have elapsed since the data were collected to generate the CCIR coefficients, and in that time the F2-region over Tromsø has fallen by 5.3 km. Further, foF2 has decreased by 0.65 MHz, an amount already bigger than the RMS model error of the CCIR model. The impact on HF systems has been assessed by ray tracing through two sets of analytical Dudeney profiles of elec- tron density (Dudeney, 1983). The first set of profiles was fit to the CCIR coefficients and the sec- ond was modified to reduce the peak height by 16 km. Based on the analyses of Tromsø iono- grams a 16 km drop in hmF2 is appropriate to the conditions in the year 2100 (provided there is no change in the long-term trends). If we had chosen to use the results of Ulich et al. (2003) then a 16 km height reduction was appropriate to the year 1990! For a 16 km drop it was found that ground range errors of ∼ 100 km per hop will be common if system designers and operators con- tinue to employ the old CCIR maps. A sample of the ray tracing results is presented in fig. 5.3, representing 8 MHz 1-hop propagation south of Tromsø near noon on a December day. Overall, the ground range errors (lower curve) are 50-150 km for the 1F2 mode of propagation (above 21° elevation), with greater errors at elevations just above that at which the E-layer is penetrated. These results demonstrate that for trends derived at individual stations their influence on ionospheric radio propagation could be much stronger than for global mean trends. This result could be expected from the broad trend histograms shown in figs. 5.1a,b and 5.2a,b. 5.3. SCIENTIFIC ASPECTS OF IONOSPHERIC TRENDS This section presents the trend results which were derived by different European groups using partly different methods. The results are separately discussed for the different ionospheric re- gions. 5.3.1. Trends in the F2-region Using a similar analysis method as described in Section 5.2.1 Bremer (1998, 2001) investigated data series of more than 100 ionosonde stations. For the elimination of the solar and geomagnetical- ly induced parts, however, the following twofold regression equation was used: X a b R c Apth $ $= + + (5.4) with the solar sunspot number R and the geomagnetic activity index Ap. Also other solar indices have been tested (F10.7, E10.7, R12), but the derived trends are not markedly influenced by this choice (Bre- mer, 2001). In fig. 5.4 the updated histograms with the individual trends of foF2 and hmF2 are pre- sented together with the corresponding median values. As in the investigations discussed in Section 5.2.1, the distributions of the individual trends for both parameters are relatively broad, and the me- dian values are very small. As shown in table 5.II the mean trends of foF2 and hmF2 are not statisti- cally significant different from zero as follow from the corresponding mean errors. In fig. 5.5 the individual hmF2-trends are presented in dependence on longitude and latitude. Nega- tive trends are marked by blue, positive trends by red symbols. Significant trends with a 95% relia- bility level are characterised by full dots, non significant trends by circles. There is no clear depend- ence of the trends on latitude or longitude, but there seem to be some regional differences. Whereas in the western part of Europe most trends are negative, in Central Asia the trends are more positive. The causes of these regional differences are not known. Similar trend analyses have been made by Ulich and Turunen (1997a) at first with hmF2 data of Sodankylä and later for more than 50 stations (Ulich and Turunen, 1997b). These authors also found similar regional differences in the hmF2 trends as shown in fig. 5.5. Ulich (2000) also investigated 1015 Long-term trends in the ionosphere and upper atmosphere parameters Fig. 5.4. Histograms for foF2 and hmF2 trends derived from ionosonde observation at different stations. Num- ber of stations are given in brackets. Table 5.II. Mean experimental trends and error limits (95%) derived from trend analyses of different ionosonde parameters at different stations. Significant mean trends are bold-faced. Parameter N Mean trend Error (95%) F2-region foF2 106 − 0.0018 MHz/yr ± 0.0025 MHz/yr hmF2 87 − 0.009 km/yr ± 0.076 km/yr F1-region foF1 51 0.0027 MHz/yr ±± 0.0011 MHz/yr E-region foE 72 0.0014 MHz/yr ±± 0.0007 MHz/yr h′E 31 − 0.040 km/yr ± 0.070 km/yr Fig. 5.5. Trends of hmF2 in dependence on latitude and longitude derived from data sets of 87 different ionosonde stations. Negative trends are blue, positive are red, significant trends are marked by full dots, non sig- nificant trends by circles. 1016 Jürgen Bremer, Lucilla Alfonsi, Pal Bencze, Jan Laštovička, Andrei V. Mikhailov and Neil Rogers different formulas to derive hmF2 from ionosonde M(3000)F2 data. He found quantitative differ- ences in the derived trends depending on the hmF2 formula used. This finding is caused by the fact that the ionisation below the F2-layer is differently considered by the different formulas. Therefore, trends in the underlying ionisation parameters my also cause differences in the hmF2 trends. Bencze et al. (1998) derived trends in hmF2 for some Asian and American stations using the same method as Bremer (2001) and Ulich and Turunen (1997a,b). He also found similar regional differ- ences with more negative trends at coastal regions and more positive trends inside the continents. Alfonsi et al. (2001, 2002a) derived foF2 trends for different stations at high and middle latitudes using monthly mean values as well as hourly values. In their analyses only such experimental data have been introduced which are not accompanied by qualification letters to reduce the uncertainty of the hourly foF2 values. The investigation of the ionospheric long-term changes follows an analytical approach similar to that shown in eq. (5.4). The original contribution comes from the use of Ap(τ) and MACAp(τ), used for investigating the geomagnetic control on the ionospheric trends. Ap(τ) and MACAp(τ) have been introduced for representing the non-linear relationship between geomagnetic activity and ionospheric behaviour (Wrenn, 1987; Perrone and De Franceschi, 1999). In particular, Ap(τ) is a time-weighted accumulation of Ap index, and MACAp(τ) is the related catalogue, intro- duced for characterising four levels of geomagnetic impact on the ionosphere ranging from 0 (quiet magnetic condition) to 3 (severely disturbed magnetic condition) (De Franceschi et al., 1999, 2002). Extracting from long data series of hourly foF2 data those values recorded under magnetically quiet conditions (MACAp(τ) = 0), Alfonsi et al. (2002a,b) found a long-term decrease of the F2 plasma fre- quency. All foF2 trends derived from data of four different stations shown in fig. 5.6 are negative with a median value of − 0.0035 MHz/yr. A totally different method for derivation of trends has been developed by Mikhailov and Marin (2000, 2001) and Marin et al. (2001). They use 12-monthly running hourly mean values of foF2 or hmF2 and of the solar sunspot number (R12) and geomagnetic activity (Ap12). For foF2 only data during three years near solar minimum and maximum are taken into consideration to avoid the hysteresis effect whereas for hmF2 data of all years have been analysed. Data gaps are filled by the MQMF2 model (Mikhailov et al., 1996). The regression model used is a modification of eq. (5.1b) according to X a b R c R d R e Apth 12 12 2 12 2 12$ $ $ $= + + + + ^ h. (5.5) The analyses have been made with or without the geomagnetic term. Relative deviations of the ob- served data have been used according to eq. (5.2b). The main results can be summarized: 1) The influence of the geomagnetic activity can only partly be removed also if Ap12 is introduced in eq. (5.5). 2) Periods with increasing and decreasing geomagnetic activity are related to different ionos- pheric trends: foF2 trends are in anti-phase, hmF2 trends in-phase with Ap trends. This behavior can be seen in fig. 5.7 for long-term observations at the station Slough. 3) The diurnal and latitudinal variations of the foF2 and hmF2 trends can be explained by the cur- rent geomagnetic storm behavior. Therefore, this method is called geomagnetic control concept by their authors. Danilov (2002, 2003) developed a method to eliminate the long-term changes in geomagnetic ac- tivity in foF2 trends. Starting with relative ∆foF2 values, calculated in a similar way as by Mikhailov and Marin (2000) but using data of all available years, he assumes that the derived foF2-trend k(obs) is the result of a linear combination of a non-geomagnetic trend k(tr) and a geomagnetically caused trend. He further assumes that the geomagnetic activity is described by annual mean values of Ap and the geomagnetic trend is proportional to the gradient of Ap, k(Ap) ( ) ( ) ( ).k k a k Aptr obs 1$= + (5.6) 1017 Long-term trends in the ionosphere and upper atmosphere parameters Here al is a scaling coefficient which includes the efficiency of the geomagnetic activity impact on foF2 and may therefore be different for different stations and local times. The details of the method to estimate the al values can be found in Danilov (2002, 2003). Using foF2 data of 23 stations with observations of more than 35 years a mean non-geomagnetic trend of − 0.012 MHz/yr was found for the period between 1958 and the mid-nineties and − 0.0075 MHz/yr for an earlier interval between 1948-1985 (here only data of 8 stations were available). Another attempt to eliminate the geomagnetic influence on the foF2 trends has been tried by Mikhai- lov et al. (2002). They start their analysis with monthly regressions according to eq. (5.1a), estimate rela- tive deviations after eq. (5.2b) and calculate 11-year running mean values from these deviations: δfoF2132. The geomagnetic activity effect in δfoF2132 is derived from 11-year running mean values Ap132 by ( ) ( )foF b b Ap t n b Ap t n2 th132 0 1 132 2 132 2 $ $= + + + +d (5.7) where n is the time shift in years of Ap132 with respect to δfoF2132 variations. From the residuals foF foF foF2 2 2obs th132 132= -d d d (5.8) linear trends according to eq. (5.3) have been calculated with the following main results: 1) The derived foF2 trends are nearly independent of latitude and of the phase of long-term changes in geomagnetic activity (basic point of geomagnetic control concept of Mikhailov and Marin, 2000, 2001). Fig. 5.6. Long-term trends of foF2 from different stations under magnetically quiet condition (MACAp(τ) = 0) and after removal of the solar cycle influence. 1018 Jürgen Bremer, Lucilla Alfonsi, Pal Bencze, Jan Laštovička, Andrei V. Mikhailov and Neil Rogers 2) Most of the residual trends are negative, often however not significant. In fig. 5.8 a significant negative trend is shown using the long-term data series observed at Slough. Summarizing it can be stated that in all analyses with data of different stations the derived mean trends in the F2-region are relatively small and not important for practical applications. At individual stations (Ulich and Turunen, 1997a; Rogers et al., 2002), however, the trends can be markedly more pronounced, resulting in broad global histograms as shown in fig. 5.4. Fig. 5.7. (Top) Long-term variation of geomagnetic activity (11-year running mean of Ap) together with rela- tive deviations of foF2 (middle) and of hmF2 (bottom) derived from ionosonde observations at Slough. 1019 Long-term trends in the ionosphere and upper atmosphere parameters Fig. 5.8. Relative deviations of foF2 after very long-term ionosonde observations at the station Slough (full line: linear trend; Kr: trend parameter). 5.3.2. Trends in the F1-region Trends in foF1 data series were derived by Bremer (2001) using observations at 51 different ionosonde stations. For the elimination of the solar and geomagnetically induced parts the twofold re- gression eq. (5.4) was used. In fig. 5.9 the histogram of the updated individual trends is shown together with the global median value (marked by an arrow). The distribution is asymmetric to zero as slightly seen by essentially more positive (grey) than negative trends (white). The estimated mean trend with 0.0027 MHz/yr is significantly different from zero as seen in table 5.II where the corresponding mean er- ror can also be found. Therefore this trend has been marked there by bold letters. 5.3.3. Trends in the E-region Trends in the ionospheric E-region were estimated by Bremer (2001) from h´E and foE observa- tions at different ionosonde stations. For the elimination of the solar and geomagnetically induced parts the twofold regression eq. (5.4) has been used. The results of the trend analyses are shown by histograms together with their global median values (marked by arrows) in fig. 5.10. The estimated mean foE trend is significant with a reliability level of more than 95%, whereas the mean h´E trend is not significant different from zero. This can be seen in table 5.II where the mean values and their mean errors are presented. Mikhailov and de la Morena (2003) investigated foE trends from different ionosonde stations with their revised trend method which was developed for foF2-trend investigations (Mikhailov et al., 2002). This method is briefly described above. Instead of eq. (5.8), they used a slightly modified equation ( )foE b b Ap t nth132 0 1 132$= + +d b (5.9) 1020 Jürgen Bremer, Lucilla Alfonsi, Pal Bencze, Jan Laštovička, Andrei V. Mikhailov and Neil Rogers where β is a fitting parameter and n again the time shift in years of Ap132 with respect to δfoE132. Their main results are: 1) Before about 1970 δfoE132 and Ap132 are negatively correlated. The authors call it «natural foE variation». 2) After about 1970 positive foE trends were detected which could anthropogenically be caused. In fig. 5.11 an example is shown of the dependence between δfoE132 and Ap132 using ionosonde data observed at the station Rome. 5.3.4. Trends in the D-region In the ionospheric D-region the available data series for trend analyses are limited, therefore we restricted these investigations to some reflection height and ionospheric absorption data observations. Fig. 5.9. Histogram for foF1 trends derived from ionosonde observations at different stations. Number of sta- tions is given in brackets. Fig. 5.10. Histogram for foE and h´E trends derived from ionosonde observations at different stations. Num- bers of stations are given in brackets. 1021 Long-term trends in the ionosphere and upper atmosphere parameters From the interference field between the ionospherically reflected sky wave and the ground wave of a far distant radio transmitter it is possible to derive ionospheric reflection heights (often called phase heights, for details of the method see Bremer and Berger, 2002, and references therein). From long-term field strength observations of a measuring path at 162 kHz in mid-latitudes trends of the reflection height have been derived for a constant solar zenith angle χ = 78.4°. In fig. 5.12 the yearly trend can be seen. Here the solar and geomagnetically induced variations have been removed ac- cording to eq. (5.4) using the solar Lyman α radiation instead of the solar sunspot number R. A sig- nificant negative trend could be derived with about − 0.028 km/yr. Such a shrinking of the lower ionosphere could also be detected from observations of reflection heights at ionospheric drift measurements during night-time (Kürschner and Jacobi, 2002). These measurements were made in the LF range (177 kHz) at short distance between transmitter and re- ceiver (170 km). Both mean trends of reflection heights near 82 km (LF phase heights) and near 95 km (drift meas- urements) as well as the trends in the E-layer peak height near 100-110 km (h´E-trends in Section 5.3.3) have negative values which may slightly increase with increasing height. The significance lev- el of the phase height trend is high (mostly higher than 99%), for the two other methods this level is however lower due to the shorter data series (drift method) or due to marked differences between in- dividual stations (ionosonde data, see the broad histogram in fig. 5.10 and the derived mean trend and error values in table 5.II). Ionospheric absorption data have been analysed for three different reflection height ranges, below about 80 km, between about 80-90 km, and above about 90 km (Bremer and Laštovička, 2002; Laštovička and Bremer, 2004). From long-term observations of the absorption of a very long measuring path at 164 kHz (distance between transmitter and receiver: 1720 km) during summer and noon conditions (Nestorov et al., 1991, and updated later by Laštovička and Pancheva, 1999) a clear positive trend could be detected. The reflection height of this measuring path is clearly below 80 km. Fig. 5.11. Variation of δfoE132 in dependence on Ap132 for foE data series observed at the station Rome with the marked increase after 1972. 1022 Jürgen Bremer, Lucilla Alfonsi, Pal Bencze, Jan Laštovička, Andrei V. Mikhailov and Neil Rogers Using, however, data from LF absorption measurements at short distances (about 180-220 km) the derived trends change their general behavior. In fig. 5.13 the absorption trends are shown for three different frequencies at the same solar zenith angle χ = 90°. Whereas the trend at the lowest frequen- cy is negative, the trend increases with increasing frequency and becomes even positive at the high- est frequency. That means that with increasing reflection height the trends become more positive. We observe the same behavior if trends have been estimated for the same measuring path but at different solar zenith angles (example can be found in Bremer, 2003). With increasing solar zenith angle the trends becomes again more positive. Investigating ionospheric absorption measurements in the MF and HF ranges at oblique incidence (distances between transmitter and receiver are about 170 km up to 550 km) only positive trends have Fig. 5.12. Long-term variation of yearly mean phase-heights at mid-latitudes (162 kHz, χ = 78.4°) after elimi- nation of the solar and geomagnetically induced variations. Fig. 5.13. Trends of ionospheric absorption for three measuring paths with different frequencies at a constant solar zenith angle χ = 90° after elimination of the solar and geomagnetically induced parts. 1023 Long-term trends in the ionosphere and upper atmosphere parameters been estimated. Some examples can be seen in Laštovička (2001, 2002) and Bremer and Laštovička (2002). The reflection heights of these measuring paths are above 90 km. The tendency detected at the LF absorption data series that with increasing reflection height the trends become positive, is therefore confirmed by the absorption measurements in the MF and HF ranges. 5.4. DISCUSSION Two main problems will be discussed here concerning a) the practical importance of ionospheric trends and b) the scientific reasons for possible trends in the different regions of the ionosphere. 5.4.1. Practical aspects Ionospheric HF propagation is mainly influenced by the behaviour of the F2-region expressed by their characteristic parameters foF2 and hmF2. As shown in Section 5.2.1 and summarised in table 5.I, possible global trends of both parameters are very small. This finding is in agreement with other investigations presented in Section 5.3.1 where the results of more complex trend analyses are shown using observations from all over the world. Therefore, in global ionospheric prediction models it is not necessary to include the effect of ionospheric long-term changes. The solar and geomagnetically induced variations are markedly more important. At individual stations, however, the trends can be essentially more pronounced. This can be seen in the relatively broad histograms of the foF2-trends in fig. 5.1a,b and of hmF2-trends in fig. 5.2a,b. The trend values of the stations Tromsø and Sodankylä used in Section 5.2.2 to investigate their impact on ionospheric radio propagation are well inside the corresponding histograms. Un- fortunately we do not know the reasons for the strong differences of the trends derived at individ- ual stations. One reason for such differences between the trends at (partly neighbouring) stations could be artificial steps in the data series caused by technical changes at the ionosonde equipment and/or changes of the evaluation algorithms. In Bremer (1998, 2001) some examples of such dis- turbances can be found. Therefore, the quality of ionosonde data sets have to be carefully con- trolled before trend analyses are carried out (Ulich et al., 2003). As long as we do not have final decisions about the quality of the ionosonde data sets we can only use mean or median values of all reasonable data sets, hoping that some errors of the derived individual trends may compensate each other. 5.4.2. Scientific aspects What are the reasons of trends in the different layers of the ionosphere? At least two candidates are often discussed: anthropogenic pollution (e.g., CO2, CH4, O3, H2O, ...) and geomagnetic trends. To investigate the influence of an increasing atmospheric greenhouse effect the trends in the ionospheric F2-, F1-, and E-regions can be compared with model calculations of Rishbeth (1990) and of Rishbeth and Roble (1992). Their theoretical results have been derived for a doubling of the at- mospheric greenhouse gases CO2 and CH4. The effective change in greenhouse gases during the last 40 years where trends of the ionosonde data have mainly been investigated is about 20% (Brasseur and de Rudder, 1987; Houghton et al., 2001). Assuming a linear dependence between the content of the atmospheric greenhouse gases and the ionospheric effect, the experimental trends of table 5.II can be extrapolated to a level of doubled greenhouse gases. These values called CO2*2 (exp) are com- pared with the corresponding model values CO2*2 (mod) of Rishbeth (1990) and Rishbeth and Roble (1992) in table 5.III. 1024 Jürgen Bremer, Lucilla Alfonsi, Pal Bencze, Jan Laštovička, Andrei V. Mikhailov and Neil Rogers Table 5.III. Mean experimental (exp) trends of different ionospheric parameters (from table 5.II) and expect- ed changes of these data assuming a doubling of the atmospheric greenhouse gases (CO2* 2). The model data (mod) are from Rishbeth (1990) and Rishbeth and Roble (1992). Parameter Mean exp trend CO2* 2 (exp) CO2* 2 (mod) F2-region foF2 − 0.0018 MHz/yr − 0.36 MHz − 0.2 ... − 0.5 MHz hmF2 − 0.009 km/yr − 1.8 km − 10 ... − 20 km F1-region foF1 0.0027 MHz/yr 0.54 MHz 0.3 ... 0.5 MHz E-region foE 0.0014 MHz/yr 0.28 MHz 0.05 ... 0.08 MHz h′E − 0.040 km/yr − 8.0 km − 2.5 km 5.4.2.1. F2-region The agreement between the experimental and model trends in foF2 and hmF2 seems to be quite reasonable looking at the data in table 5.III. However, the differences between the individual trends in foF2 and hmF2 at different stations are very strong, and the derived mean trends are not statistically significant different from zero as can be seen in table 5.II. Therefore, the agreement between model and experimental data is accidental. The reason for the strong variability in the F2-region is not quite clear. After fig. 5.5 there seem to be regional differences which could be caused by dynamical effects in the plasma of the F2-region. Such strong regional differences of the hmF2 trends were also detect- ed earlier by other authors using a more limited data volume (Ulich and Turunen, 1997b; Bencze at al., 1998; Bremer, 1998). From satellite observations (Keating et al., 2000) it is known that the ob- served long-term neutral density reduction near 350 km altitude is in good agreement with model cal- culations of an increasing greenhouse effect (Akmaev, 2002). That means that the possible greenhouse effect in ionospheric data series is superposed by unknown dynamical processes which are more im- portant for the variability of the ionized component than for the neutral gas at F2-region heights. Using differences of hmF2-h´F Bencze (2002) derived long-term trends which contain informa- tion about the plasma temperature of the lower part of the F2-region. For most of the used 20 Euro- pean stations a negative trend was detected in qualitative agreement with an expected cooling of the atmosphere due to an increasing atmospheric greenhouse effect. With the method developed by Mikhailov and Marin (2000, 2001) the estimated trends in foF2 and hmF2 are explained by the remaining influence of geomagnetic activity (geomagnetic control concept). After their investigations their derived trends depending on local time and latitude can be explained by long-term variations of geomagnetic activity. With a revised method Mikhailov et al. (2002) tried to remove this geomagnetic influence. The re- maining relative trend in foF2 noon values for the station Slough with − 2.2 10−4 year−1 can be convert- ed to a mean trend of − 0.0022 MHz/yr if a mean foF2 value of 10 MHz is assumed. This trend value agrees well with the mean global trend for foF2 in table 5.II and is thus in reasonable agreement with the model value in table 5.III. Therefore, this foF2 trend can be explained by an increasing greenhouse effect. However, an influence of the very long-lasting increase of the geomagnetic activity, as found in long-term variations of the geomagnetic aa index during the last century (Clilverd et al., 1998; Gulyae- va, 2002), on trends in foF2 should also lead to a negative foF2-trend due to the negative correlation be- tween trends in geomagnetic activity and foF2 after the geomagnetic control concept. Alfonsi et al. (2002a,b) tried to remove the influence of geomagnetic activity on foF2-trends by using only geomagnetically quiet hourly data in their analyses. Their mean foF2-trend derived from 4 different stations (see fig. 5.6) with − 0.0035 MHz/yr is slightly stronger but not very different from the expectations after the greenhouse effect. 1025 Long-term trends in the ionosphere and upper atmosphere parameters Danilov (2002, 2003) used quite another method to derive non-geomagnetic trends in foF2. His mean trend, derived from 23 different stations, is with − 0.012 MHz/yr also negative but markedly stronger than the trends in table 5.II as well as those from Mikhailov et al. (2002) and Alfonsi et al. (2002a,b). But such differences are not very surprising if we remember the strong variability between the trends detected at different stations and the fact that these authors used different stations in their analyses. 5.4.2.2. F1-region The agreement of the mean experimental trend in the F1-region (described in Section 5.3.2) with model trends in table 5.III is surprisingly good. Therefore, the global foF1-trend can be explained by an increasing atmospheric greenhouse effect. 5.4.2.3. E-region As seen from table 5.III, in the E-region the experimental and theoretical trend values agree qual- itatively with a lowering of the height h´E and an increase in foE. However the experimental trends are markedly stronger than the model values. The derived positive foE trend is also in general agree- ment with rocket mass spectrometer measurements of the ion density ratio [NO+]/[O2+] in the E-re- gion (Danilov and Smirnova, 1997). The observed negative trends of [NO+]/[O2+] cause increasing electron densities and therefore increasing foE values as the dissociative recombination coefficient of NO+ is markedly larger than that of O2+. For years after 1970 Mikhailov and de la Morena (2003) also found with their revised trend method increasing foE values, a quantitative comparison with the model trends in table 5.III is how- ever impossible. Nevertheless they propose an increasing chemical pollution due to an increasing number of rocket launchings and/or an increasing atmospheric greenhouse effect to explain the ob- served foE increase. 5.4.2.4. D-region The shrinking of the lower ionosphere as derived from LF reflection height observations (fig. 5.12) can be explained by a cooling of the strato- and mesosphere as expected from an in- creasing atmospheric greenhouse effect (Bremer and Berger, 2002). The estimated amount of the atmospheric cooling depends however also on possible changes in the density of nitric ox- ide nNO and of the effective recombination coefficient αeff near the reflection height of about 82 km. Unfortunately there are no experimental data of long-term changes in these two quantities, however there exist model results which predict a decrease of nNO (Roble and Dickinson, 1989; Beig, 2000) and an increase in αeff in the mesosphere (Chakrabarty, 1997; Beig, 2000). As shown by Bremer and Laštovička (2002) and Bremer (2003) a common interpretation of the trends of the LF reflection heights and of the absorption data at different frequencies is only possible if the trends in nNO and αeff as predicted from model calculations are taken into consideration. The mesospheric temperature trends derived from the trends in LF reflection heights and ionospheric absorption data are in reasonable agreement with other experimental trends derived from long- term lidar and rocket measurements (Bremer, 2003). As shown by model results (Bremer and Berger, 2002; Bremer, 2003) the experimental temperature trends can qualitatively be explained by an increasing atmospheric greenhouse effect. The experimental trends are however stronger than the model results. 1026 Jürgen Bremer, Lucilla Alfonsi, Pal Bencze, Jan Laštovička, Andrei V. Mikhailov and Neil Rogers 5.5. SUMMARY AND CONCLUSIONS Using data of reflection height observations in the LF range and ionospheric absorption meas- urements in the LF, MF, and HF ranges for trend analyses in the D-region as well as data of ionosonde observations for trend estimations in the E-, F1-, and F2-regions, the following short summary and conclusions can be given: General aspects – Long-term ionospheric data series have to be carefully checked concerning their homogene- ity. Discontinuities caused by different technical changes can markedly disturb the results of trend analyses. – The solar and geomagnetically induced variations of most ionospheric parameters are essen- tially stronger than the long-term trends and have carefully to be eliminated in trend analyses. F2-region – The mean global trends are unimportant for practical ionospheric prediction models. In such models the solar and geomagnetic variability are the most important external factors. – Due to a wide variability of the individual trends in the F2-region no significant global trends could be derived for foF2 and hmF2. Therefore, the relatively reasonable agreement between the mean global experimental and model results could be accidental. The regional differences of the trends hint at an unknown dynamical process which superposes a possible greenhouse effect in the F2-region. – Using different hmF2-h´F however a general cooling of the lower part of the F2-region could be derived in qualitative agreement with the greenhouse effect. – The derived trends after the geomagnetic control concept contain geomagnetically induced components causing typical local time and latitudinal trend variations. – Different approaches to eliminate geomagnetic activity changes lead to negative foF2-trends as proposed by model results of an atmospheric greenhouse effect. The experimental trends are however slightly stronger than the model predictions. If the geomagnetic influence could not be totally removed, then these negative trends could also be explained by the very long-term increase in geomagnetic activity. F1-region – The mean global trend in the F1-region (increase in foF1) agrees quite well with model re- sults of an increasing greenhouse effect. E-region – The mean global trends in the E-region (lowering of h´E, increase of foE) are in qualitative agreement with other experimental results and with model results of an increasing greenhouse ef- fect. However, the experimental trends are stronger than the model results. – Investigations with a revised trend method lead to an foE-increase after about 1970. D-region – The trends in the D-region (lowering of reflection heights and different trends in ionospher- ic absorption in dependence on frequency of the measuring path) can qualitatively be explained by an increasing atmospheric greenhouse effect. The experimental trends are however stronger than expected from model results. Summarizing, the investigation of ionospheric trends is an important task for the understanding of structure and dynamics of the ionosphere in general and for the derivation of anthropogenic influ- 1027 Long-term trends in the ionosphere and upper atmosphere parameters ences in particular. Combined experimental trend analyses (including quality control of existing da- ta series as well as introduction of additional data sets) and model activities are necessary for future trend investigations. ACKNOWLEDGEMENTS The investigations presented here were mainly carried out in the COST 271 project «Effects of the upper atmosphere on terrestrial and Earth-space communications (EACOS)». The authors are very grateful to all participants who took actively part in the analysis of trends in the Earth’s atmos- phere and ionosphere: P.S. Cannon, U.K.; G. De Franceschi, Italy; T. Gulyaeva, Russia; C.M. Hall, Norway; M. Materassi, Italy; B. de la Morena, Spain, L. Perrone, Italy; A. Poole, South Africa; G. Sole, Spain; Th. Ulich, Finland. The authors thank R. Conkright of the NGDC, Boulder, U.S.A., the U.K. ionosonde group, the World Data Centre C-1, the IPS (Australia) and IRF (Sweden) for providing ionosonde data. REFERENCES AKMAEV, R.A. (2002): Modeling of the cooling due to CO2 increases in the mesosphere and lower thermosphere, Phys. Chem. Earth, 27, 521-528. ALFONSI, L., G. DE FRANCESCHI and L. PERRONE (2001): Long term trend in the high latitude iono- sphere, Phys. Chem. Earth, 26, 303-307. ALFONSI, L., G. DE FRANCESCHI, L. PERRONE and M. MATERASSI (2002a): Long-term trends of the critical frequency of the F2-layer at northern and southern high latitude regions, Phys. Chem. Earth, 27, 607-612. ALFONSI L., G. DE FRANCESCHI, L. PERRONE and C. SCOTTO (2002b): Long-term trends of the iono- sphere at mid and high latitude regions, in Proceedings URSI, August 2002, Maastricht, The Netherlands. BEIG, G. (2000): The relative importance of solar activity and anthropogenic influences on the ion composition, temperature, and associated neutrals of the middle atmosphere, J. Geophys. Res., 105, 19841-19856. BENCZE, P. (2002): Some results referring to the long-term change of ionospheric parameters, Acta Geod. Geophys. Hung., 37 (4), 403-408. BENCZE, P., G. SOLE, L.F. ALBERCA and A. POOLE (1998): Long-term changes of hmF2: possible lat- itudinal and regional variations, in Proceedings of the 2nd COST 251 Workshop (RAL U.K.), 30- 31 March 1998, Side, Turkey, 107-113. BILITZA, D. (2002): Ionospheric models for radio propagation studies, Review of Radio Sciences, ed- ited by W.R. STONE (IEEE and Wiley), 625-679. BRASSEUR, G. and A. DE RUDDER (1987): The potential impact on atmospheric ozone and tempera- ture of increasing trace gas concentrations, J. Geophys. Res., 92, 10903-10920. BREMER, J. (1992): Ionospheric trends in mid-latitudes as a possible indicator of the atmospheric greenhouse effect, J. Atmos. Terr. Phys., 54, 1505-1511. BREMER, J. (1998): Trends in the ionospheric E- and F-regions over Europe, Ann. Geophysicae, 16, 986-996. BREMER, J. (2001): Trends in the thermosphere derived from global ionosonde observations, Adv. Space Res., 28 (7), 997-1006. BREMER, J. (2003): Mesospheric temperature trends at mid-latitudes derived from LF radio propaga- tion experiments, Adv. Space Res., 32 (9), 1653-1662. BREMER, J. and U. BERGER (2002): Mesospheric temperature trends derived from ground-based LF 1028 Jürgen Bremer, Lucilla Alfonsi, Pal Bencze, Jan Laštovička, Andrei V. Mikhailov and Neil Rogers phase-height observations at mid-latitudes: comparison with model simulations, J. Atmos. Solar- Terr. Phys., 64, 805-816. BREMER, J. and J. LAŠTOVIČKA (2002): Trends in the ionospheric D- and E-regions using results of dif- ferent radio propagation experiments, in Proceedings of the 2nd COST 271 Workshop, 2-4 October 2002, Faro, Portugal. CCIR (1966): CCIR Atlas of Ionospheric Characteristics (CCIR, Geneva). CHAKRABARTY, D.K. (1997): Mesopause scenario on doubling of CO2, Adv. Space Res., 20 (11), 2117-2125. CLILVERD, M.A., T.D.G. CLARK, E. CLARKE and H. RISHBETH (1998): Increased magnetic storm ac- tivity from 1868 to 1995, J. Atmos. Solar-Terr. Phys., 60, 1047-1056. DANILOV, A.D. (2002): The method of determination of the long-term trends in the F2-region inde- pendent of geomagnetic activity, Ann. Geophysicae, 20, 511-521. DANILOV, A.D. (2003): Long-term trends of foF2 independent of geomagnetic activity, Ann. Geo- physicae, 21, 1167-1176. DANILOV, A.D. and N.V. SMIRNOVA (1997): Long-term trends in the ion composition of the E-region, Geomagn. Aeron., 37 (4), 35-40 (in Russian). DE FRANCESCHI, G., T.L. GULYAEVA, L. PERRONE and B. ZOLESI (1999): MAC: an oriented Magnetic Activity Catalogue for ionospheric applications, URSI IRI Newsl., 6 (4), 5-6. DE FRANCESCHI, G., T.L. GULYAEVA, L. PERRONE and B. ZOLESI (2002): A statistical analysis of ionos- pheric irregularities at mid and high latitudes, Inverse Problems, 18, 67-78. DUDENEY, J.R. (1983): An improved model of electron concentration with height in the ionosphere, J. Atmos. Terr. Phys., 45, 629-640. GULYAEVA, T. (2002): Forecast of recurrent magnetic storms with a one-day lead time, Geomagn. Aeron., 42, 159-164 (in Russian). HALL, C.M. and P.S. CANNON (2001): Indication of the shrinking atmosphere above Tromsø (69°N, 19°E), Atmos. Sci. Lett., R. Meteorol. Soc., doi10.1006/asle.2001.0036. HALL, C.M. and P.S. CANNON (2002): Trends in foF2 above Tromsø (69°N-19°E), Geophys. Res. Lett., 29 (23), 2128, doi:10.1029/2002GL016259. HALL, C.M. and T.L. HANSEN (2003): 20th century operation of the Tromsø ionosonde, Adv. Polar Upper Atmos. Res., 17, 155-166. HOUGHTON, J.T., Y. DING, D.J. GROGGS, M. NOGUER, P.J. VAN DER LINDEN, X. DAI, K. MASKELL and C.A. JOHNSON (2001): Climate change: the scientific basis, in Contribution of WG 1 to the 3rd As- sessment Report of the IPCC (Cambridge University Press). JONES, W.B., R.P. GRAHAM and M. LEFTIN (1969): Advances in ionospheric mapping by numerical methods, ESSA, US Government Printing Office, Washington DC. KEATING, G.M., R.H. TOLSON and M.S. BRADFORD (2000): Evidence of long term global decline in the Earth’s thermospheric densities apparently related to anthropogenic effects, Geophys. Res. Lett., 27, 1523-1526. KÜRSCHNER, D. and CH. JACOBI (2002): Trends and periodicities in nighttime LF radio wave reflec- tion heights, Paper presented at the 27th General Assembly EGS, Symposium ST11, Nice. LAŠTOVIČKA, J. (2001): Long-term trends in the lower ionosphere, Adv. Space Res., 28 (7), 1007- 1016. LAŠTOVIČKA, J. (2002): Long-term changes and trends in the lower ionosphere, Phys. Chem. Earth, 27, 497-507. LAŠTOVIČKA, J. and J. BREMER (2004): An overview of long-term trends in the lower ionosphere be- low 120 km, Surv. Geophys., 25, 69-99. LAŠTOVIČKA, J. and D. PANCHEVA (1999): Trends in the characteristics of the annual and semiannual variations observed in the radio wave absorption in the lower ionosphere, Ann. Geophysicae, 17, 1336-1343. MARIN, D., A.V. MIKHAILOV, B.A. DE LA MORENA and M. HERRAIZ (2001): Long-term hmF2 trends 1029 Long-term trends in the ionosphere and upper atmosphere parameters in the Eurasian longitudinal sector from ground-based ionosonde observations, Ann. Geophysi- cae, 19, 761-772. MIKHAILOV, A.V. and B.A. DE LA MORENA (2003): Long-term trends of foE and geomagnetic activity variations, Ann. Geophysicae, 21, 751-760. MIKHAILOV, A.V. and D. MARIN (2000): Geomagnetic control of the foF2 long-term trends, Ann. Geo- physicae, 18, 653-665. MIKHAILOV, A.V. and D. MARIN (2001): An interpretation of the foF2 and hmF2 long-term trends in the framework of the geomagnetic control concept, Ann. Geophysicae, 19, 733-748. MIKHAILOV, A.V., V.V. MIKHAILOV and M.G. SKOBLIN (1996): Monthly median foF2 and M(3000)F2 ionospheric model over Europe, Ann. Geofis., XXXIX (4), 791-805. MIKHAILOV, A.V., D. MARIN, T.YU. LESCHINSKAYA and M. HERRAIZ (2002): A revised approach to the foF2 long-term trend analysis, Ann. Geophysicae, 20, 1663-1675. NESTOROV, G., D. PANCHEVA and A.D. DANILOV (1991): Climatic changes of ionospheric absorption of radio waves in the LF range, Geomagn. Aeron., 31, 1070-1073 (in Russian). PERRONE, L. and G. DE FRANCESCHI (1999): A correlation study between time-weighted geomagnet- ic indices and high latitude ionosphere, Phys. Chem. Earth, 24, 389-392. RISHBETH, H. (1990): A greenhouse effect in the ionosphere?, Planet. Space Sci., 38, 945-948. RISHBETH, H. and R.G. ROBLE (1992): Cooling of the upper atmosphere by enhanced greenhouse gas- es – Modelling of the thermospheric and ionospheric effects, Planet. Space Sci., 40, 1011-1026. ROBLE, R.G. and R.E. DICKINSON (1989): How will changes of carbon dioxide and methane modify the mean structure of the mesosphere and thermosphere?, Geophys. Res. Lett., 16, 1441-1444. ROGERS, N., P. CANNON and C. HALL (2002): Practical implications of long term trends in the critical frequency and the height of the F2-layer above Tromsø (70°N, 19°E), in Proceedings of the 2nd COST 271 Workshop, 2-4 October 2002, Faro, Portugal. RUSH, C., M. FOX, D. BILITZA, K. DAVIES, L. F. MCNAMARA, F. STEWART and M. POKEMPNER (1989): Ionospheric Mapping – An update of foF2 coefficients, Telecom J., 56, 179-182. SHIMAZAKI, T. (1955): World wide daily variations in the height of the maximum electron density in the ionospheric F2-layer, J. Radio Res. Labs., Japan, 2, 85-97. TAUBENHEIM, J. (1969): Statistische Auswertung geophysikalischer und meteorologischer Daten, Akad. Verlagsgesellschaft Geest und Portig K.-G., Leipzig. ULICH, TH. (2000): Solar variability and long-term trends in the atmosphere, Sodankylä Geophysical Observatory Publications No. 87, Oulu. ULICH, TH. and E. TURUNEN (1997a): Evidence for long-term cooling of the upper atmosphere in ionospheric data, Geophys. Res. Lett., 24, 1103-1106. ULICH, TH. and E. TURUNEN (1997b): Long-term behaviour of ionospheric F2-layer peak height on a global scale, paper presented at Session 2.18 of the 8th Scientific Assembly of IAGA, Uppsala. ULICH, TH., M.A. CLILVERD and H. RISHBETH (2003): Determining the long term change in the iono- sphere, EOS Trans. Am. Geophys. Un., 84 (52), 581-585. WRENN, G.L. (1987): Time-weighted accumulations Ap(τ) and Kp(τ), J. Geophys. Res., 92, 125-129.